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Recent results obtained by applying the method of self-consistent Green's functions to nuclei 
and nuclear matter are reviewed. Particular attention is given to the description of experimental 
data obtained from the (e,e'p) and (e,e'2N) reactions that determine one and two-nucleon removal 
^ I probabilities in nuclei since the corresponding amplitudes are directly related to the imaginary 

Ch ■ parts of the single-particle and two-particle propagators. For this reason and the fact that these 

amplitudes can now be calculated with the inclusion of all the relevant physical processes, it is 
useful to explore the efficacy of the method of self-consistent Green's functions in describing these 



X 



' experimental data. Results for both finite nuclei and nuclear matter are discussed with particular 



emphasis on clarifying the role of short-range correlations in determining various experimental 
quantities. The important role of long-range correlations in determining the structure of low- 
energy correlations is also documented. For a complete understanding of nuclear phenomena it 
is therefore essential to include both types of physical correlations. We demonstrate that recent 
experimental results for these reactions combined with the reported theoretical calculations yield 
a very clear understanding of the properties of all protons in the nucleus. We propose that this 
knowledge of the properties of constituent fermions in a correlated many-body system is a unique 
feature of nuclear physics. 
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7 Summary and Outlook 



1 Introduction 

Various approaches exist to deal with the nuclear many-body problem. A comparison of these ap- 
proaches has recently been given in Ref. pj. In this review only one specific method is singled out for 
discussion. This method is most intimately connected to quantum field theory since it corresponds to 
the nonrelativistic implementation of this approach to the many-body problem. While the propagator 
or Green's function method was the most important tool in the formal development of many-body 
theory |2| - IHI, only in the last ten years has it been applied to many-body problems beyond its mean- 
field implementation (the Hartree-Fock approximation). Examples of such applications have appeared 
in atomic [7| and condensed matter physics [5^. In the atomic case, the essential new development in- 
volves taking the self-consistent determination of the electron propagator from the Hartree-Fock level to 
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the fully self-consistent inclusion of the second-order self-energy. Important quantitative improvements 
over the Hartree-Fock results are obtained for standard quantities like the binding energy but other 
experimental data like spectroscopic factors are significantly improved as well [7j. In the case of the 
electron gas, the self-consistent determination of the electron propagator involves the summation of all 
ring diagrams in the self-energy which is necessary to treat the well known divergence of the Coulomb 
interaction in an infinite system. Also in this case, one obtains a much improved description of the 
energy per particle essentially in agreement with the benchmark Monte Carlo results of Ref. 

In the nuclear case, similar developments have taken place which have focused on the self-consistent 
treatment of short-range correlations (SRC) in nuclear matter by including ladder diagrams in the 
self-energy. For finite nuclei, the self-consistent treatment of long-range correlations (LRC) has been 
considered at different levels of sophistication. The nuclear case is particularly unique due to the 
availability of a magnificent set of experimental data which have illuminated in great detail the properties 
of the nucleon propagator in the medium. These experiments involve nucleon knockout reactions induced 
by electrons which are measured in coincidence with the ejected proton ^Ul - CHI- These experiments 
identify in an essentially unambiguous manner the removal probability of protons near the Fermi energy 
from a variety of target nuclei. The results of these experiments have clearly identified the need to 
consider nucleon properties beyond the mean-field since these removal probabilities have consistently 
been about 35% below the simple shell-model value of 1 for closed-shell systems. By identifying these 
removal probabilities experimentally, one gains crucial insight into those quantities which determine the 
proton propagator in the medium. It is then natural to approach the theoretical description of these 
data by employing the propagator method. 

It is the purpose of this review to document the recent development of the theoretical understanding 
of these (e,e'p) data over the last ten years. The wonderful interplay between experimental and theoret- 
ical work will be exploited and presented in some detail. The currently emerging picture suggests that 
the properties of all protons in the nucleus are within the domain of experimental scrutiny including 
those with high momenta. The accompanying improved theoretical understanding of these properties 
of protons in the nucleus has also acted as a catalyst for the renewed study of the nuclear-matter 
saturation problem. 

The paper starts in Sec. |21 with the introduction of the relevant formalism for the determination of 
the one and two-particle propagators. By employing the equation of motion formulation it becomes 
straightforward to introduce the concept of self-consistency necessary for a proper description of the 
dynamics as it is encountered in nuclei and nuclear matter. Sec. El is devoted to a brief pedagogical 
exploration of the relation between the experimental data accessible in (e,e'p) and (e,e'2N) reactions 
and the relevant theoretical quantities introduced in Sec. |21 A detailed analysis of Green's function 
calculations for nuclear matter is presented in Sec. El We limit the discussion to zero-temperature results 
and do not cover neutron matter. This section starts with a brief introduction of the material and gives 
a short review of the first-generation results for spectral functions to establish a framework for the 
assessment of the most recent fully self-consistent results. These early calculations solved the scattering 
problem in the medium by propagating mean-field (mf) nucleons in the medium which is equivalent 
to the so-called Galitski-Feynman treatment ^7]. A similar approach is adopted to generate spectral 
functions for A hyperons which are compared to the corresponding results for nucleons. Particular 
attention is given to recent developments involving the fully self-consistent inclusion of SRC. The results 
for this second generation of spectral functions are discussed together with an analysis of the scattering 
process in the medium that underlies the properties of individual nucleons in nuclear matter. A detailed 
discussion of healing properties of nucleons in the medium is presented in order to resolve the paradox 
that arises when the original discussion of healing fSl ^| is confronted with recent (e,e'p) data. The 
consequences for the nuclear matter saturation properties including these recent developments concludes 
this section. 

Section |21 is devoted to a discussion of the results for the single-particle (sp) propagator in finite 
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nuclei. After a brief introduction we start this section with an overview of the early calculations for the 
spectral strength distribution. The influence of LRC is illustrated by a direct comparison of the spectral 
strength distribution with the data. The dual role of SRC in determining the distribution of the sp 
strength is then clarified in detail. The important influence of both low-energy particle-particle (pp), 
hole-hole (hh), and particle-hole (ph) correlations in the form of collective behavior in determining 
the low-energy sp properties is emphasized to provide the basis of a new development involving the 
Faddeev technique applied to the many-body problem. This new technique is then reviewed and shown 
to incorporate the correct way to include these correlations simultaneously. State of the art results 
for ^^O demonstrate that further improvement of the description of LRC is required for a complete 
description of the experimental data. This is further illustrated by analyzing the collective phonons 
in ^^O that form the ingredients of the Faddeev summation. A comparison with the experimental 
excitation spectrum of ^^O demonstrates that further improvement of the theoretical description of this 
spectrum is necessary. A first step in this direction is presented which takes the self-consistent coupling 
of single phonons to two phonons into account. 

A discussion of recent calculations of two-proton removal cross sections follows including a com- 
parison with experimental data. Recent experimental developments will be briefly highlighted which 
elucidate the properties of all protons in the domain of the traditional nuclear mean-field. These results 
are also discussed in Sec. Altogether with recent attempts to determine the properties of high- momentum 
protons in the nucleus. In Sec. IHl the results are summarized in as far as they clarify the properties of 
protons in the nucleus. A brief comparison with other correlated many-body systems is also presented 
there. A final summary and outlook for future work is given in Sec. 13 

2 Formal development 

2.1 Single-particle Green's function 

The nuclear systems to be discussed in this review will be described as a collection of nonrelativistic 
nucleons interacting by means of a two-body interaction V that describes nucleon-nucleon (NN) scatter- 
ing data up to pion production threshold. One can generate a zero-order approximation of the system 
under study by introducing an appropriate mean-field potential U and splitting the Hamiltonian into an 
unperturbed one-body part Hq = T + U and a residual interaction Hi = V — U. The total Hamiltonian 
can then be written as 



where we have chosen to use a set of sp states {a} that diagonalize Hq with eigenvalues cj^ (cq) 

are the creation (destruction) operators of a particle in the state a, Vap,-^& the antisymmetrized matrix 
elements of V , and Ua^p correspond to the matrix elements of tf . In general, the best choice of Hq (and 
therefore of the basis {a}) depends on the symmetry properties of the system under study. For infinite 
nuclear matter translational invariance suggests the use of momentum eigenstates without any need 
for choosing an auxiliary potential, while for finite nuclei one needs to employ a potential {e.g. of the 
harmonic oscillator or Woods-Saxon type) to localize the nucleons. Obviously, the spectrum {e^^} can 
contain both discrete and/or continuum states and we therefore use the summation symbols in Eq. 
to include both a summation over the discrete part of the spectrum as well as an integration over the 
continuum part. This convention will be used throughout this section. 

Three-body forces are known to produce important effects on details of nuclear spectra [201 1^ - 
These forces can in principle be added to Eq. ^ and included in the formalism. Since the calculations 
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to be discussed in this paper do not explicitly include their effects, we will therefore not consider them 
in the following but we will discuss some related issues in Sec. 14.71 

The experimental information concerning sp properties of a many-body nuclear system is usually 
obtained by probing the ground state \^q) of a nucleus with A particles by the addition or removal of 
a nucleon as discussed in more detail in Sec. 01 The one-body (sp) propagator (or two-point Green's 
function) associated with the state I^^q) embodies this information and is defined according to - jH] 

9a,{t-t') = -.«|r[c„(t)4(0]K), (2) 

where T[...] represents the time-ordering operation and cjj(t) and Ca{t) now correspond to operators in 
the Heisenberg picture (with h = 1). For r = t — t' > 0, Eq. Q incorporates the probability amplitude 
for adding a particle in a state (3 and removing it from a state a after a time r. Similarly, a hole can 
be generated in a and annihilated from f3 after a time t' — t > 0. 

The information contained in the sp propagator gapiT) becomes more transparent after Fourier 
transformation from the time to the energy formulation. After separating the positive and negative 
time contributions in Eq. (0) one includes the completeness relations for A ± 1 nucleons to obtain 

p+oo 



27r 



— oo 



^ a;-£+ + zr/ ^ u; - ~ ir] 

which is know as the Lehmann representation of the sp propagator [22] • Here, and in the following, the 
indices n and k enumerate the eigenstates of the system with A + 1 and A — 1 particles, respectively. 
As discussed above, the summations are intended to represent sums over the discrete spectrum and 
integrals over the continuum part. The poles = E^^^ — E^ and = Eq — E^^^ in Eq. Q 
correspond to the excitation energies of the 1-body systems with respect to the ground-state energy 
Eq. Equation Q includes information on the transition amplitudes for the addition and removal of a 
nucleon in the numerator and these excitation energies in the denominator. This feature illustrates the 
wealth of information in the sp propagator that can be compared to experimental data. Experimentally, 
this information is typically obtained in the form of the one-hole spectral function which is related to 
the imaginary part of g{uj) by 



S^{uj) = - Im g^ai.^) 



c^\^^)\U{lo - el) i^<ep, (4) 



for nucleon removal. The corresponding particle spectral function is given by 

Sai^) = -- Im gaa{^^) 
TT 

= T.\('^n-'VaK)\'S{co-et) u;>et, (5) 

n 

and in this form is mostly of theoretical relevance. S'^{uj) ( 5*^(0;) ) represent the probability for removing 
(adding) a particle from (to) an (in) orbital a while leaving the residual system in a state with energy 
—UJ [uj) relative to the ground state Eq of the system with A particles. The Fermi energies Sp and Ep 
denote the minimum excitation energy needed to remove or add a particle and are given by 

Ep = E^-Ef (6a) 
4 = E^^'-E^- (6b) 
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For infinite systems with a ground state that does not exhibit pairing these quantities are equal in the 
thermodynamic hmit and will be referred to as the Fermi energy ep = £p = £p- In finite nuclei a 
considerable difference exists between Sp and Ep for the closed-shell nuclei considered in this review. 

An important experimental result, when considering the knockout of nucleons from a mf orbit, is the 
reduction of the observed cross section with respect to the prediction of the independent-particle model 
(IPM) jTT] - [in] as discussed in Sec. El This is due to the fact that these orbits are only partially occupied 
in the correlated system. From the theoretical point of view, this effect can be quantified by looking 
at the normalization of the overlap integral between the initial and the final states, {'^f~^\ca\'^o) |IE] 
(see Sec. 01 for more details). Considering the transition to a state k of the A — 1 particle system one 
may define the theoretical spectroscopic factor as 

a 

Obviously, this quantity is also contained in the hole spectral function given in Eq. 

The information on nuclear structure that is contained in Eqs. dH) and (0) can be fruitfully compared 
to the case of an uncorrelated mf system. When the residual interaction is neglected in Eq. (^, the 
wave function of the system is an eigenstate of the unperturbed Hamiltonian Hq. The corresponding 
ground state, is a Slater determinant with all the lowest sp orbitals fully occupied up to the hole 

Fermi energy Ep, while those above are empty. The unperturbed propagator associated with \^q) is 
given by 

(0), . ^ f eja-F) ejF-a) 

[UJ - E^ + 17] UO - E%-n] 

where 6{a — F) ( 0{F — a) ) is equal to (1) if a is an occupied state and it is 1 (0) otherwise. The poles 
correspond to the unperturbed sp energies and the forward and backward going contributions 
contain either orbitals that are outside or inside the Fermi sea F, respectively. We note that in order 
to consider the whole sp spectrum, one needs to consider both the particle and hole states. The hole 
spectroscopic factors associated with Eq. (jH]) are exactly one for those states that belong to the Fermi 
sea and zero otherwise, reflecting the occupancy of the orbitals in |$^). A deviation from unity occurs 
in a finite system when the center-of-mass motion is properly treated |23] ■ 

This simple mf picture is modified when one includes many-body correlations, generated by the term 
Hi in the full Hamiltonian. Comparing Eqs. Q and (|H1), one can still find a correspondence between 
the unperturbed energies of particle orbitals and the excitation spectrum {n} of the system with A + 1 
nucleons and between the hole energies and the spectrum {k} oi A — 1 nucleons. However, a larger 
number of states can appear in the correlated systems due to the mixing of sp excitations with more 
complicated configurations. As discussed in Sec. 12.21 these more complicated states are at least of two- 
particle— one-hole (2plh) or two-hole— one-particle (2hlp) character. In the case of nuclear systems, 
the mixing of these states is generated by the presence of both collective modes of the system and 
the relevance of short-range effects associated with the nature of the underlying two-body interaction 
V. This mixing generates a large number of sp fragments (or poles in Eq.Q ). At the same time, it 
reduces the total sp strength of each peak of the IPM and redistributes it among all the fragments of 
both particle and hole type, over a wide range of energies. The reduction of the spectroscopic factors ((7|) 
from 1 can therefore be taken as a direct measure of correlations in the system. 

The occupation number of a given sp particle state in the correlated ground state can be obtained 
from the hole spectral function 

= {'^t\cic.\^o) = E K^^'l^-l^o)^ = Hdu S^^iuj) , (9) 
k 

where the integral ranges over all the fragments of the orbital a that are spread over different quasihole 
energies, as discussed above. Correlations have the effect of reducing the occupation of the orbitals in 
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the Fermi sea (which would be 1 for the unperturbed state |$^)) and to partially populate states that 
were originally empty. A measure of the emptiness of the orbital is obtained by integrating the particle 
spectral function 

P + OO 

da = {^t\cAK) = E K^^'l^"l*o)r = / du S^^iuj) . (10) 

Using the fermion anticommutation relations, one obtains the sum rule na + = 1, valid for any 
orbital a. This is an important result since it connects information on the removal and the addition 
of a particle to the system which refer to quite different experimental processes. It also supports the 
suggestion that both quasiparticle and quasihole states, together, form a natural space in which to 
study sp properties 

It is also worth noting that the knowledge of the sp propagator allows the evaluation of the expec- 
tation value of any one-body operator O according to 

«l O K) = J20^,, I ^.e-'gp^{uj) . (11) 

Moreover, it is also possible to compute the energy of the ground state by means of the Midgal- 
Gahtski-Koltun (MGK) sum rule ^ 2^ 

= (vi>^|^|^^) = ]- r dcu V ({a\f\P) + u S^^ -Im gf^^ico) , (12) 

which is an exact result when only two-body forces are considered. Equation (fT^ is noteworthy since it 
includes the expectation value of the potential energy, which is a two-body quantity, solely in terms of the 
one-body propagator. This contribution is obtained by integrating over all possible separation energies 
the product of this separation energy multiplied by the corresponding spectral strength (summed over all 
relevant sp orbits). This result is very important because it shows how the details of the hole spectral 
distribution strongly influence the binding energy of the correlated ground state. The occupation 
number Ha yields information about the overall population of the state a but the energy integral in 
Eq. loses track of the energies at which these nucleons are located inside the nucleus. In the case of 
nuclear systems a sizable amount of strength is shifted down to very negative sp energies by short-range 
and tensor correlations, with important consequences for the binding energy of the system as discussed 
in Sec. 10 



2.2 Equation of motion method and Dyson equation 

The Green's function Q can be computed by applying the methods of quantum field theory. This path 
was developed in the fifties by Migdal and Galitski [23 1^ and by Martin and Schwinger Q and makes 
use of a Feynman-diagram expansion. The method has found extensive application in several areas of 
many-body physics including atomic, condensed matter, and nuclear physics. One possible approach 
consists in choosing the propagator g^^\uj) as a starting point and to generate a perturbative expansion 
in the interaction Hi. The diagram rules for calculating the contributions to g{uj) can be found in 
several text books (see for example Refs. [3 13 E])- In general, the expansion will contain diagrams 
that describe all the possible ways in which the unperturbed sp propagator interacts with particle 
and hole excitations inside the system. Due to the strongly repulsive character of the nuclear force a 
truncation of the perturbative expansion to a given power in the interaction Hi is quite inadequate 
and it becomes necessary to include the relevant physical processes by means of all-order summation 
techniques. One therefore introduces the concept of the irreducible self-energy S*(u;) as the collection of 
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Figure 1: Graphical representation of the Dyson equation (Eq. ()13|) ). The double line repre- 
sents the exact (dressed) sp propagator g^u), while the single line stands for the unperturbed 
propagator g^^\uj). 



all connected diagrams that are one-particle irreducible, that is, they cannot be separated by cutting a 
single propagator line. The sum of these contributions represents the effective interaction that a particle 
is subjected to when interacting with the nuclear medium as will be shown explicitly in Sec. El These 
contributions can be summed to all orders by the Dyson equation shown in Fig. ^ which generates the 
whole perturbation series as shown in Fig. |21 This Dyson equation is given by 

9a,H = gTp{u^) + E<^(^)^75M^?^/^H- (13) 

7,5 

An alternatve derivation of Eq. ()13|) can be obtained by either considering the equation of motion for 
g{Lo) j3] or working directly with the effective action functional [27 . The result is an exact formulation 
relating each many-body Green's function to the exact propagators of higher order. This generates a 
hierarchy of relations between the two-, three- and A-body Green's functions, of which Eq. (fT^ below 
is the first example j3J Ej- Employing the equation of motion for a Heisenberg operator dca{t)/dt = 
—i[ca(t),H] one obtains the time derivative of Eq. Q 0, 

d 

i-g^9ap{t - t') = 6{t~t')6o,p + elga,p{t-t') - ^U^^g^p{t - t') 

7 

+ ^ E W (-^) I T[ct(t)c^(t)c,(t)4(t')] K) . (14) 
This generates a term containing the 4-point Green's function 

9t-p%{hMhM) = «l T[c;3(t2)c.(tl)ct(t3)4(t4)] 10 . (15) 

By contracting this with the matrix elements of the two-body interaction, as in the last term of Eq. ()14|1 . 
one can relate the expectation value of V to the one-body propagator whence the MGK sum rule ()12j) 
follows. 

In general, (7^"^* can describe the propagation of either two-particle (pp), two-hole (hh) or particle- 
hole (ph) excitations depending of the ordering of its time arguments. From a diagrammatical point 
of view, this propagator is the sum of two contributions. The first is a diagram in which two different 
particles propagate, fully correlated, without any interaction among them. The second group defines 
the four-point vertex function r^"''* which includes all the diagrams in which two dressed particles can 
interact with each other, thereby generalizing the T matrix for the scattering of two free particles. The 
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Figure 2: Illustration of how the Dyson equation generates the expansion in which the 
irreducible self-energy S*(a;) is summed to all orders. 



corresponding contribution to (7^ ^* is simply obtained by adding external lines to this vertex. One 
therefore obtains 



(16) 



which is depicted in Fig. |3] 

The Dyson equation (fT^ can now be obtained by Fourier transformation of Eq. ()14|) . Employing 
Eq. dH)) one then finds an exact expression for the self-energy T,*{u) in terms of the sp propagator g{uj) 




4-pn zz 



A - 




Figure 3: Diagrammatic structure of the four-point Green's function in terms of the vertex 
function T^-p* (see Eq. (ITHj) ). 
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and the vertex By using Eq. ()16|) one obtains 





{l\U\S) + s. 



(17) 



where the term —{'y\lf\6) subtracts the auxihary potential which was added in the unperturbed Hamil- 
tonian ^^o- This result is shown in part a) of Fig. H}, The second term in Eq. (|17p is the Hartree-Fock 
contribution to the self energy, 



= E / ^ ^"7,/.^ 9,siu;) . (18) 



which represents the (energy-independent) interaction of the nucleon with the quasihole excitations 
inside the system. When all the higher-order terms of the self-energy are neglected, the solution of the 
Dyson equation g{uj) takes the same simple form as Eq. (jHI), which is equivalent to a Slater determinant 
description of the ground state. In this case the quasihole wave functions contained in g{uj) refer 
to completely filled orbitals and the iterative solution of Eqs. (fT^ and (fTHj) generates the standard 
Hartree-Fock approximation. 

The many-body correlations beyond the nuclear mean field are included in the energy-dependent 
part of the self-energy. A possible expansion of this self-energy, in terms of F^"''*, is given by the last 
term of Eq. (fT7j) and shown diagrammatically in Fig. EJi). One should note that Eq. (fT7j) is formally 
exact and therefore the latter term contains implicitly the effect of both pp, hh, and ph collective 
excitations and their coupling to the sp motion. This is of course at the prize of having to keep track of 
the different time arguments in F^~^*, which involves an awkward mathematical structure. When only 
the pp (hh) effects or the ph ones need to be studied, F"^"^* can be conveniently approximated by the 
ladder or the ring series, respectively (to be described in Sec. 12. 4|) . In these cases F^~^* reduces to a 
two-time quantity, i.e. depends on only one energy variable. In the cases where this simplification is 
not possible, Eq. (fT7|) may no longer correspond to the best approach and it becomes useful to go one 
step further in the hierachy of Green's functions by taking a second derivative in Eq. (jl4j) . with respect 
to the time argument t'. This will result in introducing a 6-point Green's function Ra/3-y,^iux(t — t') 
that includes explicitly the coupling of sp motion to 2plh and 2hlp propagation. Only the two-time 
reduction of this propagator is needed to compute the self-energy. After Fourier transformation one 
obtains. 



which is depicted by part c) of Fig. HJ The propagator R{uj) contains all the diagrams that describe 
the propagation of 2plh and 2hlp (and more complicated configurations) but is one-particle irreducible 
since these terms are generated by the Dyson equation. 

It is worth mentioning that the same development of Eqs. ()14|) to (fT7j) can be carried out by 
considering the time derivative of Eq. (j21) with respect to t'. This leads to an analogous result shown 
diagammatically in part b) of Fig. 0] Baym and Kadanoff have shown that an approximation chosen 
for F''"^* should be such that parts a) and b) in Fig. HI generate the same self-energy [2H] - ED] • With 
this symmetry requirement it is assured that the solution of the Dyson equation g{uj) satifies basic 
conservation laws, such as particle number, total energy, total momentum, and total angular momentum. 
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Figure 4: Various possible expansions of the irreducible self-energy in terms of higher order 
many-body Green's functions. In all cases the first two terms represent the mean-field 
contributions —If and the Hartree-Fock self-energy S^^. The dynamic part of the self- 
energy can be expressed in terms of the four-point vertex function r^~^*. Part a) represents 
the relation given in Eq. (fT7j) . while the part b) corresponds to an alternative derivation 
involving the time derivative of g with respect to t'. If r^~^* is approximated in such a 
way that the corresponding two diagrams in a) and b) are equivalent, the Dyson equation 
will satisfy appropriate sum rules. Part c) gives the expansion in terms of the one-particle 
irreducible 2plh/2hlp propagator R{uj), Eq (fT^ . Note that the full four-time dependence of 
p4-pt jg ngg(jg(j in principle, while in the R{u!) formulation one can specialize to a two-time 
quantity. 



2.3 Self-consistent approach 

In Eqs. ()17|) and ()19p. the term — (7|f/|(5) removes the auxiliary potential U included in g^^^'^u) (or 
equivalently in Hq, as discussed in Eqs. (P^jl - fjHTj) below). This makes the solution of the Dyson equation 
formally independent of the choice of U. An implicit dependence on Hq = T + U remains for the case 
of a standard perturbative expansion, where the irreducible self-energy S*(c<j) is expressed as a series of 
Feynman diagrams, in terms of g^^\u) and the vertices of the residual interaction [21 El- However, this 
is not the case when one uses the approach of the equation of motion discussed in the previous section 
since this provides for an expansion of the self-energy in terms of the exact propagator g{u!) and fully 
dressed higher-order Green's functions. The use of Bethe-Salpeter-like equations to evaluate the 4- and 
6-point Green's functions results in an expansion of the self-energy only in terms of g{uj). While this 
g is not exact for a given approximation, it is self-consistent and for this reason a propagator obtained 
from this method can be referred to as a self-consistent Green's function (SCGF). This self-consistency 
feature has the advantage of generating a formalism in which the actual excitations that propagate 
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Figure 5: Iteration scheme used to reach self-consistency between the input propagator and 
the result of the Dyson equation. At each step in the iteration loop new approximations for the 
one-body propag ator © and for T^'-p* ^ or R{uj) are obtained until self-consistency is obtained. 



in the system are employed as the basic degrees of freedom, instead of a mf approximation of the 
nuclear orbitals. Of course, this result is obtained at the price of an increased complexity in practical 
calculations. In general, one seeks an adequate approximation of the self-consistent propagator g{uj) in 
terms of the relevant degrees of freedom that are important in the system under study. 

In the nuclear case, there exist strong couplings between the sp degrees of freedom to both low-lying 
collective excitations and high-lying states, the latter coupling being due to the strong short-range 
repulsion in the nuclear force. These mechanisms generate the sizable fragmentation of sp strengths 
which is observed experimentally. In order to obtain a detailed description of the resulting nuclear 
structure it may be necessary to include the details of the fragmented one-body propagators directly in 
the construction of the nucleon self-energy. 

The practical approach of SCGF theory consists in starting with an approximation for the input 
Green's function (usually a Hartree-Fock or other IPM propagator). This can be used to construct an 
approximation for the irreducible 2plh/2hlp propagator R^u\,ai3'y{^) or the four-point Green's function 
fj^Qj-Q "vvhich one derives the irreducible self-energy S*(c(j). Then, the solution of the Dyson Eq. ()13p 
will give an improved approximation to the self-consistent sp propagator that can, subsequently, be 
employed in a next iteration step. For infinite nuclear matter the self-energy should include both pp 
and hh propagation to properly account for short-range effects and particle number conservation. This 
can be done by summing the ladder equation for r^~^* and employing Eq. ()17|) . In the case of finite 
nuclei both pp (hh) and ph collective motion can be relevant and one is forced to seek a more complex 
expansion, based on Eq. (fTIJ|) . This may require the computation of both the pp (hh) and ph propagators 
and then to couple them through a Faddeev expansion . As shown in Fig. El the whole procedure 
should be iterated until consistency is obtained for two successive solutions {i.e. between the input 
propagator and the solution that is generated). We note that the Hartree-Fock formalism, in which 
-p4-pt R{uj) are neglected, is the simplest realization of this formulation. 

The attractive feature of the SCGF method is not only restricted to the fact that the effects of frag- 
mentation are included in the calculations. Also, solutions for the one-body, two-body and higher-order 
propagators are generated all at the same time while their mutual influence is taken into account. In 
addition, the comparison with experimental data can be performed for all ingredients of the calculation 
yielding important clues for further improvements of the chosen approximation scheme. We stress again 
that the final self-consistent solutions will not depend in any way on the choice of Hq. Nevertheless, a 
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good choice for the unperturbed propagator g^^^ can provide a useful starting point for the first iteration. 

The self-consistent propagator g{uj) enters the construction of higher-order Green's functions through 
the use of Bethe-Salpeter equations (BSE). We give a brief description of this approach by considering 
the three-times polarization propagator, 

n„/3,75(tl - t',t2 - t') = -I (^^1 T[4(t2)c„(tl)c;(t')Q(t')]K) 

+^ (^^1 T[4(t2)c„(ti)]|vI/^)K|ct(Oc,(t')|^^) , (20) 

which governs the response of the system to an external perturbation. Working with the equation of 
motion method, one finds that Eq. (j^UI) is the solution of the following Bethe-Salpeter equation 

na/3,75(i^i — t' , h — t') = —igajih — t')gsp{t' — ^2) 



+ 



j dts j dti j dU j dte ^ 5-0^(^1 - h)guf3{U - ^2) 



X ^S,A,(^3, U; h) ^xnAh - 1\ h - 1') , (21) 

where the ph irreducible kernel K^^, represents all the basic interactions that a quasiparticle and a 
quasihole can undergo. The latter can be obtained from the exact irreducible self-energy in terms of a 
functional derivative, 

K,,Jt,,tM) = ^-^—^ . (22) 

In the language of Feynman diagrams this corresponds to taking the diagrammatic expansion of ^*{uj) 
in terms of the propagator g{uj) and then generating the contributions to K^'^ by removing such a sp 
line in every possible way. Baym and Kadanoff have shown that a conserving approximation of the 
polarization propagator 11 is guaranteed if both K^^ and g{uj) are derived from the same approximation 
of the self-energy (by means of Eqs. (j22I) and (fT5|) . respectively) [2H1 EOl- From a practical point of 
view, this approach is not always possible to implement since relatively simple approximations to S*(co') 
can generate a very large number of contributions to K^^. Fortunately, actual nuclear systems do not 
necessarily require the inclusion of all of these terms to obtain reasonable results. Nevertheless, the 
above construction can be used as a guideline to help in choosing the approximation that is suitable 
for the problem under study. The resulting K will involve an expansion in terms g{uj) and r^~^* 
or R{(jj). The fact that this perturbative expansion is based on the fragmented sp propagator — 
and not on its unperturbed approximation g^'^\ijj) — , acts to renormalize the interaction kernel K and 
drastically improves the convergence properties, while still allowing to keep track of the relevant physical 
ingredients in the calculation. The simplest realization of this scheme consists in considering only the 
bare interaction potential V in either the pp or ph channels. The result is the usual random phase 
approximation (RPA). This approximation is described in more detail in the next section since it is 
relevant to most of the calculations discussed in this review. 



2.4 Ring and ladder approximations 

By taking the limit t2 = t m Eq. ()20j). one is left with the usual two-time polarization propagator 

'^ap,'ys{^ ~ t')- This can be Fourier transformed to its Lehmann representation, 

where the poles correspond to the excitation energies = E:^ — Eq of the states of the system 
whereas the residues contain the corresponding amplitudes to determine transition probabilities from 
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Figure 6: Feynman diagram representation of the DRPA equation for the polarization propagator 



the ground state to these excited states. It is important to note that this information is included for all 
those states that can be excited by means of a one-body operator, therefore also those with important 
two-particle— two- hole (2p2h) (or even more complex) admixtures. In Eq. (j23p . the forward- and the 
backward-going contributions are identical apart from a time reversal transformation. 

Analogously, one can define the free ph propagator n-^(a;) by coupling two excitations of particle 
and hole type, not interacting with each other. Graphically, this corresponds to two dressed sp lines 
propagating with opposite direction in time. By applying Feynman diagram rules [Oj, one obtains the 
following Lehmann representation for n-^(ci;). 



- E 



H\cl\^t-'){H\Ca\^^~'') {^^^'\c\\H){^t-'\cs\H) 



n,k 



n,k 



UJ 



(24) 



where the forward- and the backward-going contributions describe the propagation of a particle-hole 
and hole-particle type excitation, respectively. 

In order to account for the collective excitations present in the A-body system, one needs to take 
into account the interactions between the two lines propagating in 11^ (a;). The simplest approximation 
consist in considering only the contribution of the bare two-body potential V. This corresponds to 
choosing K'^^^^ = Vas,i3-y in the Bethe-Salpeter Eq. (j^D) and results in the standard RPA approximation, 

n,^,^,(w) = Ui^^^.iu) + ^U,ui^) ^Me,.p ^pe,,s{uj) , (25) 

which is depicted in Fig. [Ul in terms of Feynman diagrams. Equation (j25p implicitly generates a series 
of diagrams in which the particle and the hole interact any number of times. The resulting expansion 
is depicted in Figs. |7|and|Hl in terms of g^^\uj) propagators. 

A simpler approximation consists in considering only the forward-going propagation in Eq. ()25|) . 
i. e. to include only the first term of Eq. This corresponds to the so-called dressed Tamm- 

Dancoff approximation (DTDA), depicted in Fig. [7| (where an explicit time ordering is assumed and 
the dressing of sp Green's functions is implied). When the full Ilf{uj) is employed in Eq. (j25j) . one 
obtains the dressed random phase approximation (DRPA). The relevant diagrammatic expansion is 
the one given in Fig. |H| Clearly, the whole TDA series is contained in the RPA one. Besides this 
contribution, the backward-going component of Il-^{uj) generates diagrams containing more that a single 
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Figure 7: Diagrammatic expansion of the ph propagator 11(0;) in TDA. An implicit time ordering 
is intended for these diagrams. Dressing ah the sp lines would result in the DTDA. 



ph excitation at each time. These refer to contributions in which additional intermediate ph states are 
present and eventually annihilated by two-body interactions. Consistent with the aim of performing 
a self-consistent calculation, Eqs. and (^3)) have been formulated in terms of the full, dressed, sp 
propagator g{u) given in Eq. (j2)). When an unperturbed, IPM, Green's function g^{uj) (jSI) is used as 
input, the approximations described here reduce to the standard TDA and RPA results [32j. 

In general, the interparticle distance between nucleons inside the nucleus is of the order of the nucleon 
size itself so that the nuclear medium can be considered a rather dense system. Besides, the nuclear 
effective interaction remains substantial. A natural consequence is that the "vacuum fluctuation" 
contributions that are included in the RPA expansion play an important role in the description of 
nuclear systems. It must be realized, however, that the intermediate n-particle-n-hole excitations 
depicted in Fig. |H1 do not represent all the possible diagrams related by Pauli exchange of two sp lines. 
The philosophy underlying the use of RPA relies on the assumption that the corrections coming from 
Pauli correlations sum up in a random way and tend to cancel each other (whence the name "random 
phase") jnn] - IHEj- As a result, the Pauli principle may only be violated slightly and the improvement 
obtained by including RPA correlations is generally more relevant. In such cases, the RPA is a useful 
approximation. It must be kept in mind, however, that in systems containing strong collective and/or 
pairing correlations the effects of Pauli violation may sum up in a coherent way, invalidating the standard 
RPA approach. Typically this situation is signaled by the appearence of a pair of complex eigenvalues 
of Eq. (j25|l together with diverging solutions for the ph spectroscopic amplitudes. 

We also note that using dressed propagators g{uj) in Eq. results in including implicitly diagrams 
that go beyond the bare Iplh configurations already at the level of n-^(a;). This is due to the self-energy 
insertions in each particle and hole line. The fragmentation of ^'(ci;) acts once again to renormalize the 
free propagator Ilf{uj) thereby adding stability to the DRPA approach, as compared the the bare RPA 
one. 

The saturation energy of nuclear matter, as well as the spectral distribution at high missing energies, 
is governed by short-range effects and two-body correlations. Moreover, some particular nuclei present 
strong pairing effects. In these cases, one is faced with two-body correlations in the system and a proper 
description of the physics calls for a direct evaluation of the interaction between two particles in the 
medium. From the experimental point of view, the ideal approach to study such effects is represented 
by addition or removal of the correlated pair (see Sec. EJ. The information required to study such 
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Figure 8: Diagrammatic expansion of the ph propagator n(w) in RPA. This contains the whole 
phTDA series plus diagrams involving the inversion of the propagation time. Dressing all the sp 
lines results in the DRPA. 



processes are also contained in the four-point Green's function, p*. One can define the two-particle 
propagator g^^{t — t') as a two-time reduction of Eq. (fT3j) 

9ii,s{t-t') = (26) 

Its Lehmann representation, 

contains the spectroscopic amplitudes for a transition from the ground state |\E'o ) to the eigenstates of 
the systems with A ±2 particles and the respective excitation energies with respect to Eq. 

As for the polarization propagator, the zero-order contribution consists of two noninteracting lines, 
propagating in the same time direction. The Lehmann representation for the dressed case is given by 



ni,n2 



J, u ^ 



where the forward and backward parts refer to the (independent) propagation of two particle and two 
hole hues, respectively. The two-particle dressed RPA (ppDRPA) equation is given by 



M= {9'JU^) - 9'p';U^)) + Wi',.H, (29) 

^ivpe 



9al3,-y5{ 



where the symmetry factor ^ is required by the Feynman rules jB]. 

Equation ()29|) is depicted in terms of Feynman diagrams in Fig. El and generates a series of ladder 
diagrams in a similar fashion as higher order terms are generated in Fig. |H| We note that this approach 
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Figure 9: DRPA (ladder) equation for the two-particle propagator g {oj). 



is equivalent to computing the two-body wave function of a pair propagating inside the medium, while 
the Pauli blocking effects are being taken into account in a similar fashion to the Bethe-Goldstone 
equation EZj . The ppDRPA has two advantages over that approximation. First the details of the 
sp particle fragmentation, included in Eq. ()28|). allows for a partial propagation below the Fermi level. 
This is due to the fact that in the correlated system the Fermi sea is partially depleted. The dressed 
result therefore corresponds to propagation with respect to the correlated ground state and not the IPM 
as in the Bethe-Goldstone equation. Second, the advantage of an RPA approach over the TDA one is 
that it accounts for the propagation of hh excitations in the ground state of the system. Without this 
latter feature, the resulting self-energy would not yield the fragmentation of the sp strength below the 
Fermi energy. As in the ph case, the corresponding ppDTDA approximation is obtained by considering 
propagation only in one time direction. This corresponds to including only the first term from the left 
of Eq. ()28|l into Eq. ()29|1 when solving for the pp part of g^^ , therefore neglecting the hh contributions. 

In general, both the two-particle and the polarization propagator represent different time orderings 
of the same 4-points Green's function, see Eqs. and (j^ . By making use of Eq. (fT^ and (fT7|) one 
can include the effects of ph or pp(hh) motion directly in the nuclear self-energy. 



3 Relation to experimental data 

In this section we will explore the connection between the information contained in various propagators 
and experimental data. Detailed work on this subject has been presented in Ref. jSH] and we refer to 
that paper for a discussion of the relation between elastic nucleon scattering and the particle part of 
the sp propagator. In this paper the focus is on the experimental properties that are probed by the 
removal of nucleons. Nevertheless, it will be necessary within the theoretical treatment to simultaneously 
consider the addition and removal aspects of the propagators as discussed in Sec. El Before making 
this connection it is useful to demonstrate that the Dyson equation (Eq. (fT!?|) ) yields a Schrodinger-like 
equation with corresponding interpretation of the removal or addition amplitudes appearing in Eq. (jS)) • 
We will discuss here the case when the spectrum for the A ± 1-particle system near the Fermi energy 
involves discrete bound states which applies to a finite system like a nucleus. An appropriate form of 
the Lehmann representation in this case is given by 

^<*j'|c„|*r'>(*^+i|ct|4'i?> /--^ ^(<I'„''|c„|*^i>(*^i|4|<I'i?> 

f-^(-) = E E^TtTn, + L E^wn^, 



+ 2^ = ■■ + / de^ ^— : , (30) 



k 



where the continuum energy spectrum for the Ail systems has been included and the corresponding 
energy thresholds are denoted by A change of integration variable was also used to obtain this form 
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of the Lehmann representation introducing the integration variables for the continuum energies in the 
form = E^~^^ — Eq and e~ = Eq — E^~^ , respectively. For the unperturbed propagator one will 
encounter sp energies associated with Hq that are different from the poles of g{uj). By exploring the 
equations of motion of the two-body propagator, it is possible to show that a Lehmann representation 
also exists for the exact self-energy and has different poles from the ones for g{oj) 0. These features 
can be used to identify the residues of discrete solutions in the Lehmann representation from the Dyson 
equation. In practice, one may proceed with taking the following limit of the Dyson equation for the case 
of the hole part of the propagator (without generating contributions from the poles in the self-energy 
or the noninteracting propagator) 

- el) [g^p{u) = g^^^iij) + <)(^) S;,(.;) gs^iu) }. (31) 
This limit process generates an eigenvalue equation of the following kind 



lim 



7,(5 



where 



4- = {^k\caM)- (33) 



Since the Dyson equation can be written in any sp basis, one may choose the coordinate representation 
with sp quantum numbers r,m for the position and spin projection, respectively, to solve Eq. (jH^ . In 
this basis one obtains 



d^l^l d^1^2 glm nmii^k) ^nmi r2m2 (^fc ) ^r2m2- (34) 

mi,m2 

Equation ()34|) can be rearranged by inverting the unperturbed propagator according to 

d^r (rW| - Hq |rm) c/^^ r^mA^k) = Sm',m^S{r' - ri). (35) 

The corresponding operation on yields 

J2 fd'r (rW| e, - Ho |rm) z',- = ^ + ^ - f/(r')| 4w> (36) 



where U is assumed to be local and spin-independent for simplicity. Combining these results yields the 
explicit cancellation of the auxiliary potential U and the following result 



^ k^ 
' 2m 

mi 



^ ^ jd 1^1 ^ rm rimi (^fc)^rim,i ^k^rmi (37) 



where the notation S'* has been used to signify that the U contribution has been removed. This 
equation has the form of a Schrodinger equation with a nonlocal potential which is represented by the 
self-energy S'*. Note that an eigenvalue can only be obtained when it coincides with the energy 
argument of the self-energy. An important difference with the ordinary Schrodinger equation is related 
to the normalization of the quasihole "eigenfunctions" The appropriate normalization condition 

is obtained by considering terms beyond the pole contribution in the Dyson equation. This result is 
most conveniently expressed in terms of the sp state (itself normalized to 1) which corresponds to the 
quasihole wave function z^^. In other words, one can use the eigenstate which diagonalizes Eq. (j2Zj) 
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with eigenvalue , to express the normahzation condition. Assigning the notation aqh to this sp state, 
one obtains with z^" = (^^"^j c„^, |^^) 

) (38) 

where the subscript qh refers to the quasihole nature of this state and the fact that for states very 
near to the Fermi energy with quantum numbers corresponding to fully occupied mean-field states the 
normalization yields a number of order 1. This result is equivalent to the definition given in Eq. (jT)). 



dr.' 



dE 



3.1 Spectroscopic strength from the (e, e'jy) reaction 

In order to make the connection with experimental data obtained from knockout reactions, it is useful 
to consider the response of a system to a weak probe. The hole spectral function introduced in Sec. |21 
can be experimentally "observed" in so-called knockout reactions. The general idea is to transfer a 
large amount of momentum and energy to a proton of a bound nucleus in the ground state. The proton 
is then ejected from the system, and one ends up with a fast-moving particle and a bound {A — 1)- 
particle system. By observing the momentum of the ejected particle it is then possible the reconstruct 
the spectral function of the system, provided that the interaction between the ejected particle and the 
remainder is sufficiently weak or treated in a controlled fashion, e. g. by constraining this treatment 
with information from other experimental data. 

We assume that the A-particle system is initially in its ground state, 

l*.) = W>, (39) 
and makes a transition to a final A-particle eigenstate 

l*/)=4l*n-'>' (40) 

composed of a bound {A — l)-particle eigenstate, |\E'^~^), and a particle with momentum p. 
For simplicity we consider the transition matrix elements for a scalar external probe 

A 

Pil) = ^exp{iq-rj), (41) 

i=i 

which transfers momentum qr to a particle. Suppressing other possible sp quantum numbers, like e.g. 
spin, the second- quantized form of this operator is given by 

p(g) = Yl ■ ^) 4v = Yl (42) 
p,p' p 

The transition matrix element now becomes 

(^/i kq) i^i) = Yl <^""'l «p4'V-<j \^t) = Y <^""'l v,p«p'-<j + «j,'«p'-<jap |^^> 
p' p' 

^ (^n">p-,W>- (43) 

The last line is obtained in the so-called Impulse Approximation, where it is assumed that the ejected 
particle is the one that has absorbed the momentum from the external field. This is a very good 
approximation whenever the momentum p of the ejectile is much larger than typical momenta for the 
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particles in the bound states; the neglected term in Eq. PHj) is then very small, as it involves the removal 
of a particle with momentum p from I'^o)- 

There is one other assumption in the derivation: the fact that the final eigenstate of the A-particle 
system was written in the form of Eq. ()4U|) . i.e. a plane- wave state for the ejectile on top of an {A — 1)- 
particle eigenstate. This is again a good approximation if the ejectile momentum is large enough, as 
can be understood by rewriting the Hamiltonian in the A-particle system as 

= E ^ + E ^(^' = + ^ + E (44) 

1=1 i<j=l i=l 

The last term in Eq. (j44p represents the Final State Interaction, or the interaction between the ejected 
particle A and the other particles 1 . . . A — 1. If the relative momentum between particle A and the 
others is large enough their mutual interaction can be neglected, and Ha ~ Ha~i + p\/2m. The result 
given by Eq. (pSj) is called the Plane Wave Impulse Approximation or PWIA knock-out amplitude, for 
obvious reasons, and is precisely a removal amplitude (in the momentum representation) appearing in 
the Lehmann representation of the sp propagator (see Eq. ^ ). 

The cross section of the knock-out reaction, where the momentum and energy of the ejected particle 
and the probe are either measured or known, is according to Fermi's golden rule proportional to 

da ^Y,^{u: + E,- Ef)\ f\ p{q) |^,) (45) 

n 

where the energy-conserving 5-function contains the energy transfer uj of the probe, and the initial and 
final energies of the system are Ei = Eq and Ef = E^~^ +p^/2m, respectively. Note that the internal 
state of the residual A — 1 system is not measured, hence the summation over n in Eq. (j45j) . Defining 
the missing momentum Pmiss and missing energy E^iss of the knock-out reaction as^ 

Pmiss =P-q (46) 

and 

Emiss =p^/2m-uj = E^ - Et\ (47) 
respectively, the PWIA knock-out cross section can be rewritten as 

da ~ 5^5(E_,-i?o^ + E^i)|(vI/^i|ap_|OP 

n 

ShijPmiss^ E^niss) • (48) 

The PWIA cross section is therefore exactly proportional to the hole spectral function defined in Eq. 
This is of course only true in the PWIA, but when the deviations of the impulse approximation and 
the effects of the final state interaction are under control, it is possible to obtain precise experimental 
information on the hole spectral function of the system under study. Although the actual (e,e'p) 
experiments involve more complicated one-body excitation operators than the one considered here in 
this simple example, the basic conclusions are not altered |lUj . 

An even more realistic description of the proton knockout reaction again identifies the external 
electron probe and the corresponding virtual photon as represented by a one-body excitation operator. 

d = Y,{oi\0\l3)alap. (49) 

a/3 

^We will neglect here the recoil of the residual A — 1 system, i.e. we assume the mass of the A and A — 1 system to 
be much heavier than the mass m of the ejected particle. 
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The response of the system to such an external excitation operator is described by the polarization 
propagator which has a Lehmann representation given by Eq. (j23|) . The transition probability induced 
by O from the ground state to an excited state is given by 

(^n I O K> ' = 5^ 5^ (7l O \b) a\a, K> («| O W {<\ |^^>* • (50) 

a/3 75 



This result demonstrates that the numerator of the first term in Eq. (jzHl) contains the relevant transition 
amplitudes for a given state n to evaluate this transition probability. In general, there are important 
correlations between the ph states, in particular at low energy where collective surface vibrations and 
giant resonances occur. At higher excitation energy and momentum transfer, these collective coherence 
effects tend to disappear. In this domain, one may therefore write the polarization propagator as the 
product of two dressed sp propagators as in Eq. ()24|) . Replacing 11 by this dressed II-^ one neglects 
contributions where the particle and hole interact but includes all other correlations associated with 
the full dressing of the removed particle in terms of the corresponding removal amplitude (or spectral 
function) and the corresponding dressing of the particle that will ultimately be detected in the (e,e'p) 
experiment. The wave function of this dressed particle corresponds to the addition amplitude in co- 
ordinate space and can be associated with an optical model wave function [SH]. From this analysis it 
becomes therefore clear that one may use empirical information associated with the elastic scattering 
of protons in terms of optical potentials to describe this wave function of the outgoing proton. Clearly, 
it is this interpretation that clarifies the assumptions that underly the standard analyis of the (e,e'p) 
reaction ^T] - JHl- In the practical analysis of an (e,e'p) experiment it is conventional to find a local 
potential well (mostly of Woods-Saxon type) which will generate a sp state at the removal energy for 
the transition that is studied. This state is further required to provide the best possible fit to the 
experimental momentum dependence of the cross section (with proper inclusion of complications due to 
electron and proton distortion) ^T] - ^Hl- The overall factor necessary to bring the resulting calculated 
cross section into agreement with the experimental data, can then be interpreted as the spectroscopic 
factor corresponding to the "experimental" quasihole wave function according to Eq.()38p. 

The resulting cross sections obtained at the NIKHEF facility are shown for four different nuclei 
in Fig. ^1 It is important to realize that the shapes of the wave functions in momentum space 
correspond closely to the ones expected on the basis of a standard Woods-Saxon potential well (or more 
involved mf wave functions). This is itself an important observation since the (e,e'p) reaction probes 
the interior of the nucleus, a feat not available with hadronically induced reactions. 

While the shapes of the valence nucleon wave functions correspond to the basic ingredients expected 
on the basis of years of nuclear structure physics experience, there is a significant departure with regard 
to the integral of the square of these wave functions. This quantity is of course the spectroscopic factor 
and is shown in Fig. ^2 for the data obtained at NIKHEF ^1] . The results shown in Fig. ^2 indicate 
that there is an essentially global reduction of the sp strength of about 35 % which needs to be explained 
by the theoretical calculations. This depletion is somewhat less for the strength associated with slightly 
more bound levels. An additional feature obtained in the (e,e'p) reaction is the fragmentation pattern of 
these more deeply bound orbitals in nuclei. This pattern is such that single isolated peaks are obtained 
only in the immediate vicinity of the Fermi energy whereas for more deeply bound states a stronger 
fragmentation of the strength is obtained with larger distance from sp. This is beautifully illustrated 
by the (e,e'p) data from Quint (see Sec. EJ. Whereas the 3si/2 orbit exhibits a single peak, there 
is a substantial fragmentation of the 1/ strength as indicated in this figure. Additional information 
about the occupation number of the former orbit is also available and can be obtained by analyzing 
elastic electron scattering cross sections of neighboring nuclei jlU]. The actual occupation number for 
the 3si/2 proton orbit obtained from this analysis is about 10% larger than the quasihole spectroscopic 
factor 11111121 and therefore corresponds to 0.75. All these features of the strength need to be explained 
theoretically. This will be attempted in the material covered in later sections. 
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Figure 11: Spectroscopic factors from the (e,e'p) reaction as a function of target mass. Data 
have been obtained at the NIKHEF facihty ;14i| . 
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3.2 Information from the (e,e'2N) reaction 

Suggestions to explore SRC in two-nucleon emission reactions go back to the work of Gottfried |43j . 
More recently, theoretical work has focused on the possibility to utilize the (e,e'2N) reaction to probe 
nucleon-nucleon correlations ^3] - [IH]- Practical descriptions of this reaction have been developed by 
the Pavia group ^Tj - [SUI- Proceeding in a similar vein as in the analysis of the (e,e'p) reaction, which 
yields information about the one-nucleon (removal) spectral function, one may hope to learn about the 
two-nucleon (removal) spectral function in two-nucleon emission processes. The emission of two protons 
is particularly promising for studies of SRC since the effect of meson-exchange currents and isobars is 
not expected to dominate the cross section under suitable kinematic conditions |1H]- 

Experiments have been carried out for ^^C jHI] and ^^O (discussed in more detail in Sec. 15. 7|) to 
explore the feasibility of gaining insight into nucleon-nucleon correlations in finite nuclei using the 
(e,e'pp) reaction. Triple coincidence measurements involving protons with large initial momenta seem 
particularly suitable to provide information on SRC. The scattered electron is then expected to transfer 
a virtual photon to one of these two protons which have large and opposite momenta and therefore a 
relatively small center-of-mass momentum. This strong correlation results from hard collisions due to 
the strong repulsive core of the NN interaction. When one of the protons is removed by the absorption 
of the virtual photon, its partner will also leave the nucleus under the assumption that the energy 
transfer is mainly to the hit pair (the residual nucleus stays at a low excitation energy) jJTj. It is 
therefore hoped that, if the coupling of the virtual photon to one nucleon is the dominant mechanism, 
the (e,e'pp) process may be exploited as a useful tool to investigate these short-range correlation effects 
(see also Ref. [S3)- 

One of the goals of triple coincidence measurements is to illuminate the features of the interaction 
between two nucleons before their knockout of the nucleus and their subsequent detection. Early ex- 
periments on ^^C employing the (e,e'pp) reaction already assumed in the analysis of the data obtained 
at NIKHEF jH^l EZj that the virtual photon is coupled to one of the detected protons. This approxi- 
mation can be understood by considering the transition matrix element of the nuclear charge operator 
in momentum space 

P(g) = e J] aj,,^qap/ (51) 

(with spin implicit in the summation) between the initial state \'^q} and an approximate final state of 
the form 

|vI/^w) = <,at,Jvl>;^-2). (52) 

This final state contains two plane wave protons and an exact state |\E'^~^) for the system with two 
protons removed. In calculating the matrix elements of the transition charge operator one obtains 

(^^na/|p(g)|^o) = (^n'>p?,Vi-<jl^o) + (^n I «P^-<7 Vj ^0 ) , (53) 

assuming that one of the detected protons absorbs the momentum q. To obtain the contribution to the 
cross section one requires the square of this matrix element. This yields four terms, each representing 
a particular term of Eq. (jMj) for the two-nucleon spectral function. An important ingredient in the 
description of the two-nucleon knock-out reaction is the two-hole spectral function defined by 

n 

where may denote the (0^) ground state of the target system (for example ^^O) and |\E'J^~^) the 
n-th excited state of the residual nucleus (^^C). In Eq. fl3^ . aj, (op) represent the addition (removal) 
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operators of a nucleon with momentum p (spin and isospin are implicit). It is therefore clear that the 
appropriate expression for the transition probability is to consider the quantity 



where the vector q is the momentum of the virtual photon. 

More realistic implementations of this reaction model naturally require the consideration of the 
distortion effects of the outgoing particles. Nevertheless, this simple representation of the reaction 
clarifies the intrinsic importance of two-nucleon spectral functions in describing two-nucleon knockout 
reactions. More details related to actual comparison of calculations with experimental data will be 
presented in Sec. 15.71 

4 Theoretical calculations for nuclear matter 

The study of nuclear matter has remained a prominent field of study in recent years. We will focus in 
this section on recent developments related to self-consistent Green's functions and put these results 
in perspective with regard to other many-body approaches to the relevant quantities of interest where 
appropriate. All many-body techniques that depart from the free nucleon-nucleon interaction are re- 
quired to treat the short-range repulsion of this interaction in an appropriate manner. In perturbative 
methods one must sum the infinite set of ladder diagrams that fulfills this requirement. This was re- 
alized by Brueckner long time ago jHS]- Since then many different implementations have been studied 
with perturbative methods to treat short-range correlations. Some of this work has been reviewed in 
Refs. |54[ P . The main emphasis in the present review is on the progress that has been made recently 
in implementing the self-consistency condition on the single-particle (sp) propagators in consort with 
the summation of ladder diagrams for the effective interaction. According to the Lehmann represen- 
tation of the sp propagator (see Eq. © ) knowledge of the propagator is equivalent to knowledge of 
the particle and hole spectral functions given in Eqs. (0)) and (0). The first nuclear-matter spectral 
functions in SCGF theory were obtained for a semirealistic interaction by employing mean-field (mf) 
propagators in the two-body scattering equation [SHI EE]- Self-consistency was limited to the sp spec- 
trum obtained from the real part of the on-shell self-energy. The corresponding ladder equation includes 
both particle-particle (pp) and hole-hole (hh) propagation and is sometimes called the Galitski-Feynman 
equation J7j. Hole spectral functions using correlated basis function theory (CBF) were first obtained 
in Ref. jHZ! while particle spectral functions were reported in Ref. [HB]- Spectral functions for the full 
Reid potential jHS] were obtained in SCGF by employing a self-consistent gap in the sp spectrum [^01 - 
to avoid pairing instabilities in the ^5*0 and ^Si-^Di channels The first solution of the ladder 
equation using fully dressed sp propagators jMI was obtained by employing a parametrization of the 
spectral functions |65j . Important consequences of treating this dressing of the sp propagators include 
a strong reduction of the in-medium cross section and the disappearance of the pairing instabilities at 
normal nuclear- matter density jHE]- The latter conclusion was recently confirmed in detail in Ref. |(i7j . 

Several approaches towards the goal of self-consistent calculations of nuclear-matter Green's func- 
tions have recently been reported in Refs. [HH] - [ZH]- In Ref- [H] the average effect of the spectral 
strength distribution was included in the determination of the sp spectrum thereby avoiding pairing 
instabilities. In addition several realistic interactions were employed to study the dependence of spectral 
functions on the choice of the nucleon-nucleon (NN) interaction. The Ghent group has concentrated on 
a discrete representation of the the sp propagator which yields a scattering equation containing several 
discrete poles but otherwise similar to the usual one based on mean- field propagators [7T1 [731 US] • This 
type of approach will likely generate the same results for the energy per particle as the continuous 
version implemented in Refs. [72] and [HZl [77j for an appropriate choice of the discretization scheme. 



+ S'\p[,p',-q,p[ 



q,p'„E) + S^'^p', - q,p'2,p[,P2 - q,E) 
q,p'„E) + S^\p[,p', - q,p[,p', -q,E), 



(55) 
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Initial results from the various groups indicate an important change of the saturation properties with 
respect to Brueckner Hartree Fock (BHF) calculations with a continuous choice for the sp spectrum. 
These results have been obtained for an updated version of the Reid potential in Ref. [73]. In 
Refs. [THl EZ] a separable representation of the Paris interaction (HOI |HI] has been used to generate 
self-consistent propagators including short-range correlations by employing a direct discretization of 
the spectral functions and self-energy as a function of energy. Also in this work ^| a repulsive effect 
is observed for the energy per particle leading to less binding as compared to the conventional BHF 
approach with a continuous sp spectrum. These results therefore also generate a corresponding sub- 
stantial reduction of the saturation density. Several realistic interactions were studied in Ref. [H2]- In 
this recent work a new argument about the importance of separating long and short-range contributions 
to binding in nuclear matter was introduced. The hope is that such new developments will lead to new 
insights into the long-standing problem of nuclear-matter saturation. 

A discussion of the formalism necessary to accomodate a self-consistent treatment of short-range 
correlations in nuclear matter is presented in Sec 14.11 We begin the discussion of results by reviewing 
some early calculations of first-generation spectral functions in Sec. 14.21 A useful digression is made in 
Sec. 14.31 where spectral functions for the A hyperon in nuclear matter |j83| are discussed. These spectral 
functions are also based on the solution of the corresponding ladder equation in the medium employing 
mf intermediate propagators. A useful comparison between nucleon and A spectral functions is then 
possible, illuminating the simililarities and differences between these nuclear constituents. The difficulty 
of propagating fully dressed particles in the ladder equation is associated with the necessity to treat all 
off-shell aspects of the scattering process in the medium. In Sec. 14. 41 a discussion of the scattering process 
of dressed particles [HI] is presented that has wider implications for the conceptual understanding of 
the properties of nuclei in terms of reconciling their violently interacting constituents with aspects of 
the simple shell-model picture. Results that compare typical quantities that characterize scattering 
for dressed, mf, and free particles are presented in Sec. 14.51 Results that illustrate the properties of 
fully self-consistent Green's functions in compariosn with those of the first generation are presented in 
Sec. 14.61 The last subsection 14.71 concerned with nuclear matter properties is devoted to the discussion 
of the consequences of self-consistent treatments of short-range correlations for the understanding of 
nuclear saturation properties. 



4.1 Formalism for self-consistent Green's functions in nuclear matter 

The formalism for the improved determination of the sp propagator in nuclear matter using a self- 
consistent scheme which includes the full treatment of short-range correlations will be outlined below. 
In symmetric nuclear matter one may take advantage of the invariance properties of the system to write 
the sp propagator in the following way 

r,.'^%^,r,.'^#^. (56) 

where Sp{k,uj) {Sh{k,u)) represents the probability distribution to add (remove) a nucleon with mo- 
mentum k to (from) the system leaving it at an energy u. The latter functions are also referred to 
as the particle and hole spectral functions, respectively, and were introduced in Sec. 12.11 by consid- 
ering the Lehmann representation of the sp propagator in Eq. Q. The propagator can be obtained 
diagrammatically by considering its perturbation expansion [2 IHl or, alternatively, by considering the 
equation of motion which relates it to the two-body and, subsequently, higher-body propagators [Ij as 
discussed in Sec. 12.21 In the latter approach a natural introduction of self-consistent propagators is 
obtained. The link between one-body and two-body propagators is made through the self-energy which 
can be studied at various levels of approximation. The one including SRC is obtained by describing the 
two-body propagator in the same way as one would proceed in free space, i. e. sum all ladder diagrams. 



25 



Figure 12: Diagrammatic representation of the ladder equation [Eq. ()59|) ] in the medium. 
Note that fully dressed propagators are iterated between successive interactions. Only in 
this diagram the subscript "L" is used to identify the ladder summation. 



This procedure leads to the in-medium equivalent of the T matrix in free space and therefore properly 
accounts for short-range correlations. Using a notation with q, q' for relative and K for the conserved 
total momentum, the resulting two-body propagator can be written as 

g''{q,q';K,n) = g'/ {q, q'; K ,n) + g'/ {q; K ,n) J d'q" {q\V \q") g'' {q" , q'; K ,n) (57a) 

= gj\q, q'- K, Q) + g'/ [q- K, Q) {q\ T{K, Q) \q') gf{q'- K, Q), (57b) 

where 

g'/{q, q'; K, = 6{q - q') g^iq; K, (58) 

is the noninteracting but dressed two-particle propagator which conserves the relative momentum as 
expressed by the 6 function in Eq. (j^Hj) . The presence of exchange terms in Eqs. ()57p and ()58|) and 
possible summations over spins and isospins is hereby acknowledged but suppressed in the presentation. 
Eq. ()57b|) links the two-particle propagator with the four-point vertex function F shown explicitly in 
Fig. El 

The four-point vertex function or effective interaction F then contains the summation of all ladder 
diagrams in the following way 

{q\ r{K, Q) \q') = {q\ V \q') + J d'q" {q\ V \q") g'/iq"; K, Q) {q"\ F(K, Q) \q') (59a) 

= {q\V\q') + {q\AT{K,n)W). (59b) 

The ladder equation is shown diagrammatically in Fig. ^1 It is convenient to decompose the effective 
interaction in two parts as in done in Eq. ()59b|) which separates the energy-dependent part from the 
energy-independent NN interaction V. Employing this separation, one obtains the nucleon self-energy 
in the form shown diagrammatically in Fig. |^) with U discarded and the full vertex function replaced 
by the corresponding ladder approximation. The diagram with the closed loop corresponds to the 
energy-independent contribution involving the NN interaction 

Sy(fc) = J ^ {kk'\ V \kk') n{k'), (60) 

with 

rep 

n{k) = I duj Sh{k,uj) (61) 



oo 



representing the correlated momentum distribution. Note that explicit summation over spin and isospin 
projections have been suppressed and a notation with individual particle momenta k and fc' is used for 
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the two-body interaction in Eq. (jUUI) . Both the effective interaction and the self-energy fulfill important 
dispersion relations which are helpful in devising strategies to achieve a numerically tractable self- 
consistency procedure. For the energy-dependent part of the effective interaction one has 

= {q\AT^{K,n)\q') + {q\AT^iK,n)\q'). (62b) 

The decomposition of the effective interaction into forward and backward-going contributions as ex- 
hibited in Eq. (jHH) is essential to obtain a proper construction of the self-energy |S3]. Indeed, using 
Eqs. fl56|) and ()62a|) one obtains the remaining contributions to the self-energy in the following way 



/A3u pep 
J duo' {kk'\ AT^iuj + u') \kk') Sh{k', uj') 

r f]3u poo 

- / / duo' {kk'\AT^{u + uo')\kk')Sp{k\uo') (63a) 

J (27r)^ Jep 

= AS|(A;,c<;) + AS|(A;,cj). (63b) 

As for Eq. (|6(J|) the notation for the two-body matrix elements of AF involve individual momenta. The 
result of Eq. also implies the following dispersion relation for the self-energy 

^ ^ n^ ^ T , M ^{k,uo') 1 ^ ,ImS(fc,u;) 
E(k,uo) = Tjv(k) / duo' — - + - / duo' ^-J—L (64a) 

= Sy(A;) + ASj,(A;,cj) + AS|(A;,cj). (64b) 

Using this self-energy in the Dyson equation one obtains 

gik,uo) = g^'^\k,uo) + g^^\k,uo)J:{k,uo)g{k,uo) (65a) 
uo — e{k) — Re S(fc, a;) + zim S(A;, uo) (ccu\ 
[uo - e{k) - Re S(A;, uo)f + (Im S(A;, uo)f ' 

where the explicit form of the noninteracting propagator 

^'*:*;>. + (66) 

(jj — e[k)+irj u) — e[k) — IT) 

with 

e{k) = f- (67) 
2m 



is used to obtain the result of Eq. (j65bj) . The diagrammatic version of the Dyson equation is shown in 
Fig. n Although it is now possible to solve the ladder equation with mean-field propagators without 
an angle-averaging procedure 86| , it is necessary for a practical implementation using fully dressed 
propagators to make this approximation. To obtain the ladder equation in a partial wave basis one 
therefore proceeds by an angle-averaging procedure for the noninteracting but dressed two-particle 
propagator 

g'/{q; K,n) = ^ j duo g{K/2 + q, + uj)g{K/2 - q, n/2 - u) 

Sj,{K/2 + q,uj,)Sp{K/2-q,uo,) 



-oo 

duOl / duO' 



2 n 

\l — UOi ~ UO2 + IT] 



-00 J —00 



Q — uoi — UO2 — ii] 

27 



Eq. (j68|) shows that the angle-averaging is confined to the numerators. It is therefore possible to consider 
the imaginary part of for this purpose. For Q > 2eF one has 

1 f°° 

Im gy{q; K,n) = / du Sp{K/2 + q, n/2 + uo) Sp{K/2 - q, - u) (69) 

and for Q < 2eF the corresponding result is given by 

Im gYiq; K,Q) = - [ ^ du Sh{K/2 + q, Vl/2 + uj) Sh{K/2 - q, Q/2 - u) (70) 

The convolution integrals in Eqs. and ()7()|1 also provide a practical route to obtain the real part of 
g^ by means of the following dispersion relation 



The cycle of steps necessary to perform a self-consistent calculation of the sp propagator is now complete. 
Starting from a reasonable set of spectral functions, the construction of the noninteracting but dressed 
two-particle propagator first requires the convolution integrals given in Eqs. (jU^ and (f7n|) . In a next 
step the dispersion integrals in Eq. (f7T|l can be used to obtain the real part of the this propagator. 
After an angle- averaging procedure the resulting propagator only depends on the absolute values of the 
relative and total momenta. It is therefore permissible to solve for the ladder equation in a partial wave 
basis 

{q£\T^^^{K,n)\q'e') = {qe\V^^^ \q'i') (72) 

POO 

+ E / ^^PP' lO a'/ip-^ n) {pi"\ T^'^iK, n) \q'i') . 

According to Eq. ()63b|) one only needs the diagonal elements of the imaginary part of F (properly 
adding the different partial wave contributions) and the original input spectral functions to obtain the 
imaginary part of the self-energy. The dispersion integrals in Eq. ()64j) together with Eq. ()60|) then allow 
the construction of the complete real part of the self-energy. The solution of the Dyson equation is then 
straightforward. One finally obtains the new spectral functions from 

bJk,uj) = 5 ^ (73) 

for energies above Sp and from 

bh[k,uj) = ^ ^ (74) 

for energies below Sp. At this point the cycle can be repeated until by some numerical criteria the input 
sp propagator is sufficiently close to the output sp propagator. One of the most important quantities 
to emerge from the spectral functions is the energy per particle [221 

E 2 f d^k f e \ ^ n ^ /-..N 
duJ h UJ bh[k,uj), (75) 



A pj (27r)3 7_^ \2m 
where p is the density both related to the Fermi momentum by 
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and the self-consistent momentum distribution 



2 

~2 



dkk^ n{k) (77) 







since self-consistency guarantees particle number conservation j2Hll2n]- One may also separately consider 
the kinetic energy per particle 

rf n poo 7,2 

dkk^ —n{k). (78) 



A 7r2 7g 2m 

The potential energy can then be obtained by considering the difference between Eqs. fl75|) and ()78|). 
This general outline requires further details related to a suitable numerical implementation. Some of 
these aspects will be considered in Sec. 14.61 where results for fully self-consistent spectral functions are 
discussed. 



4.2 Spectral functions obtained from mean-field input 

The first generation spectral functions were based on the solution of the scattering equation (Eq. (|72|l ) 
employing mf propagators in the construction of with corresponding 5 function spectral functions. 
The corresponding result for the self-energy given in Eq. (jfiHj) also simplifies and is given by 



SAr(A;,a;) = J -0^ {kk'\ ATf\Lu + sik')) \kk') e{kp - k') 

{kk'\Arf{uj + e{k'))\kk') eik'-kp), (79) 



{2ny 

where a superscript has been attached to F to identify that it was obtained from a scattering equation 
with mf propagators. A serious difficulty arises when propagating mf propagators corresponding to 
a continuous sp spectrum in the ladder equation when a realistic NN interaction is employed 
Indeed, so-called pairing instabilities arise in the ^Si-^Di and ^5*0 partial wave channels in certain 
density regimes jHS]- Self-consistency with unperturbed propagators in the ladder equation can be 
achieved while solving the above difficulties and still maintaining a connection between the sp potential 
and the full S by making the following choices Below kp, let e{k) be the first moment of Sh{k,uj) 
with respect to energy (the contribution to the total energy at that momentum) divided by the number 
of particles in the system with that momentum. One cannot extend this intuitive prescription to k > kp 
by replacing Sh with 5*^ since Sp contains a significant fraction of strength up to GeV energies and would 
result in unrealistically large values for e{k) when k > kp. Therefore, above kp, Sp is replaced with 
the quasiparticle (qp) contribution to Sp. This is equivalent to setting e{k) equal to the qp energy 
(discussed in more detail below). One then obtains 

£{k) = \ ^ an\ k<kp (80) 



/ du Shik, 



and 



The above prescription for the sp spectrum requires calculation of the full energy dependence of S in 
each iteration step and is therefore more computationally intensive. However, it naturally contains a 
gap large enough to prevent the pairing instabilities. More details related to this prescription of self- 
consistency can be found in Ref. The qp energy alone has no appreciable gap, but near the Fermi 
momentum there is considerable spreading of hole strength to energies below the qp energy. 
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Figure 13: Imaginary part of the self-energy as a function of energy below the gap for three 
momenta, k = 0.01 (solid), 0.51 (dotted), and 2.1 fm~^ (dashed). 



Some results from this scheme will be reviewed now to establish a benchmark for the second- 
generation spectral functions discussed in Sec. 14.61 One characteristic feature of the resulting imaginary 
part of the self-energy is shown in Fig. ^| Comparing different values of k for these energies below the 
Fermi energy indicates that ImE^ becomes weak while its energy range enlarges in a smooth manner 
as k increases. The energy range for this ImE^ is determined by using the analysis of Ref. [22] • This 
analysis shows that for each k there is a minimum energy above which 2hlp states (of mf type) can 
mix in the self-energy. So momentum conservation and the constraint on the location of 2hlp energies 
due to their mf character are responsible for the energy range observed in Fig. ^1 This wide range of 
energies will be available to the hole spectral functions discussed below. 

For momenta near kp one finds that the energy dependence of Sh and Sp is dominated by the qp 
peak characterizing the extent to which noninteracting features are maintained. Each momentum has 
an associated qp energy which is the solution of 

egp{k) = ^ + ReS(fc, Sgpik)). (82) 

The spectral function displays a peak at Sqp because of the vanishing term in the denominator of Eq. (|7S|) 
or (fTIj) . The qp-peak itself is represented by 

Sqp{k,uj) = -— tryn.\,x,(u\\'? (83) 



TT 



~ eqp{k)f + {Z{k)W{k)y 



(9ReS(fc,w)' ■ ^ 



where W{k) = Im T,{k,eqp{k)) and 

Z{k) = {l-{ -—-^^-^ \ } (84) 

is the strength contained in the peak. These first-generation SCGF calculations show that near kp 
the qp peak contains around 70% of the total sp strength, while an extra 13% is contained in the 
background composing the remainder of hole strength. The background is uniformly distributed across 
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Figure 14: Illustration of the decreasing width of the quasiparticle peak in the spectral 
function for three momenta below kp given by k = 0.48, 0.79, and 0.93 fm'^ The vertical 
lines indicate the position of the gap in the sp spectrum that results from employing Eqs. (|80|) 
and (jHD). 



several hundred MeV below ep corresponding to the range of Im S but depends significantly on the value 
of k as shown in Fig. El The final 17% of the strength has moved to energies greater than ep including 
a significant high-energy tail in Sp discussed momentarily. Farther below kp this picture breaks down 
as the qp peak melts into the background resulting in hole strength which is spread over a much wider 
range of energies. To the extent that spectral functions are described by the qp approximation, the 
excitations in nuclear matter are like those of a Landau Fermi-liquid as illustrated in Fig. El This figure 
contains several hole spectral functions for momenta below the Fermi momentum at kp = 1.36 fm~^. 
Notice that as k ^ kp this peak becomes extremely sharp due to the vanishing ImS in Eq. (ff^ . The 
infinite-lifetime character of such excitations is made possible by the loss of phase-space available to 
the states in S near the Fermi energy. This is essentially the same argument used by Landau in more 
general terms to develop the microscopic foundations of Fermi-liquid theory jHHj • 

In Fig. El the particle spectral function is plotted for three different momenta, k = 0.79, 1.74, and 
5.04 fm~^ as a function of energy. All momenta below kp have the same high-energy tail as the dotted 
curve for k = 0.79 fm~^ in Fig. El For momenta larger than kp a quasiparticle peak, which broadens 
with increasing momentum, can be observed on top of the same high-energy tail. The results therefore 
display a common, essentially momentum independent, high-energy tail. The location of sp strength 
at high energy simply means that the interaction has sufficiently large matrix elements to compensate 
energy denominators encountered in the ladder equation. For this particular interaction a significant 
amount of strength is found at high energy. This result was of course already anticipated a long time 
ago [Hi. 

A quantitative characterization of the missing sp strength for = 0.79 fm~^ shows that the integrated 
strength accounts for 17% of the sp strength. This is in agreement with the sum rule since the integrated 
hole strength provides 83% of sp strength. The strength in the interval from 100 MeV above the Fermi 
energy to infinity amounts to 13% with 7% residing above 500 MeV. To understand the influence of 
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Figure 15: Particle spectral functions at = 0.79 (dotted), 1.74 (solid), 5.04 fm~^ (dashed). 
All three spectral functions converge to the same tail at high energy. 



the tensor force on this distribution, a calculation of the ladder equation was performed in which the 
tensor coupling in the ^Si-^Di coupled channel was switched off. In this case, the integrated sp strength 
amounts to 10.5% and should be regarded as resulting from pure short-range correlations. In Ref. [nOl 
it is shown that the tensor force moves the additional 6.5% of strength to the first 1000 MeV above 
the Fermi energy. This is consistent with CBF calculations of the momentum distribution which show 
depletions of a similar size due to tensor correlations |90|j. 

Figure ITHl exhibits the graph of the occupation probability, n(A;), or the number of particles in the 
ground state of the system with sp quantum number k. Near A; = 0, n{k) becomes fairly constant with 
a value of 0.83. From the discussion of Sp, roughly 1/3 of the 17% depletion is due to the effect of 
tensor correlations in the ladder equation. Another 1/3 is due the to high-energy tail in Sp at energies 
above 500 MeV. Results from other many-body methods such as Brueckner theory [211 1221 ^^"^ CBF 
theory jHTj using other realistic interactions give very similar occupation near A; = 0. In Refs. [21] and 
|92j . 0.82 is reported for the Paris potential. Older CBF calculations for the Urbana vu interaction give 
0.87 whereas more recent CBF results [201 §i^^ 0.83. All these calculations for different interactions 
using different methods give a strikingly similar result for n(0). This was encouraging since it implies 
that non-relativistic many-body calculations are under control in the region where one would like to 
compare to finite nuclei. 

In contrast to n(0), the occupation at kp varies significantly between methods. This variability is 
most easily expressed in terms of the discontinuity in n{k) at kp, Z{kF) = n{kp) — n{kp). For the 
result shown in Fig. Zlkp) = 0.72 and for the CBF calculation Zlkp) = 0.70. However, for the 
Paris interaction, Z{kp) = 0.35 and Z{kp) = 0.47 has been obtained in Refs. [21] and respectively. 
The extra depletion in n{k) as k —>■ kp arises from the enhanced ability of the sp state to couple to low 
lying 2plh excitations as its energy approaches the 2plh states. Also, the discontinuity depends on the 
level at which pairing correlations are included in the calculation. In the case of a paired system Z{kp) 
would be zero. 
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Figure 16: Occupation probability for nuclear matter at equilibrium density obtained by 
integrating hole spectral functions obtained with mean- field propagators as input. 



4.3 Spectral functions for A-hyperons obtained from mean-field input 

Hypernuclei, especially those with one A hyperon, have been studied for a long time [Hniini]- When a 
A hyperon is placed in a nucleus or nuclear matter it will interact with the nucleons in its environment. 
As a result of these strong interactions, the A becomes correlated with nucleons in the medium. The 
study of the properties of the A hyperon in an environment of nucleons aims to answer a number of 
fundamental questions related to the properties of strange particles in the nuclear medium. Considerable 
attention has been given to the potential energy the A experiences in the nucleus and the corresponding 
single-particle (sp) energies. From this experimental work it becomes clear that the A hyperon is less 
strongly bound to nucleons than either a proton or a neutron. Such sp properties have been studied 
theoretically for finite nuclei by several groups [HH] - [^Hl- The general conclusion from the experimental 
work is that the A hyperon experiences a potential well in the nucleus that has the familiar Woods- 
Saxon shape with a depth of about 30 MeV for a wide range of heavier nuclei. Global sp properties of 
the A hyperon can be studied in nuclear matter. Results of such calculations have also been reported 
by several groups jSH] - |lU5j . 

The correlations of the A hyperon in nuclear matter have typically been studied at the level of its 
average binding or sp energy. Full propagation of the A including the determination of its complex 
self-energy has not been reported sofar. In view of the relevance of the properties of a strange particle 
in a nuclear system, it seems timely to study the properties of a A hyperon when it is embedded in such 
a nuclear system. As in the case of NN interactions, typical hyperon-nucleon (YN) interactions |lU(jj 
- incorporate substantial repulsion at short distance. As discussed above, the consequences of 

this strong interaction can be accounted for in the framework of the Green's function formalism by 
including the proper treatment of these short-range correlations (SRC) in the form of ladder-diagram 
summation for the hyperon-nucleon interaction in the medium [G matrix). The effects on the dynamical 
single-particle properties of the A can then be explored by evaluating the complex self-energy of the 
A in nuclear matter. The solution of the Dyson equation for the A then yields again information on 
the net binding of the A in nuclear matter, but also determines the distribution of spectral strength 
for its addition to the nuclear-matter ground state as a result of SRC. Such a spectral function for a 
A hyperon is shown in Fig. [T7| The addition of a strange test particle opens the door for quantitative 
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Figure 17: Spectral function for a A with k = 100 MeV/c |HHllll2j . The vertical dashed line 
indicates the position of a delta-function spectral distribution for the limiting case of a free 
particle. Because of the 30 MeV binding for a A in nuclear matter, it is convenient to shift 
the horizontal axis by 40 MeV for plotting on a log scale. 



comparisons with spectral functions obtained for nucleons discussed in the previous section. The weaker 
YN interaction is expected to result in similar but less extreme modifications to the spectral distribution. 
However, the presence of the A hyperon also requires consideration of its heavier sibling, the isospin 
one S hyperon. The two hyperons have a small enough mass difference that a coupled-channel problem 
must be solved. In the case of the A propagator some qualitative changes occur in comparison with the 
case of the nucleon propagator. The result corresponding to Eq. (|56|) now becomes 



Since no Fermi surface for A hyperons is considered here, it is not possible to remove a A from the NM 
ground state. As a result, the propagator only contains the probability amplitude for adding a A particle 
with momentum k at an energy u'. The energy threshold at which it becomes possible to add a A to the 
NM ground state is denoted by e^. One expects this energy to be accessible only for a A at rest (/c = 0). 
If the A were added to a free Fermi gas NM ground state and would not be correlated otherwise, this 
threshold energy would simple be the kinetic energy of a A with zero momentum and would therefore 
correspond to zero energy. Based on previous work by other groups jlOU I103j one expects the actual 
value of the threshold energy at normal NM density to be around —30 MeV indicating the substantial 
attraction a A experiences in NM. Since the A can only propagate as a particle, the spectral strength 
above the threshold energy must integrate to 1. This condition on the A spectral function is given by 



This result is quite different from the corresponding nuclear case where one has to consider the strength 
in the particle domain to supplement the contribution given by Eq. ()6H) . As in the case of nucleons one 
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Figure 18: Imaginary part of the A self-energy for k = 100 MeV/c. The broken curves 
represent contributions to the overall self-energy from the ^Si (dash) and ^Sq (dot) partial 
wave channels. 



can obtain the spectral function from the imaginary part of the propagator 

SA{k;u) = Im g\{k;u!) u! > e^. (87) 

TT 

The sp spectral function as shown for example in Fig. El involves the overlap between the simple 
physical state given by 

aLl^o) (88) 

and the complicated eigenstate 

l^^(E)). (89) 

which includes all interactions between the A and the nuclear medium. The extent to which there is 
overlap illustrates how well the A sp state survives intact in the medium. For the case of no interactions 
between the A and the nucleons, the overlap is perfect, since the state of Eq. (jHHj) is an eigenstate in 
this situation. This is evidenced by the 6 function spectral distribution appropriate for a free particle, 
as indicated by a dashed line at the kinetic energy in Fig. El Interactions between the A and nucleons 
are responsible for the transition from the simple 6 function structure to the more complex distribution 
of sp strength realized in NM. The mechanism behind the spreading of sp strength can be understood 
as the mixing of a sp state at a given energy with 2plh states which span a continuum of energies. 

Although the sp state is no longer an eigenstate of the many-body Hamiltonian, its quantum numbers 
are still conserved by the interaction. The total strength associated with the original sp state, though 
fragmented, is fixed. This is reflected in the sum rule of Eq. (jHBjl . Details of the strength distribution 
are determined by the density of 2plh states which increases with energy and the strength with which 
the interaction couples these states to the unperturbed sp state. This information is summarized in 
the imaginary part of the self-energy, shown for example in Fig. El Note that the decomposition in 
partial wave contributions emphasizes the dominance of the ^5*1 channel. A nuclear matter calculation 
for nucleons similar to this one jHOj yields a particle spectral function shown in Fig. El for a momentum 
just above kp- The Z-factor obtained from this calculation is Z^ikF) = 0.72, which is substantially 
reduced compared to .^a(O) = 0.87 for a similar A qp state. These two momentum values are compared 
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Figure 19: Nucleoli particle spectral function (solid) for k = 316 MeV/c with lambda spectral 
function (dashed) at A; = 60 MeV/c for comparison. 



because each qp sits at the lowest possible excitation energy for a cjp in the respective systems. In 
Ref. [nn]; this depletion of the qp strength is explained in terms of couplings to 2hlp states, which 
moves approximately 10% of the sp strength to energies below ep, and coupling to 2plh states, which 
distributes another 18% to higher energies in the particle domain. The corresponding fraction of sp 
strength in the particle domain is 13% for the A, compared to 18% for nucleons. Turning off the AN-T,N 
coupling in the Nijmegen soft core (NSC) potential reveals |109j that tensor effects are responsible for 
almost half of the reduction in the A qp strength consistent with the nucleon case discussed in Sec. 14.21 
A value of Z^ikp) = 0.94 is obtained when coupling to T,N states is cut off. 

Away from the qp peak, at high-energy, the size and structure of the spectral function is primarily 
determined by two factors. The density of 2plh states increases like cu^/^ at high energy. This growth 
in spectral strength with energy is moderated by the strength of the coupling to these high-energy 
states. A A with a reasonably low momentum couples to a nucleon hole state only with a low relative 
momentum. The high-energy AN two-particle states couple most strongly to high relative momentum 
and the strength of the potential matrix elements between these two states depends on the short-range 
characteristics of the two-body interaction. A harder core allows a stronger coupling between states and 
correspondingly more spectral strength at high energy. The fact that structure in the high-energy region 
of the spectral function is primarily determined by the short-range behavior of the two-body interaction 
should be tempered by the knowledge that the short-range part of baryon-baryon interactions are poorly 
known. Typical potentials are designed, within whatever model, to fit only low-energy experimental data 
which does little to constrain the details of the repulsive core. This situation can be taken in two ways. 
On the one hand, the high-energy tail of the spectral function is just as uncertain in detail as the core of 
the interaction from which it is derived. On the other hand, it is also just as experimentally inaccessible 
and any observable which can be related to the detail of the tail in the spectral strength distribution 
could be used to gain insight into the behavior of the bare two-body interaction at short-range. In Fig. QUI 
the similarities of the tail of the spectral strength for different momenta is illustrated. We note that 
this common feature of the strength is identical to the case of nucleons illustrated in Fig. The only 
difference is in the distribution of the strength which reflects the hardness of the underlying interaction. 
In conclusion, one may state that the properties of a A hyperon in nuclear matter are qualitatively 
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Figure 20: Spectral function for three values of /ca show the fc- independence of the high 
energy tail for k = 10 MeV/c (sohd curve), 110 MeV/c (dashed), and 210 MeV/c (dash- 
dot). 



quite similar to those of a nucleon. The physical consequences of the underlying interactions with their 
strong repulsive cores are also quite similar. Naturally, there are quantitative differences which indicate 
that a nucleon is somewhat stronger correlated with the nuclear medium than a A. Nevertheless, a 
common high-energy tail for the sp strength and a similar role for the tensor force in reducing the qp 
strength also emphasize the common features of these strongly interacting particles in the medium. 

4.4 Healing properties of nuclear wave functions 

The conceptual understanding of strongly interacting nucleons yielding a mf shell-model (Fermi-gas) 
picture has relied heavily on the concept of the healing of the relative wave function to the nonin- 
teracting one as discussed in Refs. ^HIE]- Experimental evidence based on the (e,e'p) reaction J3] 
has demonstrated that nucleon sp motion must be understood in terms of Landau's quasiparticle de- 
scription ISH] . In turn, this requires a substantial modification of the simple shell-model or Fermi-gas 
picture. The conventional Bethe-Goldstone propagator used to determine the effective interaction in the 
medium leads to the healing property of the relative two-nucleon wave function and takes proper care 
of short-range correlations fHl ^] • However, the corresponding effective interaction is not sufficient to 
generate a nucleon self-energy which realistically describes the sp strength distribution below the Fermi 
energy both in nuclear matter and finite nuclei. Inclusion of additional terms involving hole-hole propa- 
gation as in a Galitski-Feynman propagator is essential to achieve a realistic spectral function |Mllll3j . 
The main problem with the Bethe-Goldstone propagator is that it generates no self-energy diagrams 
which couple to 2hlp states. As a result, the experimentally observed fragmentation of the sp strength 
can not be described in this approach. Inclusion of hole-hole propagation in the ladder equation solves 
this problem of describing the fragmentation of the sp strength below the Fermi energy. However, this 
inclusion destroys the healing property of the relative function since it produces a nonvanishing phase 
shift in the relative wave function for energies below twice the Fermi energy |114j when mf propagators 
are employed. The resolution of this puzzle requires the consideration of the consequences of the dress- 
ing of the nucleons for the description of the scattering process in the medium jH3j. In order to clarify 
this issue it is therefore necessary to delve in to the details of this scattering process and elucidate the 
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essential differences with the corresponding process involving free or mf particles. 

For the purpose of the present discussion it is sufficient to consider the two-time two-particle prop- 
agator as given by Eq. (jSTj). The main objective of this development is to study the consequences of 
propagating dressed particles for the description of the scattering process. All subsequent discussion 
will be based on the solution of a Lippmann-Schwinger-type scattering equation for the effective in- 
teraction which is equivalent to summing the ladder diagrams for a particular choice of noninteracting 
two-body propagator and two-body interaction. In order to clarify the difference between the conven- 
tional discussion of scattering in free space and the one necessary for the medium, it will be useful to 
cast the description completely in the language of the two-body propagator. While this is certainly not 
necessary for the scattering of free particles, it will provide a simple way to clarify the changes that are 
required to extend the description to dressed particles in the medium. 

It is important to realize that the usual results from scattering theory are obtained in the coordinate 
representation. The relevant double Fourier transform of the two-particle propagator in Eq. ()57|) is given 

by 

/^(r,r';fi) = ^ Jd'q Jd'q' e'^-^g''{q,q';n) e~'^'-^\ (90) 

where without loss of generality we consider the case that the total momentum K is zero. The transform 
of the noninteracting propagator only involves one integration due to the presence of the 6 function in 
Eq. dSHl) 

gyir,r';n) = ^ jd'q jd\' e'^^-g^ (q, q';n) e'^^'^' = ^ jd'q e^^^^^^^'^ gf {q-^). (91) 

The result for Eq. (jS7j) can then be transformed to yield 

g'\ry-n) = gy{r,r'-n) + jdhr Jd'r2 g'/ (r^r.-.n) {r^\V\r2) g'\r2,r';n) (92a) 
= g'/ir,r';n) + Jd\, ^rfV^ (r, n; fi) (ri|r(r])|r2) (r^, r'; ^]). (92b) 

Employing these equations one can arrive at an asymptotic analysis and resulting definition of the cross 
section which is equivalent to a standard analysis involving the equation for the wave function in the 
case of scattering in free space or mf particles in the medium. Before this result is developed, it is useful 
to summarize the propagator equations in the partial-wave representation. These allow the introduction 
of the phase shift which contains the relevant information to obtain the asymptotic propagator or wave 
function in this representation. 

A partial wave decomposition of the two-body propagator in Eq. (j37j) yields the corresponding 
integral equation and the relation between the propagator and the vertex function 

~ '^^h,,9)\q: n) + g^iq; fi)(g£|r^^^(fi)|gr)^f (A;'; , fi), (93b) 



-.2 



where the vertex function or effective interaction F can be obtained from the numerical solution of 
the ladder equation in a partial- wave momentum representation as given in Eq. ()72|1 . The energy Q is 
conserved and must be viewed as a variable upon which the propagator depends (it also depends on 
the total momentum in the case of the medium). The noninteracting propagator is again denoted by 
g^ and may include the dressing of the individual particles when the scattering takes place in matter. 
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The coordinate space version of Eq. (j^^ is obtained by a double Fourier-Bessel transform 

r) fOO POO 



gYsAri, r'i'; n) = - dq dq' 4' Mqrhe'iq'r') gVsT^q^. gT; fi). (94) 

The corresponding result for the noninteracting part of the propagator, represented by the first term in 
Eq. (|93p. reduces to one integral on account of the 5 function which conserves relative momentum 

2 /""^ 

gy^,{r, r'- ^) = ' dq q' Mqr)Mqr')gy (q; fi). (95) 
The Fourier-Bessel transform of Eq. ([HH]) then has the following form 

noo 
dnrl dr^rlgy^r, n; n){ni\V''^\r2n9jsAr2r , r'i'; n) 



= 6e,egy/r,r';n)+ / rfnrfrfrsr^^ (7;;,(r, n; fi)(ri^|r-^^^ (fi)|r2£>;;,,(r2, r'; fi). (96) 
Jo Jo 

When the bare two-body interaction V is local in the relative coordinate, only one integral in the first 
equality remains. The second equality can be used to study the asymptotic behavior of the propagator 
outside the range of the interaction. 

The subsequent discussion for the scattering of dressed particles in the medium requires the con- 
sideration of the two-body propagator in the medium. For this reason it is advantageous to collect the 
conventional results for scattering in free space using this language. In the case of free particles the 
noninteracting propagator in momentum space is given by 

gyiq;n) = —-^^-—-. (97) 

\l — h q'^/m + trj 

Defining the on-shell momentum by 

n = ^, (98) 

m 

one can perform the relevant Fourier-Bessel transform of the noninteracting propagator in Eq. (j95|) 
analytically (see e.g. |115p with the well-known result 

TTt 

9f/i'^, r'\ qo) = -iqo-^jeiqor<)he{qory). (99) 

The coordinate argument in the spherical Hankel function must be the larger of r and r' and is denoted 
by r> while the argument of the spherical Bessel function is the smaller and denoted by r<. For the 
current analysis it will be assumed that the interaction has a finite range, (r|V|r') = for r, r' larger 
than some tq. Substituting Eq. ()99p in the second part of Eq. (jUUj) in the case of an uncoupled channel 
for r' > r and r' > Tq yields 

9usTir,r';qo) = -iqo^je{qor)hi{qor') 

h 

+ j^j^ drirl dr2rl gfj,{r, ri; qo){ri\T^^'^ {qo)\r2) (~^^o^) jeiqor2)hi{qor') 

= -iqo^'4'f^{r\qo)hMor'), (100) 

where 

foa poo 

\JST(^. ^ \ _ ■ , / / Ji Ji M ^ . ^ \rrJST 



i^i {r;qo) = jeiqor) + / drir^ dr2r2 gf/r,ri;qo){ri\T^ {qo)\r2)j£{qor2), (101) 

Jo Jo 
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and the conventional notation T instead of T has been introduced together with the replacement of 
by go- This result demonstrates that under the given conditions the propagator separates as a product 
of a function of r and a function of r'. This result can be substituted into the first part of Eq. (j96|) to 
obtain the relevant integral equation for the wave function ip (under the condition that r' > ro) 



ipf^ir; go) = iMo^) +11 ^^i^i ^'^'i'^l gj/ir, ri, go)(ri|l^/*"^'|r2)V'/'^"^'(r2; go) 




rJST\ 



•JST( 



(102) 



^0 



which can be found in standard textbooks (see e.g. |115j for the case of a local potential). It is derived 
here to demonstrate the relation between the propagator and the wave function since for the case of 
dressed particles one has to start with the formulation in terms of propagators. 

The asymptotic analysis of the propagator can be performed by using Eq. (j^^ in Eq. (j^Bj) under 
the assumption that the propagator will be considered for r < r' while both these coordinates are 
larger than ro, the range of the interaction. Values of ri and r2 in Eq. (I96|) larger than tq yield no 
contributions to the integral. As a result, the effective interaction, T, has a range similar to the one of 
the bare interaction V. Using the relation between spherical Bessel and Hankel functions 



(103) 



one obtains the asymptotic behavior of the propagator for the case of an uncoupled partial wave channel 
from the second part of Eq. (|96|) in the following form 



9VjsTir,r';qo) 



qohe{qor')< K^q^r) + /i^(gor) x 



m 

^ /.00/.00 

1 - 2i-^qQ I / dri dr2 {ri\T/''^'^{qo)\r2)jeiqori)jeiqor2) 
n Jo Jo 



m 



-i^%hi{qor') { h}{qor) + hf{qor) 



i-2vr. (go|r/^^(go)|go) 



(104) 



In the last step of Eq. ()104|1 one can return to the on-shell matrix element of the T-matrix in momentum 
space which completely determines the outcome of the scattering process. The term in square brackets 
corresponds to the S'-matrix element in terms of which one can define the phase shift 



{qo\Sf^iqo)\qo) 
This result can be represented by 



mgo\ 



l-2vr^( ^ 1 (go|T/^^(go)|go) 



(105) 



tan 5 



jsT ^ Im (go|T/^^(go)|go) 
Re (go|T/^ngo)|go) 



(106) 



which explicitly shows that a nonzero imaginary part of the effective interaction is required to obtain a 
nonvanishing phase shift. In turn, this imaginary part of the interaction only appears for energies where 
the noninteracting propagator has a nonvanishing imaginary part. For the scattering of free particles 
this corresponds to all positive energies. By substituting the explicit form of the spherical Hankel 
functions for £ = in Eq. ()104j) one can construct the asymptotic propagator for the ^5*0 channel 
explicitly 

9'e^Soiry;io) - -^^^e^^^^^'+^^^o) sin(gor + (107) 
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The standard result for the asymptotic wave function is contained in this equation and the imaginary 
part of Eq. (jl(J7|) is simply the product of these wave functions as a function of r and r', respectively. 

To obtain the relation between the cross section and the propagator it is necessary to return to 
Eq. ()92|) and perform the Fourier transform of the noninteracting propagator (Eq. (jlTj)) in Eq. ()91|). 
This Fourier transform is given by the well-known result (again replacing the energy Q by the on-shell 
momentum go) 

gy(r,r';qo) = --^'- -. (108) 

Airh \r — r'\ 

A similar procedure as used for the asymptotic analysis in the partial wave basis can be employed to 
obtain the corresponding result for Eq. (j^^ . Whereas in the former analysis the separable form of the 
noninteracting propagator in Eq. (j^^ is valid without constraint on r and r', this is not the case here 
since Eq. fll()8j) only becomes separable in the case r' ^ r or vice versa. In the former case one can 
write Eq. (fmH|) as 



,f(r,r';go)--^^e--'-. (109) 

Substituting this result in the second part of Eq. for both r' ^ r and r' ^ r2 demonstrates that 
g^^ is separable and can be written as 

fry) zD^^O^ 

/^(r, r'; q,) = - l^^(r; go) (110) 

in the asymptotic domain. By substituting this result in turn in Eq. (j92j) one obtains the standard 
integral equation for the wave function and the appropriate formulation for the asymptotic wave function 
to obtain the scattering amplitude 

V;(r-;go) = e'^'''''-'-^ + jd'n jd^r^ g'/{r,r,-qo){n\V\r2Mr2-,qo) 

-mr'-r^ (r,n;go)(ri|r(go)|r2)e-^''«^'-^^ (111) 



e 



One may identify the origin of the motion in the direction of the negative z-axis, meaning that r' points 
in that direction, so that q = —qov' points into the positive ^-direction. If one assumes that r is also 
much larger than the range of the potential, and, therefore much larger than any contributing value of 
ri, one can use Eq. ()109p again in the second part of Eq. piip to identify the coefficient multiplying 
the outgoing spherical wave e*^'""/r as the scattering amplitude (while double Fourier transforming the 
T-matrix element back to momentum space) 



2m7r^ 



UOA) = -^:^{q'\nqo)\q), (112) 



where 6, denote the angles associated with the direction of r and q' = q^r corresponds to the momen- 
tum of the detected motion in the direction f with the same absolute value go as the initial state. The 
differential cross section for the direction (^, 0) is then simply the square of the scattering amplitude as 
given by Eq. dTT^ 

^ = l/.o(^, </')!'• (113) 

The present formulation is closely tailored to the conventional experimental situation where a collimated 
beam propagates along the ^-axis characterized by a given energy or momentum toward a target situated 
at the origin. Detection then takes place in a particular direction away from the origin characterized by 
angles 9 and 0. The only difference is that the present formulation is appropriate for the corresponding 
center-of-mass system {K = 0). 
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To obtain the phase shifts and cross sections for particles propagating in the medium with mf sp 
energies one can proceed in similar fashion. A useful reference is the work of Bishop et al. [114j where the 
introduction of the phase shift for the case of hole-hole propagation is discussed. While some details are 
different in this case, the main points of the previous analysis remain valid when both particle-particle 
and hole- hole propagation is taken into account. The corresponding mf propagator in the medium, also 
known as the Galitski-Feynman propagator, is given in momentum space by 



where, again, the case of zero center-of-mass momentum is considered. The sp energy e{q) can deviate 
from the simple kinetic energy spectrum and therefore yield a different relation between the energy Vl 
and the on-shell momentum go 



nevertheless the uniqueness of go for ^ given energy is still preserved. Although one can no longer 
evaluate the noninteracting propagator in coordinate space completely analytically from Eq. (jHSj), the 
separability of the propagator is maintained for the contribution of the pole term as in Eq. ()99|) (with a 
different constant prefactor), while the remaining term vanishes asymptotically for r sufficiently different 
from r' . A discussion of a similar result for the Fourier transform of the mf propagator given in Eq. ()114j) 
can be found in Ref. for the Bethe-Goldstone propagator. As a result, one preserves the integral 
equations for the wave function either in a partial wave basis as in Eq. (llOlj) or for the wave function 
in coordinate space as in Eq. (jlllj) in the case of mf propagators. The only difference with the free 
scattering case involves the use of the mf equivalents of the noninteracting propagators in coordinate 
space in Eqs. (I1U2|) and This result is due to the uniqueness of the on-shell momentum at a 

given energy which guarantees that the noninteracting wave function is a plane wave or spherical Bessel 
function (in a partial wave basis). One can therefore proceed with a similar asymptotic analysis as 
for free particles yielding a corresponding definition of the phase shifts as in Eq. ()1U5|) in terms of the 
on-shell scattering matrix. The result of Eq. ()106p also remains valid in this case. The presence of a 
nonvanishing phase shift is therefore continued to be linked to the nonvanishing of the imaginary part 
of the noninteracting propagator. In the case of mf scattering the corresponding energy domain resides 
above 2e{q = 0) which corresponds to the lowest energy of two hole states. 

These modest modifications of the quantities that characterize the scattering process as compared 
to the case of free-particle scattering, are related to the continued one-to-one relation of the energy 
with a unique relative momentum for which the noninteracting propagator has an imaginary part. This 
on-shell momentum emerges as the momentum that characterizes the plane wave (or spherical Bessel 
function) function describing the relative motion. The plane-wave character of the wave function allows 
for a conventional interpretation of the scattering process as in free space. 

Since the correlated wave function does not heal to the uncorrelated one when hole-hole propagation 
is included, the usual discussion of healing must be indeed be modified as indicated above. The standard 
interpretation of the validity of the shell model is couched in terms of the healing of the relative wave 
function to the noninteracting one. This interpretation was put forward in Ref. ^H] and is based on 
the use of the Bethe-Goldstone propagator in describing the effective interaction. Since this propagator 
excludes the propagation of two holes in Eq. ()114j) . the correction to the relative wave function in 
Eq. ()102j) due to the strong interaction potential heals within a characteristic healing distance to the 
spherical Bessel function provided that the energy is less than twice the Fermi energy. In scattering 
language this simply states that there is no phase shift for energies less than 2e^ when the Bethe- 
Goldstone propagator is employed since no corresponding imaginary part of the propagator exists |19j . 

An apparent contradiction then arises with this interpretation when one realizes that it is not permit- 
ted to neglect the propagation of the hole- hole term in Eq. ()114|) since it is essential for the understanding 




(114) 



n = 2£(go), 
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of the fragmentation of the sp strength below the Fermi energy Inclusion of hole-hole propagation 
yields a nonvanishing phase shift below 2ep |114j which is at odds with the healing interpretation of the 
relative wave function. On the other hand this contribution to the effective interaction is responsible for 
the presence of an imaginary part of the nucleon self-energy which is required in order to describe the 
experimental situation in nuclei as obtained from the (e,e'p) reaction TT. Evidently it is not possible 
to propagate mf nucleons which can generate a realistic sp strength distribution, and, at the same time, 
obtain the healing of the relative wave function to the noninteracting one which supposedly underlies 
the success of the mf picture. 

While recent (e,e'p) experiments have sharpened the range of the validity of the sp picture in terms of 
the more appropriate Landau quasiparticle description jHH] which is adequately described by microscopic 
theory jHH I113j . the paradox at the level of the relative wave function remains. A clue to the solution 
of this puzzle is provided by noting the inconsistency of the description of the sp strength in terms of a 
realistic spectral function and the construction of the effective interaction by means of a mf propagator. 
Clearly if the dressing effect of the nucleon is substantial and experiment |14| indicates it is, then one 
is forced to consider the construction of the effective interaction in terms of dressed nucleons. The 
consequences of this extension for the description of the scattering process in matter and the resolution 
of the healing paradox will be taken up in the following. 

The propagation of dressed nucleons requires a different treatment of the description of the scattering 
process. The main ingredient for this change is the form of the noninteracting propagator which is given 
in the medium by Eq. (jUHj) and for sp momenta k and k' reads 

k', n) ^ n^MkM',^') _ r,MkM',^i 

^ Je, Jep n-uj-uj' + tr] n-u-u'-iri 

The particle and hole spectral functions, Sp and Sh respectively, describe the distribution of the sp 
strength for a given momentum over the energy. These distributions are continuous and have sizable 
peaks either above or below the Fermi energy, corresponding to a momentum state above or below 
kp-i at the so-called quasiparticle energy. For kp = 1-36 fm~ , corresponding to normal density, the 
strength contained in the peak for momenta close to kp is typically only 70% [HniEHI- From the rest of 
the strength we refer for more details to the discussion in Sec. 14.21 

It should be noted that the noninteracting propagator in Eq. becomes the familiar mf Galitski- 
Feynman propagator (see Eq. pi4|) ) when mf spectral functions are inserted which are characterized by 
a 5 function peak of strength 1 at a sp energy either above the Fermi energy {k > kp) or below {k < kp). 
The difference between the Galitski-Feynman propagator and the dressed propagator is qualitatively 
different for the imaginary part and quantitatively for the real part. In Fig. |^ both the real and 
imaginary part of the dressed propagator (Eq. are shown as a function of the relative momentum 

for zero total momentum while the energy corresponds to an on-shell momentum go = 0.5 fm~ in the 
Galitski-Feynman case using a sp spectrum from Ref. jHIj. In the latter case the imaginary part of 
the propagator contains a S function only at a single momentum corresponding to 0.5 fm~^. It arises 
on account of the vanishing of the denominator in the hole-hole term at the on-shell momentum. The 
full curve in Fig. |^ corresponds to the imaginary part of the dressed propagator shown here as a 
function of momentum up to kp. There is a substantial spreading in the imaginary part containing 
even small high- momentum components (not shown in Fig. I?!]) . This spreading is a critical feature 
which completely alters the conventional picture of the scattering process. In Fig. |^ the real part of 
the dressed and mf propagator are compared for momenta below kp for an energy corresponding to an 
on-shell momentum of 0.8 fm~^. The dotted line corresponds to the mf propagator and also indicates 
the pole present at 0.8 fm~ . The dressed propagator shows a less dramatic momentum dependence 
and is in general substantially reduced from the mf propagator except for high momenta where both 
coincide. 
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Figure 21: Real (dotted) and imaginary part (full line) of the noninteracting propagator of 
dressed particles at an energy corresponding to an on-shell momentum of 0.5 fm~^ for the 
mf propagator. 




Figure 22: Comparison of the real part of the mf (dotted) and dressed (full line) propagator 
at an on-shell momentum of 0.8 fm~^ indicated by the vertical dotted line. 
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In order to perform the analysis of the scattering process it will be illustrative to use an analytical 
approximation to the noninteracting propagator of dressed particles (Eq. which contains the 

essential new features. In addition, only the case of zero total momentum of the propagating pair will 
be considered in the following. As a result, the noninteracting propagator contains equal and opposite 
momenta for the two particles (holes). Since the spectral functions do not depend on the direction of 
the momentum one can rewrite Eq. ()116p for the present purposes as 



' 7., 7„ n-uj-ij' + M, a-ui-^'-i,^ 

The momentum k not only corresponds to the absolute value of the sp momenta but also represents 
the relative momentum for the case of zero total momentum. 

Introducing a two-body self-energy term for purely practical reasons, one can attempt to write this 
noninteracting propagator as 

«f('=-") = n-En(,.o) - 

where the sign is determined by whether the energy Vt is above (+) or below (-) 2eF- By assuming that 
this ad hoc self-energy has a slowly varying imaginary part as a function of the relative momentum 
k one can expand the self-energy at the momentum k^ for which 

^] = Re S^^(A;o,^^). (119) 

Noting that the expansion is in the square of the momentum, one obtains a complex pole approximation 
(CPA) to the propagator by only keeping the real and imaginary part of at k^ and the first derivative 
of the real part. The resulting propagator has the form 



where the constant c is obtained from 



m I dk"^ 



(121) 



and 7 from 



7 = |lmS^^(A;o,^^) 



(122) 



Typical values of c at low energies correspond to 0.5 whereas it rises slowly to 1 for higher energies. 
This feature is closely related to the pattern of the distribution of the sp strength. The quasiparticle 
pole strength at kp = 1.36 fm~^ is about 0.7 (see Sec. I4.2|l so for a two-particle propagator close to 
these energies a factor of (0.7)^ is expected. For higher momenta the strength in the peak grows back 
to 1 yielding a propagator which is more of the mf or even free-particle kind. It is also apparent that 
this factor of about 0.5 can be identified from Fig. |221 It should be noted that the CPA is obtained 
after first numerically calculating the noninteracting propagator of the dressed particles. In Fig. 1221 the 
quality of this CPA to the propagator can be judged by comparing it to the numerically calculated 
result for the imaginary part of the propagator at an energy below 2ep. 

The CPA result for the propagator cannot be used to solve any of the integral equations for the 
effective interaction for example since it is only a good approximation to the noninteracting propagator 
close to the peak of the imaginary part. The full solution also requires an accurate representation of 
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Figure 23: Comparison for the imaginary part of the noninteracting propagator of dressed 
nucleons between the complete (numerical) result given by the full line and the simple CPA 
given by the dotted line for momenta below kp. 



the high-momentum components of the propagator in order to properly include the effect of short- 
range correlations in the interaction or wave function. The CPA result does provide a reasonable 
representation of the long-range part of the propagator and therefore can be profitably used to discuss 
the asymptotic analysis of the scattering process. Indeed, using the CPA to the dressed propagator one 
can repeat the steps involved in the Fourier-Bessel transform leading to Eq. ()99|) . For free or mf particles 
the integral in Eq. yields the product of a spherical Bessel function and one of the spherical Hankel 
functions with as argument the real on-shell momentum (see Eqs. fl98|) and ()115|) ). This on-shell 
momentum is real since the corresponding noninteracting propagators (Eqs. (j97jl and ()114|) ) can only 
have a vanishing denominator for a real momentum. Since Eq. ()95j) can be calculated by a contour 
integral for Eqs. (jHZI), ()114j) (at least for the long-range part), as well as for Eq. ()12()j) . it is clear that 
the presence of a nonvanishing imaginary part for the pole of Eq. fjl20|) . due to a the nonvanishing of 
7, will lead to a complex on-shell momentum which will be denoted by kq- Using the CPA (Eq. ()120j) ) 
for Q < 2eF one obtains from Eq. (j^H) 



9i!cPA{r,r';n) = -ic-^Koje{Kor<)h*f.[Kory^^ 



(123) 



The momentum argument of the spherical Bessel and Hankel functions, kq, is now complex, its real and 
imaginary part, Kq and Kq, respectively, are easily obtained from /cq and 7 (see Eqs. (I119|) and (jl22j) ) 
by determining the zeros of the denominator of Eq. (jl2(Jj) . Eq. (I123j) contains the Hankel function /i^ 
due to the different boundary condition associated with hole- hole propagation for energies below 2eF- 
This leads to a pole in the upper half of the complex fc-plane in contrast to the case of particle-particle 
or free-particle scattering. As a result, Kq is negative for fl < 2eF and its magnitude can become as 
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Figure 24: Comparison between the analytical approximation and the complete numerical 
result for the dressed noninteracting propagator in coordinate space. Shown is the imaginary 
part for an energy below 2eF for the CPA propagator (dotted line) and the complete result 
(full line). The corresponding propagators in momentum space are shown in Fig. |221 



large as 0.2 to 0.3 fm . The resulting propagator for i = can be written as (for r < r') 

glU,cPA{r. r'- n) = ^ . 124 

A comparison between this analytical result and the numerical Fourier-Bessel transform of the dressed 
noninteracting propagator which it approximates, is very successful in the domain where the propagator 
has its maximum as shown in Fig. |23 This is also true for the real part. This demonstrates that the 
propagator for dressed particles is radically different from the noninteracting or mf one (see Eq. (jHSl)) 
due to presence of damping terms related to the nonzero value of Kq. As noted before, there is no 
longer a unique on-shell momentum. Indeed, the complex pole at kq in the CPA propagator is just a 
simple (and approximate) representation of this feature. As a consequence, the relative wave function 
of the dressed particles contains a spread in momentum states. This, in turn, must yield a localization 
of the corresponding wave function in coordinate space. This is exhibited in the propagator Eq. ()124j) 
which represents the probability amplitude for removing a pair with relative distance r' and adding 
the pair after propagation at r (without interaction between the pair included yet). This amplitude 
peaks at r = r' (see Fig.l24p and is exponentially damped with the decay constant \kq\. This feature has 
interesting physical consequences since it means that if the separation distance between the scatterers is 
too large there is little probability that they will actually interact because this requires a small relative 
distance. Indeed, taking r' to much larger values than in Fig. |21] yields a negligible contribution to 
the noninteracting propagator near small r where the interaction will act to modify the wave function. 
Clearly this effect is governed by the size of the imaginary part of the pole of the CPA. It should be 
noted that this value will only approach zero when the scattering energy approaches 2eF- Just as in the 
case of the sp motion, this means that the noninteracting wave function will tend to a plane wave again 
only in this limit. The corresponding result for Fig. I^then yields a simple sine wave characterized 
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Figure 25: Decomposition of the numerical propagator in coordinate space into in- (dashed) 
and outgoing wave (dotted hne). Also shown is the sum of both contributions (full line). 



by Kq which approaches kp. For all other energies damping does occur sufficiently rapidly to warrant 
the following observation: since only the part of the wave which returns from the scattering can be 
affected and this part always decreases with increasing r, only that part of the noninteracting wave can 
be influenced by the scattering which is exponentially damped in r. This is illustrated in Fig. EHl where 
a numerical decomposition of the propagator is presented in terms of the incoming and outgoing wave. 
For values of r' outside the range of the interaction as in Fig. 123 this implies that even a substantial 
modiflcation of the outgoing wave will hardly affect the total propagator and the wave function must 
automatically heal, according to the equivaqlent value of Kq in CPA approximation, to the noninteracting 
one. 

The above observations allow for the resolution of the paradox related to the healing properties of 
wave functions in the medium. This property has been considered the physical justiflcation of the mf-like 
properties observed in nuclei in the presence of strong short-range interactions. There is overwhelming 
experimental evidence ^1] that sp motion in nuclei must be described in terms of dressed nucleons 
with substantial fragmentation of the strength. The original discussion of the healing properties of the 
relative wave function of particles in the medium jT^ used a Bethe-Goldstone propagator involving 
mf nucleons to arrive at the healing property of the relative wave function. The present discussion 
explores the consequences of scattering dressed particles and demonstrates that this dressing of nucleons 
automatically leads to a localization of the relative wave functions in coordinate space. The results for 
the CPA propagator analysis indicates that even with sizable phase shifts the localization of the wave 
function leads to the desired healing property of the wave function since the part of the wave function 
affected by the scattering event is exponentially damped. Also for the complete numerical propagator 
the same features are observed. Even in the presence of strong interaction processes the resulting 
picture of the nuclear medium is a tranquil one in which the dressed particles no longer remember 
their scattering event beyond some flnite distance and their wave functions heal to the corresponding 
noninteracting ones. This appears to be a satisfactory picture of a correlated medium in which particles 
do not carry the information of their interaction indefinitely around unlike a description of scattering 
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using a mf Galitski-Feynman propagator. 



4.5 In-medium cross sections and phase shifts 

For various reasons the study of the interaction between nucleons in the nuclear medium has retained 
a considerable amount of interest over the years. Of particular interest is the exploration of the "in- 
medium" interaction in the context of transport-theory descriptions of heavy-ion reactions jll7j . Typ- 
ical analyses simulate the dynamics of a heavy-ion reaction on the basis of kinetic equations like the 
Boltzmann-Uehling-Uhlenbeck (BUU) equation |118j . An essential ingredient in these BUU calcula- 
tions is the nucleon-nucleon cross section in the medium. Attempts to include medium-modified cross 
sections in such calculations have been described in Ref. jll9j . Calculations of the cross section be- 
tween nucleons in nuclear matter have been reported in Refs. |12Uj - [129\ . Some recent issues that 
have emerged from this work include the enhancement of the cross section at finite temperature due 
to the vicinity of a pairing instability [126J, the sensitivity of the cross section to the choice of the 
single-particle (sp) spectrum at zero temperature |127j . the density and energy dependence |128j . and 
temperature dependence of the cross sections |129j . 

All results obtained in Refs. [120j - |129j have been generated under the assumption that the sp 
motion of the nucleon in the medium is that of a mean-field (mf) particle. Under this assumption 
the scattering process in the medium takes place between nucleons which at most have a sp spectrum 
different from free space but are otherwise unaffected by the presence of other nucleons except for the 
Pauli principle related to a mf Fermi gas. This assumption has been contradicted unequivocally in 
recent years for finite nuclei by the careful analysis of the (e,e'p) reaction. It is therefore fair to say 
that both on the basis of experimental results as well as theoretical calculations for nuclear matter, it 
is prudent to consider the sp dressing of the nucleons in matter and the effect this has on the scattering 
process. One of the consequences of employing dressed particles is the localization of the two-body 
propagator in coordinate space as discussed in Sec. 14.41 severely limiting the range of the propagation 
for most energies. This feature is due to the presence, at a given energy, of a range of momenta which 
determine the relative wave function of the propagating particles. While for mf particles a unique (on- 
shell) momentum characterizes the relative wave function which therefore corresponds to a plane wave 
(or spherical wave), the presence of different momentum components implies that the determination of 
phase shifts and cross sections requires some kind of folding procedure over these momenta in the case 
of dressed particles. In Ref. jHl] a set of expressions have been introduced to characterize the scattering 
event of dressed particles involving such a folding procedure. Results for phase shifts and cross sections 
obtained from these expressions can therefore be compared with the corresponding results propagating 
free and mf particles using the Reid soft core potential [S^ . 

As proposed in Ref. |H1I and followed up in Ref. [HE] this folding of the effective interaction takes 
the following form for the S'-matrix element in the case of coupled channels 



This result reduces to the conventional results for free or mf particles. In the case of coupled channels 
Eq. (jl25|) can be used to follow the procedure to obtain phase shifts by diagonalization of the S'-matrix. 
For an uncoupled channel, the matrix element in Eq. ()125|) directly yields the phase shift in the form 



A consequence of the present approximation is that the phase shifts S'^^'^ remain real As a result, 
the phase shifts can be fruitfully compared with results for mf or free particles. Equation ()125|) is exact 




(125) 




(126) 



49 



for noninteracting or mf particles and for dressed particles includes the physically reasonable expectation 
that the distribution over the momenta as contained in the imaginary part of the propagator will feature 
in determining the scattering process. An example of such a distribution was shown in Figs. |^ and 
1221 While this approximation does not make sense at large distance scales as discussed in the previous 
section, it provides, locally, a very reasonable generalization of the phase shift. The corresponding 
"short-distance" approximation to the scattering amplitude yields the following result 

= 47r^ ^ /(-z)^y,^,(f)F;^,(i)(£m,5m,|JM)(fm^5m;|JM) 

U'J mm'M 

POO 

/ dkk lm{gy{k; fi)} {k{iS)J\T{n)\k{i'S)J), (127) 
Jo 

where a coupling to total spin S and projections ms,m'g for initial and final spin states has been 
included together with the usual decomposition in partial waves. For the total cross section (in the 
neutron-proton case) one obtains 

(Ttot = vr ^ (2 J + 1) 

SU'J 

Eq. (|128|) demonstrates that a sensible cross section will be obtained in the case of dressed particles 
at all energies for which a nonvanishing imaginary part of the propagator exists. For two particles 
deep in the Fermi sea, for example, Eq. p28|) avoids the divergence associated with the k^"^ term in 
the corresponding expression for free or mf particles. The formulation of the cross section in terms of 
Eq. ()128|) provides a reasonable way to assess the strength of the interaction between dressed particles 
in the medium in terms of the square of the relevant transition matrix element (F) multiplied by an 
appropriate measure of the density of states represented by the imaginary part of the noninteracting 
propagator. 

The main results for the phase shifts for some of the more important partial wave channels are 
summarized in Fig. 1221 A comparison is made between phase shifts for free particles (solid line), mf 
particles at kp = 1.36 fm~^ (dashed line), and dressed particles (short-dashed line) at the same density 
for the ^5*0, '^S'l, ^Pi, and ^Di channels (corresponding to the different panels in Fig. 1221) a function of 
the on-shell momentum. In general one finds that the dressed phase shifts suggest weaker interactions 
since in essentially all cases they are either less repulsive or less attractive than the mf result. By 
studying the individual contributions of the real and imaginary part of the effective interaction one may 
gain more insight into this issue (see below). For the two S'-wave channels the most striking feature 
of the dressed phase shift is the disappearance of the pairing signature for the ^5*0 channel and the 
enormous reduction of the signal in the ^5*1 case. While the dressed ^5*0 phase shift is essentially zero 
at kp, it is still clearly attractive at this momentum for the ^Si channel. The actual calculation of 
the phase shift for this channel displays a slight kink close to kp suggesting that the phase shift may 
actually rise very rapidly to vr very close to kp. This implies a tremendous reduction in the strength of 
the pairing correlations in this coupled channel as compared to a mf treatment. Gaps of the order of 10 
MeV have been obtained for this channel in Refs. jHSl 11301 1131j . Clearly, the dressing of the nucleons 
has a strong influence on pairing. While one would expect to obtain a gap using dressed nucleons 
based on the attractive effective interaction at the Fermi surface, its magnitude is drastically reduced 
as suggested by the phase shift calculation shown in Fig. (20 This latter feature has recently been 
confirmed in Ref. [^. The main ingredient in this reduction is the decrease in the density of states at 
2ep when dressed nucleons are propagated. This reduction is essentially the square of the strength of 
the quasiparticle pole at kp leading to a reduction factor of about 0.5. Since pairing correlations are 
particularly sensitive to this density of states, it is not surprising that the strength of the pairing is 
substantially diminished when dressing is taken into account. It is also noteworthy that one observes a 
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Figure 26: Comparison of phase shifts for free particles (sohd), mf particles (dashed), and 
dressed particles (short-dashed lines) for different partial waves. The density of the medium 
corresponds to kp — 1-36 fm~^. 
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Figure 27: Total neutron-proton (top) and neutron-neutron (bottom panel) cross sections 
for free (solid), mf (dashed), and dressed particles (short-dashed line) as a function of the 
on-shell momentum. 



smaller negative phase shift for both S*- waves at higher energy as compared to the mf result. A similar 
conclusion may be drawn by inspecting the phase shifts for the ^Pi and ^Di channel in the bottom panels 
of Fig. ESI Also for these partial waves which represent repulsive effective interactions, one observes 
a reduction of the magnitude of the phase shift when dressing is considered. It is also important to 
note that in the case of mf propagation the results in Fig. ISEl show that the corresponding results tend 
to those of free particles at high energy, whereas this is not the case for dressed particles. This latter 
result indicates that the effect of the dressing extends to a large energy domain. This observation is 
not too surprising since the spreading of the sp strength due to short-range and tensor correlations 
takes place in a very large energy domain 60, 58, 61j and is quite different from a local (in energy) 
spreading of the strength as would be obtained by a complex quasiparticle energy. The results for the 
phase shifts obtained for dressed particles lead to the expectation that the corresponding total cross 
sections are substantially reduced compared to the mf results. The results for neutron-proton (np) and 
neutron-neutron (nn) total cross sections displayed in Fig. EZI confirm this expectation. Results have 
been obtained for free (solid), mf (dashed), and dressed particles (short-dashed line) by including all 
partial wave channels of the Reid potential with J < 2. 

Results for mf particles were generated with a realistic sp energy spectrum and are similar to the 
corresponding results obtained e.g. in Ref. |127j . The effect of the pairing correlations on the cross 
sections yields a cusp-like behavior around kp reminiscent of the enhancement of the cross sections 
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obtained by the Rostock group at finite temperature jll9t I126j . As the phase shifts for mf particles 
suggest, the corresponding cross sections in the medium become essentially identical to the one in free 
space at high energy. Both for the np and nn total cross sections the effect of dressing the nucleons 
is quite dramatic leading to a substantial reduction of the total cross section at all energies. Indeed, 
on average a cross section of only about 10 mb is obtained. While this may seem a small number, 
it should be kept in mind that this by no means implies that the effective interaction in the medium 
has become insignificant In addition, one should recall that the concept of asymptotic flux in the 
medium representing preserved information of a scattering event deep in the medium is not a realistic 
consideration when the dressing of the nucleons is significant [HII, as it is at kp = 1.36 fm^^. The 
main ingredient representing the dressing is the two-particle density of states, its reduction for dressed 
particles is to a large extent responsible for the reduction of the cross section. In this respect it should 
be noted that the results of Ref. ||126j also show a substantial reduction of the total cross section at high 
energy with increasing temperature. This result implies that with increasing temperature which also 
means a larger depletion of the Fermi sea due to thermal excitations, one obtains a reduction of the total 
cross section similar to the one obtained here at zero temperature when the dressing of the particles is 
incorporated (and therefore the depletion of the Fermi sea due to correlations is included). While no 
results are shown in Fig. EH below 0.5 fm~^ in order to avoid the large value of the total cross sections 
for free particles, it should be noted that the cross sections for dressed particles smoothly go to zero 
when expression (jl28p is used at lower energies. This expression avoids the problem associated with the 
expression for free or mf particles which would yield an infinite cross section for the on-shell momentum 
going to zero. In addition, the latter expression does not yield a cross section for energies that do not 
yield a solution for the on-shell momentum, i. e. for energies deep in the hole-hole continuum. 

A related quantity to the two-body interaction is the two-body propagator given by Eq. in 
momentum space and Eq. (j94j) in coordinate space. Since a Lehmann representation exists for this 
two-particle propagator, it is possible to define a spectral density for this quantity either in momentum 
or coordinate space. We will consider here 



which represents the removal or addition amplitude for a pair of particle (for zero total momentum) 
depending on whether the energy Q is below or above 2eF, respectively. In Fig.EHlthis quantity is plotted 
for r = r' in the case of the ^5*0 partial wave for an energy corresponding to —100 MeV. This diagonal 
form of Eq. ()129|) corresponds to the two-particle spectral function in this chosen representation. The 
curves in this figure correspond to the full propagator (labeled S^^), the contribution of the first term 
from Eq. representing two dressed but noninteracting nucleons (Sj^^^), and the contribution of the 
second term of Eq. (j93p representing the vertex function F (Sj^). In addition, the result for two mf 
particles is shown by the curve labeled Sq^ . Notice that Sj^^^ has a large positive value at the origin and 
Sj^ has a large negative value at the origin. There is a precise cancellation which is illustrated in Fig. |2H1 
This cancellation is of course associated with the effect of the short-range core of the Reid interaction 
which suppresses the relative wave function at these disctances. Integrating Eq. ()129p over the allowed 
two-particle removal energies yields the corresponding contribution to the two-particle density matrix 



This quantity also characterizes the behavior at short-range for a chosen partial wave and can be 
compared to the corresponding quantity using dressed but noninteracting propagators. It is attractive 
to cast this comparison in the form of the ratio 
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Figure 28: Two-nucleon spectral function for the ^Sq channel in nuclear matter for zero total 
momentum and a removal energy corresponding to —100 MeV. 



In Fig. 1221 we plot the correlation function calculated for ^5*0 nucleons with center of mass momen- 
tum zero. Note that this correlation function shows the characteristic depression at r = and an 
enhancement somewhat above r = 1 fm. This correlation function has been employed in the analysis 
of two-proton removal experiments as discussed in Sec. 15.71 In Ref. |132j a related study has been 
reported for the probability for the removal of pairs with the actual properties of the deuteron yielding 
information about the deuteron distribution in nuclear matter. 

4.6 Results for self-consistent Green's functions in Nuclear Matter 

The solution to the scheme outlined in Sec. 14.11 has been studied by different groups. Aspects of self- 
consistency have been treated in Refs. [HHI - [ZH] • Complete self-consistency has been implemented for a 
realistic interaction in terms of a discrete representation of the propagator in Ref. [7S]. A representation 
of the spectral functions on a momentum and energy grid was made self-consistent in Re. [7B] for a 
separable version of the Paris potential. Both of these strategies work better when softer interactions 
are used. Indeed, it appears that yet another method is required to deal with interactions like the Reid 
potential jSH]- The main reasons for potential problems is the need to consider a very large energy 
domain since the peaks in the imaginary part of the self-energy and effective interaction occur around 
10 GeV [nni IM] for a potential like Reid's. This feature puts stringent demands on the sizes of grids to 
describe the spectral strength as has been found in Ref. 170] . 

The essential difficulty in the self-consistency procedure is therefore the handling of the information 
which describes the complete energy and momentum dependence of the spectral functions over the 
relevant domains. For interactions like the Reid potential the energy and momentum range needed 
to store the relevant information precludes a straightforward numerical discretization of the spectral 
functions. For this reason earlier attempts by the St. Louis group have involved a parametrization 
of the spectral functions themselves |64,, 65J. Although manageable for the calculation of the effective 
interaction at zero total momentum ;66J, it is quite difficult to manage the strong energy dependence 
of these functions in a complete self-consistency scheme. In particular it is necessary to represent the 
spectral strength by at least two functions, one of which maybe sharply peaked. This leads to some 
difficulties in generating the real parts of the propagators from dispersion relations. As a result, the 
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Figure 29: Correlation function for ^5*0 nucleons with total momentum zero. 



self-energy that can be obtained from the propagator by inverting the Dyson equation does not behave 
sufficiently smoothly at all energies for a given momentum. 

We therefore outline here another method that is able to deal with such interactions as the Reid 
potential [721 EH] implying that it will certainly work for softer interactions. First, it is noted that the 
sp spectral function is completely determined by the Im S. The spectral function can be written in 
terms of the self-energy as in Eqs. (j73|) and (j74|) . The self-energy can be written according to Eq. (j64j) 
which demonstrates that the spectral function can be completely expressed in terms of Im Il{k,uj), 
Re ASj^(/c,a;), Re AE|(/c, tu), and T,Y{k). The contributions Re AS|(/c,ci;) and Re AT,-^{k,uj) can be 
obtained from the Im E(A;,u;) by performing the dispersion integrals given in Eq. ()64a|) . The corre- 
lated Hartree-Fock contribution, Sy(A;), can be determined in a given iteration step as follows. The 
quasi-particle (qp) spectrum usually identifies the location of the peak of the spectral function and 
is determined by Eq. ()82|1 . The momentum distribution given by Eq. 1)6 HI from a previous iteration, 
say the {i — l)th one, can be used to determine the correlated Hartree-Fock self-energy from Eq. ()60|) . 
This allows the construction of a new Fermi energy and allows the calculation of the ith iteration for 
the spectral functions closing the self-consistency loop. This analysis indicates that it is sufficient to 
accurately represent the imaginary part of the self-energy for a complete determination of the spectral 
functions. Since this imaginary part is very well-behaved as a function of energy and changes also 
smoothly as a function of energy, a strategy has been developed to represent Im S in terms of a limited 
set of gaussians j73 EH] • Two gaussians below the Fermi energy and three above appear to provide 
sufficient flexibility for a complete representation of the imaginary part of the self-energy. 

Some results of this type of self-consistent calculation will now be presented here. We start by 
discussing the results for the momentum distribution at a density corresponding to kp = 1.36 fm~ . 
The corresponding result is shown in Fig. EOj This flgure should be compared to the first generation 
results shown in Fig. for a calculation that can be considered as a first iteration step in this self- 
consistency procedure since it is based on mf propagators as input for the spectral functions. One notices 
immediately the similarity of the occupation of momenta below kp. Indeed, only a slight increase in the 
occupation of at most 3% is observed. This implies that the corresponding depletion due to short-range 
correlations is still about 15% for nucleons deep in the Fermi sea. The slightly increased occupation 
below kp results in a slight reduction of momentum states above kp- The occupation of momenta just 
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Figure 30: Momentum distribution in nuclear matter obtained for self-consistent Green's 
functions using the Reid potential. 



above kp, however, hardly changes. This is further corroborated by the result for the quasiparticle 
pole at the Fermi momentum since it changes from 0.72 for the first-generation result to 0.75 in the 
self-consistent case. 

Some results for spectral functions are shown in Fig. |^ around the Fermi energy as a function of 
u! — ep- A comparison with the results reviewed in Sec. 14.21 vields a striking difference for the strength 
distribution below the Fermi energy. The spectral strength in Fig. |^ is essentially identical at large 
negative energies for all relevant momenta. This feature reflects an essentially similar energy distribution 
of the imaginary part of the self-energy for these momenta. This could have been expected since the 
constraint imposed by mf sp energies no longer applies and all momenta have similar energy domains 
associated with the imaginary part of the self-energy. In the latter case, this imaginary part of the 
self-energy has a fixed lower bound depending on the momentum as discussed in Sec. 14.21 and shown 
in Fig. ^1 The appearance of sp strength at large negative energies has important implications for 
the binding energy of nuclear matter as will be discussed shortly. The results for spectral functions at 
high energy are shown in Fig. |221 One continues to encounter a common tail of the sp strength at high 
energies above ep for all momenta which is quantitatively similar to the one shown in Fig. ^Jfor the 
first-generation result. 

The special role of short-range correlations in obtaining saturating behavior of nuclear matter is 
illustrated in Fig. EHl In this figure we plot the integrand corresponding to both terms in Eq. (fTSj) as 
a function of momentum after performing the energy integral over the spectral function for densities 
corresponding to kp = 1.36 fm~^ and kp = 1.45 fm~^, respectively. At kp = 1.36 fm~^ the high- 
momentum components still provide attractive contributions whereas for kp = 1.45 fm~^ a changeover 
occurs suggesting that at an even higher density these high-momentum terms will provide only repulsion. 
From this analysis it is clear that the expected relevance of SRC in obtaining reasonable saturation 
properties of nuclear matter is fully confirmed. It remains to relate this observation to the vast body 
of work on the nuclear-matter saturation problem. 
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Figure 31: Self-consistent spectral functions for three different momenta at kp = 1-36 fm 
corresponding to (full), 1.36 (dotted), and 2.1 fm""*^ (dashed) as a function of a; — Sp. 
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Figure 32: Self-consistent spectral functions for three different momenta at kp = 1.36 fm~^ 
corresponding to (full), 1.36 (dashed), and 2.1 fm~^ (dotted) but plotted on a logarithmic 
energy scale as a function of lu — Sp. 
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Figure 33: The high-momentum contribution to the energy per particle for kp = 1.36 fm^^ 
(sohd) and 1.45 fm~^ (dashed). This result illustrates the source of the saturation process 
when short-range correlations are considered self-consistently. 



4.7 Saturation of nuclear matter from short-range correlations 

A correct description of the saturation properties of nuclear matter has remained an unresolved issue 
for a very long time. The Brueckner-Bethe-Goldstone (BBG) expansion jTHj supplies a converged result 
for the energy per particle in the relevant density range, for a given realistic interaction, at the level 
of three hole-line contributions |133| I134j . Such calculations fail to reproduce the empirical saturation 
properties which require a minimum in the equation of state at a density corresponding to a Fermi 
momentum, kp, of about 1.33 fm~^ with a binding energy of about 16 MeV. The authors of Ref. |134j 
obtain for the Argonne f 14 interaction |135j a saturation density corresponding to 1.565 fm~^ with about 
the correct amount of binding. This corresponds to an overestimation of the empirical density by about 
60% but appears completely consistent with corresponding variational calculations |136j for the same 
interaction. 

Several different remedies for this serious problem have been proposed over the years. The intrin- 
sic structure of the nucleon and its related strong coupling to the A-isobar inevitably requires the 
consideration of three-body (or more-body) forces. When three-body forces are considered in varia- 
tional calculations it is possible to achieve better saturation properties only when an adhoc repulsive 
short-range component of this three-body force is added |137t 1138] . It has also been suggested that 
a relativistic treatment of the nucleon in the medium using a Dirac-Brueckner approach provides the 
necessary ingredients for a better description of saturation |139j - |142j . 

All many-body methods developed for nuclear matter have focused on a proper treatment of short- 
range correlations (SRC) without the benefit of experimental information on the influence of these 
correlations on the properties of the nucleon in the medium. This influence can now be clearly identified 
by considering recent results from (e,e'p) reactions ^T] - ^3] and theoretical calculations of the nucleon 
spectral function in nuclear matter jSH] - |SHl ISO] • 
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A recent analysis of the (e,e'p) reaction on ^osp]^ ^ wide range of missing energies and for missing 
momenta below 270 MeV/ c yields information on the occupation numbers of all the deeply-bound proton 
orbitals (see also Sec. l5.9|) . These data indicate that all these orbitals are depleted by the same amount of 
about 15% |143tll44j . These occupation numbers are associated with the orbits which yield an accurate 
fit to the (e,e'p) cross section. The properties of these occupation numbers suggest that the main effect 
of the global depletion of these mean-field orbitals is due to SRC. Indeed, the effect of the coupling of 
hole states to low-lying collective excitations only affects occupation numbers of states in the immediate 
vicinity of the Fermi energy |145j . In addition, nuclear matter momentum distributions display such an 
overall global depletion due to short-range and tensor correlations [SIl IHH 110] as discussed in Sec. 14.21 
and confirmed in Sec. 14.61 for the fully self-consistent results. The latter results formed the basis of the 
now corroborated prediction |3H I146| I147j for the occupation numbers in ^osp^-, \14:3\ . 

Most of this depleted sp strength is located at energies more than 100 MeV above the Fermi en- 
ergy |60t 1^ 153]. This appearance of strength at high energy is another important aspect of the 
influence of short-range and tensor correlations. Yet another characteristic feature of these SRC is that 
this depletion of the sp strength must be compensated by the admixture of a corresponding number 
of particles with high-momentum components. These high-momentum components have not yet been 
unambiguously identified but are currently studied experimentally [148^ as discussed in Sec. 15.91 Sohd 
theoretical arguments |149j and calculations clearly pinpoint this strength at high excitation energy in 
the hole spectrum both for nuclear matter [20] and finite nuclei |150j (see Sec. 15. 3|) . Indeed, experiment 
confirms that no substantial admixture of these high-momentum components is observed in the vicinity 
of the Fermi energy |151j as discussed in Sec. 15.21 

We now present an argument showing that short-range correlations are the dominant factor in 
determining the empirical saturation density of nuclear matter. We recall that elastic electron scattering 
from 208p]^ |l52j accurately determines the value of the central charge density in this nucleus. By 
multiplying this number by A/Z one obtains the relevant central density of heavy nuclei, corresponding 
to 0.16 nucleons/fm^ or kp = 1.33 fm~^. Since the presence of nucleons at the center of a heavy nucleus 
is confined to s-wave nucleons, and, as discussed above, their depletion is dominated by SRC, one may 
therefore conclude that the same is true for the actual value of the empirical saturation density of 
nuclear matter. While this argument is particularly appropriate for the deeply bound lsi/2 and 2si/2 
protons, it continues to hold to a large extent for the 3si/2 protons which are depleted predominantly 
by short-range effects (up to 15% ) and by at most 10% due to long-range correlations [121 112] • These 
considerations demonstrate clearly that one may expect SRC to have a decisive infiuence on the actual 
value of the nuclear-matter saturation density. 

High-momentum components due to SRC also have a considerable impact on the binding energy 
of nuclear matter. This result can be inferred from the energy sum rule given in Eq. ()75p. Eq. ()75|) 
illustrates the link between the energy of the system and the hole spectral function, Sh{k, uj) illustrated in 
Fig. ESI Results for the momentum distribution and true potential energy based on the spectral function 
show that enhancements as large as 200% for the kinetic and potential energy over the mean-field values 
can be obtained for both nuclear matter |61j and finite nuclei |150j (see Sec. 15. 3p . These large attractive 
contributions to the potential energy of nuclear matter are mainly from weighting the high-momentum 
components in the spectral function with large negative energies in Eq. (f7S|) . The location of these high- 
momentum components as a function of energy is therefore an important ingredient in the determination 
of the energy per particle as a function of density. So far, the determination of this location has relied 
only on quasiparticle properties in the construction of the self-energy. A self-consistent determination of 
the spectral function including the location of these high-momentum components therefore includes the 
dominant physics of SRC in the description of nuclear matter and is consistent with the experimental 
observations of the nucleon spectral function in nuclei. 

Such a determination requires the solution of the ladder equation for the effective interaction in the 
medium as discussed above. The implementation of this self-consistency scheme is numerically quite 
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Figure 34: Nuclear matter saturation points calculated with various realistic NN interac- 
tions. The open symbols refer to continuous choice Brueckner-Hartree-Fock results. The 
filled symbols refer to self-consistent results and represent saturation points calculated in 
the discrete scheme, except for the old Reid (Reid68) interaction where the binding energy 
at two densities is shown in the continuous scheme. 

involved and has been attempted by several groups [HE] - [ZH] as discussed above. As reported in Ref. [H2] 
two different approaches have been used to generate results for different interactions. In the continuous 
scheme, a representation of the imaginary part of the self-energy in terms of four gaussians is used to 
completely describe the imaginary part of the self-energy and thereby the sp propagator as discussed 
in the previous section. The parameters of these gaussians are then determined self-consistently jTH] 
for the Reid potential ISH] by solving Eq. ()72j) with the convolution of spectral functions in Eq. ()68|) as 
input, constructing the self-energy, and then solving the Dyson equation (jHSj)- In the discrete scheme a 
representation of the propagator in terms of three discrete poles [ZSIIZSI is used avoiding a full continuum 
solution of Eq. (f?^ . The latter approach is equivalent to a continuous version as far as the energy per 
particle is concerned, since it requires a reproduction of the relevant energy-weighted moments of the 
hole and particle spectral function [731 ESI- This is substantiated by comparing the results of this 
discrete scheme with the results of the continuous self-consistency scheme used in Ref. jHIj for the 
Mongan-type separable interaction |153j and in Ref. [77] for the separable Paris interaction [SUl IHT| . 
It is found that the binding energies correspond to within 5% over the relevant kp range around the 
minimum, and moreover that the location of the minimum agrees to within 3%. 

In Fig. 123] the saturation points obtained within the discrete scheme of Ref. |731 [7S] for the updated 
Reid potential (Reid93), the Nijml and Nijmll interaction [THj and the separable Paris interaction 
[SUl IHT| are shown. The results demonstrate an important and systematic change of the saturation 
properties with respect to continuous choice Brueckner-Hartree-Fock (ccBHF) calculations, leading to 
about 4-6 MeV less binding, and reduced values of the saturation density, closer to the empirical one. 
Such a trend is entirely consistent with the observations in Ref. [77| made for a separable NN interaction, 
and is now extended to more realistic (non-separable) interactions. 

As discussed in the previous section, the discrete scheme [7^1 [THj could not be used with the original 
Reid (Reid68) potential because of its slow decay in momentum space, but some results are available 
in the continuous scheme of Ref. [7^ 175] . In Fig. [^ the binding energy is shown at two densities 
{kp = 1.33 and 1.45 fm~^); the error bars are an estimate of the remaining uncertainty due to incomplete 
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convergence and the non-selfconsistent treatment of some higher order partial waves [7^ l78j. The results 
again seem to indicate a substantial shift in the saturation density for the Reid68 potential, from the 
ccBHF value of about 1.6 fm~^, to a value below 1.45 fm~^, without seriously underbinding nuclear 
matter. 

The present self-consistent treatment of SRC (scSRC) differs in two main aspects from the ccBHF 
approach, the latter giving saturation results essentially equivalent to converged three hole-line calcu- 
lations |134j . Firstly, hole and particle lines are treated on an equal footing, thereby ensuring thermo- 
dynamic consistency [67 . Intermediate hole-hole propagation in the ladder diagrams is included to all 
orders. This feature provides, compared to ccBHF, a substantial repulsive effect in the k < kp contri- 
bution to Eq. (frSj) . and comes primarily from an upward shift of the quasi-particle energy spectrum as 
a result of including ut < ep contributions to the imaginary part of the self-energy. The effect increases 
with density, and is the dominant factor in the observed shift of the saturation point. Secondly, the 
realistic spectral functions, generated from the self-consistent procedure outlined in Sec. 14.11 and used 
in the evaluation of the in- medium interaction F and self-energy E, are in agreement with experimental 
information obtained from (e,e'p) reactions. For the Reid93 interaction at kp = 1.37 fm~^ we find 
z = 0.74 for the quasiparticle strength at the Fermi momentum, whereas the hole strength for k = 0, 
integrated up to 100 MeV missing energy, equals 83%; similar values are found for the other interactions 
and confirmed by the results discussed in the previous section for the Reid potential. The depletion 
of the quasiparticle peaks is primarily important to suppress unrealistically large pairing instabilities 
around normal density. The improved treatment of the high-momentum components does affect the 
binding energy, through the k > kp contribution to Eq. (f73j) . This feature, studied in [731 ES], provides 
a sizeable attraction, but is smaller than the afore-mentioned repulsive effect. 

The inclusion of /i/i-propagation in scSRC also leads to a somewhat stiffer equation-of-state than in 
ccBHF. A recent analysis of the giant monopole resonance in heavy nuclei |154j yields an experimental 
estimate Knm = 210 ± 30MeV for the nuclear matter compression modulus, 

(132) 

At the saturation points in Fig. |^ we find ccBHF values Knm = 154 MeV for Reid93 and Knm = 148 
MeV for the separable Paris interaction, which are enhanced to Knm = 177 MeV and Knm = 216 
MeV, respectively, in our scSRC calculation. These values agree reasonably well with the experimental 
estimate. Note that reasonable values for Knm imply that the Reid68 energies in Fig. ED may still 
deviate by 1-1.5 MeV from numerically exact scSRC values, as indicated by the error bars. 

The present results indicate that a sophisticated treatment of SRC lowers the ccBHF saturation 
densities, bringing them closer to the empirical one. It remains to be understood why apparently con- 
verged hole-line calculations |134j yield higher saturation densities. The three hole-line terms obtained 
in Ref. |134j indicate reasonable convergence properties compared to the two hole-line contribution. One 
may therefore assume that these results provide an accurate representation of the energy per particle 
of nuclear matter as a function of density for the case of nonrelativistic nucleons and two-body forces. 
At this point it is useful to identify an underlying assumption when the nuclear-matter problem is 
posed [751 1Z2] • This assumption asserts that the influence of long-range correlations in finite nuclei and 
nuclear matter are commensurate. It has been suggested in Ref. that this underlying assumption 
is questionable. 

Three hole-line contributions include a third-order ring diagram characteristic of long-range correla- 
tions. The effect of long-range correlations on nuclear saturation properties is sizeable, as shown by the 
results for three- and four-body ring diagrams calculated in Ref. [155j (see also Refs. |133| 1134] ). The 
agreement of three hole-line calculations with advanced variational calculations jl36j further supports 
the notion that important aspects of long-range correlations are included in both these calculations. 
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Figure 35: The curves labeled 2N, 3N, and 4N correspond to the 2-, 3-, and 4-hole-line 
results from Ref. |133j for the Reid interaction. Added to the 4N results are the contribution 
of three- and four-body ring diagrams involving at least one A-isobar configuration. These 
terms can be considered as three- and four-body force contributions when only nucleons are 
considered are therefore labeled with 3N- and 4N-force, respectively. 



This conclusion can also be based on the observation that hypernetted chain calculations effectively in- 
clude ring-diagram contributions to the energy per particle although averaged over the Fermi sea |156j . 
The effect of these long-range correlations on nuclear saturation properties is not small and can be 
quantified by quoting explicit results for three- and four-body ring diagrams |155j . These results for the 
Reid potential, including only nucleons, demonstrate that such ring-diagram terms are dominated by 
attractive contributions involving pion quantum numbers propagating around the rings. Furthermore, 
these contributions increase in importance with increasing density. Including the possibility of the cou- 
pling of these pionic excitation modes to A-hole states in these ring diagrams leads to an additional large 
increase in the binding with increasing density |155j . This result is illustrated in Fig. where these 
contributions are added to the hole-line expansion results from Ref. |133j . Alternatively, these terms 
involving A-isobars can also be considered as contributions due to three- and four-body forces in the 
space of only nucleons. These A-isobar terms have a very different behavior from the phenomenological 
three-body forces that have been used to describe the properties of nuclear matter 1X37111381 . The latter 
three-body forces turn repulsive at higher density. The importance of the pionic long-range contribu- 
tions to the binding energy shown in Fig. 1221 is related to the possible appearance of pion condensation 
at higher nuclear density |157t I158j . These long-range pion-exchange dominated contributions to the 
binding energy appear because of conservation of momentum in nuclear matter. For a given momentum 
q carried by a pion around a ring diagram, one is able to sample coherently the attractive interaction 
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that exists for values of q above 0.7 fm~ . All ring diagrams contribute coherently when the interaction 
is attractive and one may therefore obtain huge contributions at higher densities which reflect the im- 
portance of this collective pion-propagation mode |159j . This is clearly illustrated by the results shown 
in Fig. ESI 

No such collective pion-degrees of freedom are actually observed in finite nuclei. A substantial part 
of the explanation of this fact is provided by the observation that in finite nuclei both the attractive and 
repulsive parts of the pion-exchange interaction are sampled before a build-up of long-range correlations 
can be achieved. Since these contributions very nearly cancel each other, which is further facilitated 
by the increased relevance of exchange terms |16Uj . one does not see any marked effect on pion-like 
excited states in nuclei associated with long-range pion degrees of freedom even when A-hole states 
are included |161j . It seems therefore reasonable to call into question the relevance of these coherent 
long-range pion-exchange contributions to the binding energy per particle since their behavior is so 
markedly different in finite and infinite systems. One may consider the salient difference of the ratio of 
spin-longitudinal and spin-transverse response function in nuclear matter and finite nuclei as another 
indication of the relevance of this suggestion |162j . We also like to point out that experimental infor- 
mation of these response functions |163j - |165j suggests no characteristic enhancement of the (pionic) 
spin-longitudinal response as expected on the basis of nuclear matter calculations. 

Clearly, the assertion that long-range pion-exchange contributions to the energy per particle need 
not be considered in explaining nuclear saturation properties, needs to be further investigated. In 
practice, this means that one must establish whether pion-exchange in heavy nuclei already mimics the 
corresponding process in nuclear matter. If this does not turn out to be the case, the arguments for 
considering the nuclear-matter saturation problem only on the basis of the contribution of SRC will 
be strenghtened considerably. Furthermore, one would then also expect that the contribution of three- 
body forces |l(j6j to the binding energy per particle in finite nuclei continues to be slightly attractive 
when particle number is increased substantially beyond ten particles |167j . This point and the previous 
discussion also suggest that there would be no further need for the ad-hoc repulsion added to three-body 
forces used to fit nuclear-matter saturation properties |137j . 

5 Theoretical calculations for finite nuclei 

The structure of finite nuclei presents different features from the case of nuclear matter due to the fact 
that nucleons are localized inside the system. This requires the consideration of single particle basis 
states that have good total angular momentum (and parity) and correspond to discrete set of levels, 
as in the IPM. The coupling to low-lying collective excitations in the system tends to fragment the sp 
strength of these orbitals and to spread it over a wide range of energies associated with the energy scales 
of these excitations. Along with these effects, short-range correlations can move a sizable fraction of 
nucleons to high-momentum states as in the case of nuclear matter discussed in the previous section. 

The influence of SRC in finite nuclei has been studied theoretically by calculating the momentum 
distribution in the ground state of selected nuclei |168j - |173j . These results clearly show that for 
momenta above 400 MeV/c short-range and tensor correlations completely dominate the momentum 
distribution. Although the momentum distribution is a sp quantity, and thus allows only an indirect 
measurement of correlations, a careful study can give some indications on the importance of these effects. 
More information is contained in the hole spectral function Sejr{k,E), which gives the distribution of 
nucleons both in momentum and energy for particles with orbital angular momentum i, total angular 
momentum j, and isospin r. In nuclei, substantial correlations beyond the IPM approach are present 
and a full understanding of the many-body system requires, both in experimental as well as theoretical 
studies, the knowledge of S£jr{k, E) for a large range of energies E as will be discussed in Sec. 15.21 

Microscopic calculations of the spectral function and the momentum distribution of finite nuclei 
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based on realistic nuclear hamiltonians have been performed for very light nuclei {A < 4) in Refs. |174j 
- |177j . Results for heavier nuclei are typically derived from investigations of nuclear matter assuming 
a local density approximation (LDA) |172l 1178] - |18Uj . Other variational calculations for ^^O yield the 
momentum distribution |173j and the p-shell quasihole wave functions |18H I182j but not the complete 
energy dependence of the hole spectral functions. The latter was studied by Benhar et al. |179j . for 
different nuclei using nuclear matter results and LDA, and in Refs. |183t llSUt I184j for ^^O using the 
Green's function formalism. These calculations suggest that most of this strength is found along a ridge 
in the momentum-energy plane (k-E) which spans several hundreds of MeV/c (and MeV). Moreover, 
although the distribution of this strength accounts only for a small fraction of the total number of 
nucleons, it is responsible for as much as 60% of the total binding energy of the system |15Uj . As a 
consequence, the present understanding of nuclear correlations would gain important information from 
the experimental knowledge of this distribution. Such experiments require a high energy transfer from 
the probe to the nucleus and a careful choice of the kinematics in order to minimize the final state 
interactions |185j . For this reason it has only recently been possible to obtain reliable experimental 
data in this correlated region which are still being analyzed |186t I187j . 

Several studies have been carried out to obtain nuclear spectral functions at small and medium 
missing energies, mostly by computing the many-body Green's function j l45j . |188j - |196j . The quali- 
tative features of the strength distribution at small and medium missing energies can be understood by 
realizing that a considerable mixing occurs between hole states and 2hlp configurations. This mixing 
leads to different strength patterns depending on the location of the orbital under consideration. An 
orbital in the immediate vicinity of the Fermi energy will still yield a large fragment near its original 
position since the 2hlp states are quite far in energy. For the same reason, only small components of 
this orbit appear at 2hlp energies. These features are observed for the 3si/2 proton orbit in ^osp]^ 
shown in Fig. EHl An orbit which is surrounded by many 2hlp states will yield a strongly fragmented 
pattern, the width depending on the strength of the mixing interaction. These features are all observed 
experimentally as shown in Fig. 1221 The actual values of the spectroscopic strength shown in this figure 
have not been corrected for the effect of the Coulomb distortion on the electron. This effect is sizable 
in heavy nuclei and leads to an increase of the spectroscopic factor for the 3si/2 strength to 0.65 from 
the values of 0.5 shown in this figure |13j . 

Calculations for medium-heavy nuclei like ^^Ca and ^°Zr gave a fairly good description of the frag- 
mentation pattern of valence hole states |145j in terms of the coupling between quasiparticles and 
TDA/RPA collective modes. The spectroscopic factors obtained in that work typically overestimates 
the experimental summed strength by about a 15% if only these long-range correlations are included. 
The remaining surplus of strength is presumed to be moved to very high missing energies by the effects 
of short-range and tensor correlations as discussed in Sec. |3] This is consistent with the reduction of 
the occupation of low momentum states in nuclear matter also discussed in that section. This result 
from Ref. is supported by the observation that all the removed sp strength for such small momenta 
in nuclear matter ends up at very high energies in the particle domain associated with the repulsive 
core and tensor components of the nucleon-nucleon interaction jnOI (see also Sec. 14. 2 j) . Applying a 
corresponding adhoc reduction of the theoretical strength in medium-heavy nuclei like ^^Ca yields an 
overall satisfactory description of the sp strength distribution as discussed in Sec. 15.11 

The conclusion that all is well is, however, premature. The effects of SRC on the spectroscopic 
factors of ^^O have been computed in Refs. |181j - |183j . All these works yield spectroscopic factors of 
90% for the knockout of a proton for the pi/2 and P3/2 shells in the (e,e'p) reaction before the effects 
of the spurious motion of the center-of-mass are taken into account. Such center-of-mass corrections 
are known to raise the theoretical spectroscopic factor by about 7% |^ I197j . These results are in 
disagreement with the experimental value of about 60% for these orbitals |198j . The inclusion of the 
coupling to 2plh and 2hlp states at low energy has been considered in Refs. jl94 U195] without leading 
to an explanation of this discrepancy. In particular, in Ref. |195j a successful attempt to combine the 
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Figure 36: Spectroscopic factors from the (e,e'p) reaction on ^osp]^ different valence hole 
orbits. The spectroscopic factor for the main fragment of the 3si/2 orbit must be changed 
to 0.65 from the values in the figure to reflect the analysis given in Ref. [TB] . 



65 



treatment of low-energy long-range and short-range (and tensor) correlations was achieved. The final 
result yielded a P3/2 spectroscopic factor of 76%, while still neglecting center-of-mass corrections. This 
leaves a discrepancy with the experiment of about 15-20% (taking into account the quoted uncertainties). 
The calculations of Refs. |194[ I195j take into account the interaction between the hole and particles 
propagating in the system at the TDA level and therefore cannot go beyond the 2hlp level. The 
calculation of Ref. |145j can account for the coupling of quasiparticles to RPA collective modes and 
therefore goes beyond 2hlp contributions, although it is limited by the difficulties in including the 
coupling to collective excitations in both the pp(hh) and ph channels. The effects of fragmentation 
(or self-consistency) were only taken into account in Ref. jl89j were a fully self-consistent calculation 
was obtained with a satisfactory representation of the continuum. However, this calculation does not 
contain the coupling to collective excitations and uses a phenomenological effective interaction. A 
procedure to include both ph and pp(hh) excitations in the construction of the self-energy was recently 
proposed and implemented in Refs. ^31^ 196._ . It was pointed out in this work that the RPA description 
of the excitation spectrum for ^^O is quite inadequate at present and may be partly responsible for the 
unsatisfactory situation for ^^O. 

We begin this section by reviewing some early calculations for finite nuclei in Sec. 15.11 These results 
show that the properties of the spectral fragmentation can be understood in terms of the coupling of 
single particle motion to collective pp(hh) and ph motion, except for the reduction of the spectroscopic 
factors that is due to SRC. A method that allows the evaluation of the effects of SRC on the spectral 
function is reviewed in Sec. 15.21 Calculations for the particular case of ^^O confirm the expected 
reduction of spectroscopic factors due to short-range effects. This calculation also gives quantitative 
information on the distribution of the correlated spectral strength. The location of this strength and the 
effects on the binding energy of nuclear systems are discussed on section The formalism necessary 
to include the effects of both pp(hh) and ph phonons in the nuclear self-energy, while describing them at 
least at the RPA level, is discussed in Sec. 15.41 We present recent results for the spectral function of ^^O 
at low missing energy in Sec. 15.51 and compare them to the available experimental data. The relation 
between the details of the spectral distribution and the results for the excitation spectrum allows one to 
identify relevant missing ingredients necessary to understand the spectroscopic factors of this nucleus. 
A first attempt to improve the description of the spectrum |199j is discussed in Sec. 15.61 

The two-hole spectral function can be computed, including the fragmentation effects, by solving the 
pp(hh) DRPA equation. The high momentum components -excluded form the model space- can be 
included afterwards by computing the relevant defect functions. Some results for two-nucleon emission, 
that clearly show the signature of SRC are discussed in Sec. 15.71 The issue of how to deal with final 
state interactions in reactions at high momentum transfer with regard to the extraction of spectroscopic 
factors is discussed in Sec. 15.81 We conclude this section by discussing two recent experiments that give 
the first direct experimental information of the occupation numbers of deep-lying orbitals in ^osp]^ g^j^^ 
on the spectral strength distribution at high momenta. Sec. 15.91 

5.1 Early results for the spectral distribution 

The earliest calculations of the spectral strength in nuclei were reported in Ref. |188| 1189j . The study of 
the spectral strength distribution starting from a realistic NN interaction is only possible when short- 
range correlations are properly included. These early calculations either relied on the construction of a 
local G matrix interaction obtained from nuclear matter [200J to represent SRC |188j or relied on phe- 
nomenological forces |189j . The latter approach actually generated the first self-consistent calculation 
of the spectral strength in nuclei based on the inclusion of the second-order self-energy diagram. The 
former calculations also used a second-order approach but included the G matrix effective interaction 
without achieving self-consistent results. Useful insights were obtained from these calculations which 
described certain features of the spectral distribution obtained from the ''^Ca(e,e'p) reaction j2Ulj . The 
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Figure 37: Second-order self-energy terms represented by Goldstone diagrams a and b. Ex- 
tensions are made by including TDA or RPA correlations in the particle-particle (diagrams 
c and d) or particle-hole channel (diagrams e and /). 



inclusion of SRC on the distribution of the strength was not considered in this work (see Sec. 15. 2|) . 
As a general conclusion of this work on ^^Ca, it is found that the total strength in the experimentally 
accessible domain in the (e,e'p) reaction is overestimated by about 10-15% although the shape of the 
strength distribution is already well described as shown in Fig. EH These results for the sp strength 
distribution for proton removal from ^^Ca in £ = 2 states are compared to the corresponding experi- 
mental data |2Ulj . The theoretical results are labeled by "2nd Order" in the figure. The normalization 
in the figure is such that if a single state carried all the strength, the peak height at the corresponding 
energy would be 2j + 1. The first peak corresponds to 0^3/2 removal whereas most of the rest of the 
strength corresponds to ^5/2 removal. The difference between the fragmentation of these two sp orbitals 
characterizes the basic features which are found in experiment. For an orbital which in mean field is 
very near the Fermi energy, like the 0^3/2 orbital, fragmentation of strength to all 2plh and lp2h-like 
states takes place in such a way that the distribution is characterized by a single large fragment and 
hundreds of tiny contributions spread out over the whole energy domain covered by the 2plh and lp2h 
states. This background contribution to the hole strength, i.e. the strength which is not contained in 
the area of the main peak, is of the order of 10%. 

The comparison between the second-order results and experiment in Fig. |2H1 also shows that stronger 
fragmentation in the theoretical calculation at low energy is required. Improvement of the description 
of the intermediate states in the self-energy was investigated in Ref . |145j . In this paper various possible 
descriptions of the intermediate 2plh/2hlp in the self-energy were considered. This approach can be 
summarized by considering Fig. |S7| which exhibits the diagrammatic content of these approximations. 
These different approximations correspond to the standard versions of the TDA and RPA for both 
ph and pp(hh) forms as discussed in Sec. El and are represented by diagrams c through / in Fig. ITTI 
These approximations include the effect of long-range (low-enery) correlations on the fragmentation 
pattern at low energy. It should be noted that three independent types of possible collective excitations 
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Figure 38: Distribution of sp strength for 1 = 1 proton removal from ^^Ca. Experimental 
data are from Ref. |2Ulj . Theoretical results are displayed for various approximations to the 
self-energy which are discussed in the text. 



occur, i.e. neutron particle - proton hole, proton particle - neutron hole, and a mixture of neutron 
particle - neutron hole and proton particle - proton hole with isospin components T2=l, -1, and 0, 
respectively. The T2=0 collective low-lying 2+ and 3^ phonons are especially relevant as their coupling 
to the sp motion is an important source of fragmentation of spectral functions at low energies. While 
collective effects are stronger in the RPA, it is possible to encounter serious problems with this method. 
They arise, especially in phRPA, when the RPA correlations are so strong that the lowest solution for 
some angular momentum and parity, has imaginary eigenvalues and the amplitudes cannot be properly 
normalized. This happens for the 3~ phonon, which is crucial for the fragmentation of the | strength 
in "^^Ca. As the imaginary solution can not be included it has to be discarded and the phRPA method 
is unsatisfactory in this case |145j . For the same reason the lowest 1^ states in ^^K en ^°Sc become 
unstable in ppRPA and therefore were also discarded in that approach. These problems are expected to 
disappear when a self-consistent approach would be adopted, which has not been implemented for this 
nucleus sofar. The comparison of the various approximations in Fig. EH to the experimental data shows 
that increased fragmentation at low energy is obtained when collective effects are considered. The sp 
strength farther away from the Fermi energy is largely unaffected. 

The summed sp strength below the Fermi energy generates the occupation numbers for the orbits 
considered as given by Eq. ©. Experimentally, the hole strength has been determined only within a 
limited energy region. For the valence i = and i = 2 shells the main portion of the total strength 
falls within this range of energies. In all theoretical approaches discussed here the calculated extra hole 
strength at higher energies is only 5-10% of the 2j + 1 sum rule for the 2s shell and only 10-15% for 
the Id shell. For the major shells above the Fermi level the occupation probability is about 5-10% and 
for more remote shells only 1-2%. A meaningful quantity is the jump in occupation numbers at the 
Fermi level. This jump in occupation number for the ppTDA result is 0.74. In Fig. EHlthese occupation 
numbers in the ppTDA approximation are multiplied by a factor of 0.9 to simulate the effect of SRC. 
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Figure 39: Proton shell occupation probabilities deduced from a comparison of the present 
calculation and the (e,e'p) data |2Ulj . The figure displays the calculated values with the 
ppTDA correlations in the self-energy, multiplied by a factor 0.9 to simulate the effect of 
short-range correlations. 



The resulting occupation numbers show a flat behavior for deeply bound shells while they bend down 
when the orbits approach the Fermi energy. This feature also emphasizes that LRC are most important 
in the immediate vicinity of the Fermi energy. 

5.2 Depletion due to short-range and tensor correlations 

The treatment of short-range correlations requires the consideration of the details of the two-particle 
wave function at small internucleon distances (or high relative momenta). This suggests that it is 
appropriate to work directly in coordinate or momentum space. The approach introduced in Refs. |2()2| 
I183[I15U] computes the self-energy for finite nuclei in terms of a G-matrix which is obtained as a solution 
of the Bethe-Goldstone equation for nuclear matter 

{M\GsjsKLT\k'e') = {ki\v\k'e') 

+ [ dk\k"f {ki\ V \k"i") ^^\f ^ {k"f\ GsjsKLT we) , (133) 

where k, k', and k" denote the relative momenta between the two nucleons, i, i', and i" the orbital 
angular momenta for the relative motion, K and L the corresponding quantum numbers for the center- 
of-mass motion, 5* and T the total spin and isospin, and Js is obtained by coupling the orbital angular 
momentum of the relative motion to the spin S. We note that Eq. ()133|1 is analogous to the ladder 
equation for the vertex function F, Eq. ()59aj) . except that it describes the intermediate propagation 
of two free nucleons (instead of dressed propagators), taking into account only the effects of Pauli 
blocking with regard to the free Fermi gas. This equation differs from the usual T-matrix of the 
two-body scattering problem only because of the inclusion of the Pauli operator Q that prevents the 
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Figure 40: Graphical representation of the Hartree-Fock (a), the 2-particle 1-hole (2plh, b) 
and the 2-hole 1-particle contribution (2hlp, c) to the self-energy of the nucleon. 



two nucleons to scatter into states inside the Fermi sea (which are already occupied). Equation (|133|) 
therefore generates an appropriate solution of two-body short-range dynamics. The choices for the 
density of nuclear matter and the starting energy ujnm are rather arbitrary. The calculation of the 
corresponding Hartree-Fock term is not very sensitive to this choice [203]. Furthermore, this nuclear- 
matter approximation is corrected by calculating the 2plh term displayed in Fig. EUT b) directly for the 
finite system. This second-order correction, which assumes harmonic oscillator states for the occupied 
(hole) states and plane waves for the intermediate unbound particle states, incorporates the correct 
energy and density dependence characteristic of a finite nucleus G-matrix. To evaluate the diagrams 
in Fig. EHl one requires matrix elements in a mixed representation of one particle in a bound harmonic 
oscillator while the other is in a plane wave state. Using vector bracket transformation coefficients |204j . 
one can transform matrix elements from the representation in coordinates of relative and center-of-mass 
momenta to the coordinates of sp momenta in the laboratory frame in which the two particle state is 
described by 

\hid1k2i2J2JT) , (134) 

where fcj, ii and ji refer to the momentum and angular momenta of particle i whereas J and T define 
the total angular momentum and isospin of the two-particle state. Performing an integration over one 
of the ki, one obtains a two-particle state in a mixed representation 

POO 

\n1i1j1k2i2J2JT) = dk,klRn„i,ik,) \k1i1j1k2i2J2JT) . (135) 
Jo 

Here Rm/i stands for the radial oscillator function and the oscillator length is chosen to be 6 = 1.72 fm^^ 
to achieve an appropriate description of the bound states of ^^O. Using the results of Eqs. (jl33p - (|135|) 
the Hartree-Fock approximation for the self-energy can be obtained in the momentum representation 

^iSih, k[) = ^ / Yl + + 1) iki^iJi^2i2j2\ GjT \k'Ann2i2j2) . (136) 

The summation over the oscillator quantum numbers is restricted to the states occupied in the IPM of 
^^O. This Hartree-Fock part of the self-energy is real and does not depend on the energy. 

The terms of lowest order in G which give rise to an imaginary part in the self-energy are represented 
by the diagrams displayed in Figs. EHT b) and I40r c). refering to intermediate 2plh and 2hlp states, 
respectively. The 2plh contribution to the imaginary part is given by 



Wlf{k^,k[;E)= ^ E EE / ^'^^ / K'dK{2J+l){2T + l) (137) 

^ -^^ ^ n2£2i2 JST ^ 

X (A;i£ijin2£2j2| GjT \kiSKL) {kiSKL\ Gjt \k[iijin2i2j2) tt 5 ( E + 6^,1,^, - ^ - , 

\ 4m m / 
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where the average experimental quasihole energies en2i232 ^'^^ used for the hole states (—47 MeV, —21.8 
MeV, —15.7 MeV for si/2, P3/2 and pi/2 states, respectively), while the energies of the particle states 
are given in terms of the kinetic energy only. The plane waves associated with the particle states in 
the intermediate states are properly orthogonalized to the bound sp states following the techniques 
discussed in Ref. [202 . The 2hlp contribution to the imaginary part W'l^^^{pi,p'i] E) can be calculated 
in a similar way |202j. 

The choice of pure kinetic energies for the particle states in calculating the imaginary parts of W"^^^^ 
(Eq. p37|) ) and W"^^^'^ may not be very realistic for the excitation modes at low energy. Indeed a sizeable 
imaginary part in W^^^^ is obtained only for energies E below -40 MeV. However, when the primary 
interest is in the effects of short-range correlations, the choice appears appropriate since these involve 
excitations of particle states with high momenta. A different approach is required to treat the coupling 
to the very low-lying 2plh and 2hlp states in an adequate way, as discussed in Sees. 15. II and 15.41 The 
2plh contribution to the real part of the self-energy can be calculated from the imaginary part W"^'^^^ 
using a dispersion relation |205j 

V r Wf''.^(h,k\;E') 
v.f^^'ih, k[-E) = L-j ^ E' Ie ^^^^^ 

where V stands for a principal value integral. A similar dispersion relation holds for K^^^p and VT^'^^p. 

Since the Hartree-Fock contribution S^^ has been calculated in terms of a nuclear matter G-matrix, 
it already contains 2plh terms of the kind displayed in Fig. EUfb). In order to avoid this overcounting 
of the pp ladder terms, one subtracts from the real part of the self-energy a correction term ( V^) , which 
just contains the 2plh contribution calculated in nuclear matter. Summing up the various contributions 
one obtains the following expression for the self-energy 

S = S"^ + AS = S"^ + (\/2plh _y^j^ y2hlp^) ^ (^^2plh ^ ^2hlp^ _ ^^gg) 

The Dyson equation ()13p for this self-energy can be solved in momentum space. In Ref. |150j the 
integrals were discretized by considering a complete basis within a spherical box of a radius -Rbox- The 
calculated observables are independent of the choice of -Rbox, if it is chosen to be around 15 fm or 
larger. By taking into account the different normalization of these basis functions from the case of 
plane waves in the continuum, one can express the matrix elements of nucleon self-energy in the "box" 
basis. Energies and wave functions of the quasihole states can be determined by diagonalizing the 
Hartree-Fock Hamiltonian plus AS in this "box basis" [see Eq. fl37p ]: 



s 

n=l 



{h\ A.<5,„ + Sf + AS,,(i? = at,,) \kn) (fc„|^"->,^. = et,, (^.k""),, • (140) 



In this approach AS only contains a sizable imaginary part for energies E below the quasihole eigen- 
values s'^n^f^j {n— identifying discrete solutions). Therefore the solutions of Eq. (|140|) are separated in 
energy from the continuum contribution to the spectral function. The eigenvector corresponding to 
these discrete states yields the quasihole wave function {k\z'^~)^ - in momentum space, which still needs 

to be normalized by the spectroscopic factor Z'^_f^- by means of Eq. (jHEI)- The diagonal part of the 
quasihole contribution to the spectral function is given, in the box basis, by 

^f-.,(fcn; E) = Zl% |(fcn|^"->,/ 5{E - et,,) . (141) 

In the calculations described in Refs. |150| I184j . the Bethe-Goldstone equation ()133|) was solved 
by employing for V the one-boson-exchange potential Bonn-5 developed by Machleidt in Ref. 
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Figure 41: Square of the quasihole wave function for the pi/2 state in ^^O (full curve), 
normalized to the spectroscopic factor according to Eq. (|38|) . compared to the Hartree-Fock 
result (dashed curve). 



(Tab. A. 2) and the Pauli operator Q was approximated by the so-called angle-averaged approximation 
for nuclear matter with a Fermi momentum kp = 1.4 fm~ . This roughly corresponds to the saturation 
density of nuclear matter. The starting energy unm for computing the G-matrix, Eq. ()133|) . was chosen 
to be -10 MeV. 

The square of the quasihole wave function for the pi/2 state (normalized to the spectroscopic factor, 
"^opi/2) shown in Fig. |^ as a full line. For comparison the result for the Hartree-Fock wave function 
is shown by the dashed line. From the comparison one can infer that at the quasihole energies no 
substantial change in the wave functions occurs and that the Hartree-Fock wave function is a good 
approximation. It should be further noted that the wave function of a Woods-Saxon potential, which 
is constructed as the local equivalent of the Hartree-Fock potential |2U2j . is indistinguishable from the 
Hartree-Fock wave function. This suggests that the explicit inclusion of short-range correlations does not 
lead to the strong suppression of the wave function in the interior of the nucleus as has been suggested 
by Refs. |2U7l I38j . These findings are in line with the observation that the quasiholes representing mean- 
field orbitals can be adequately described by a local potential |25- It should be emphasized that this 
is not necessarily the case for other overlap functions since the coupling to collective excitations 
can lead to significant changes to the quasihole wave function |^ l2U8j . These effects are unlikely to 
generate changes in the high-momentum components discussed here. It is important to note that these 
results also show that the influence of high-momentum components in the quasihole wave function is of 
minor importance, so the enhancement of the momentum distribution has to come from excitations at 
higher missing energies. 

We show in Fig.|321an example of the reduced cross section for the (e,e'p) reaction on ^^O leading to 
bound quasihole state of ^^N at an excitation energy of —6.32 MeV |2U9j . In this picture the effects of 
FSI have been computed in the distorted wave impulse approximation (DWIA) [16j and the data points 
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Figure 42: Reduced cross section for the ^^0{e,e'p) reaction in parallel kinematics leading 
to the 3/2~ state at —6.32 MeV of the residual nucleus ^^N. Results of the Green's func- 
tion approach (solid line) are compared to those obtained in the variational calculation of 
Ref. |181j (dashed line) and the experimental data |198j . A spectroscopic factor of 0.537 
was required for the Green's function result, while .^op3/2 = 0.459 has been used to for the 
variational calculation. 



have been obtained at NIKHEF for the so-called parallel kinematics |198j . Using the quasihole part of 
the spectral function computed from Eq. p40|) but adjusting the spectroscopic factor for the quasihole 
state contribution ^op3/2 experimental data, one obtains the solid line of Fig. |121 Comparing 

this result with the experimental data one finds that the calculated spectral function reproduces the 
shape of the reduced cross section as a function of the missing momentum quite well. The absolute value 
for the reduced cross section can only be reproduced by employing a spectroscopic factor ^op3/2 ~ 0.537, 
a value considerably below the one of 0.914 calculated using Eq. p40|) |150j . An analogous result holds 
for transition to the pi/2 ground state of ^^N for which the same analysis yields a spectroscopic factor 
^opi/2 ~ 0.644. For comparison, one may note that the use of phenomenological Woods-Saxon wave 
functions adjusted to fit the shape of the reduced cross section require spectroscopic factors ranging 
from 0.61 to 0.64 for the lowest Opi/2 state and from 0.50 to 0.59 for the 0p3/2 state, respectively, 
depending upon the choice of the optical potential for the outgoing proton |198j . 

Figure |321 also contains the results for the reduced cross section derived from the variational result 
for overlap wave function of Pieper et al. |181j . who employed the Argonne fi4 potential for the NN 
interaction |135j . Also in this case the shape of the experimental data is globally reproduced with a 
slightly better agreement for small negative values of pm but with a clear underestimation at larger pm- 
The overall quality of the fit is somewhat worse than for the Green's function approach and the required 
adjusted spectroscopic factor is Z^^^^^ = 0.459, below the value of 0.537 for the latter approach. It is 
not clear, however, whether the differences in the calculated reduced cross section are due to the use of 
different interactions or to the methods employed in calculating the spectral function. 

Both the variational calculations of Refs. |18H I182j and the Green's function approach give spec- 
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trocsopic factors for the p-shell orbitals of about 0.90. This value should be compared to the ~0.63 
obtained form the NIKHEF experimenf for these shells. The additional contribution due to the proper 
consideration of the center-of-mass motion raises the theoretical value to 0.98, worsening the agreement 
with data |197j . Nevertheless, the observation that the same results are obtained with two independent 
many-body methods but including the same short-range physics suggests that the effects of SRC on the 
quasihole strength are well under control. At the same time, the discrepancy with the experiment is 
partly due to the emphasis on the accurate treatment of short-range correlations only and one should 
view the quasihole strength that has been discussed here to be due only to the influence of short-range 
correlations |202j . It is clear that a considerable renormalization of the strength is to be expected 
due the coupling of the quasihole states to the low-lying collective excitations, as will be discussed in 
Sec. EH 



5.3 High-momentum components at high missing energies 

The continuum part of the hole spectral strength is found at higher energies and stems from the coupling 
to the continuum of 2hlp states. In this region it is useful to proceed from the Hartree-Fock propagator 
with states |a) that diagonalize the corresponding self-energy 

9S'(^-^E)= I (142) 

where the sign in front of the infinitesimal imaginary quantity if] is positive (negative) depending on 
whether e^^^ is above or below the Fermi energy. The Dyson equation can then be solved by iterating 
the AS component in Eq. (jl39|) of the self-energy to generate the reducible self-energy 

{a\ E^fiE) 1/5) = {a\ AS,,(E) + ^ {a\ AS,,(E) I7) ^^g^T; E) (7I ^fiE) \P) (143) 

7 

and obtain the propagator from 

gej{a, (3; E) = 6^,^ gf{a- E) + gf^{a- E) {a\ i:^,f{E) |/5) gf{^- E). (144) 

Using this representation of the Green's function one can calculate the spectral function in the "box 
basis" from [see Eq.(jl))] 

S'',^{kn,E) = - Im \Y,{kn\a),^ giAa,f3;E) {f3\K)J . (145) 

This spectral function is different from zero for energies E below the lowest sp energy of a given Hartree- 
Fock state (with ij) only due to the imaginary part in E'"^''. This contribution involves the coupling to 
the continuum of 2hlp states and is therefore nonvanishing only for energies at which the corresponding 
irreducible self-energy AS has a non-zero imaginary part. The 2plh contribution to the self-energy is 
responsible for the depletion of strength, which in mean field is located below the Fermi energy, to 
high energy. The 2hlp term instead, is essential for the accumulation of sp strength below the Fermi 
energy from states (in particular those with high momenta) which are empty in the mean field. The 
continuum contribution of Eq. ()145|) and the quasihole parts of Eq. ()14H) . which are obtained in the 
basis of box states, can be added and renormalized to obtain the complete spectral function in the 
continuum representation at the momenta defined by the box basis from 

S,Akf,E) = ^ l^^{S^^^{kf,E) + J2ste,{h-,E)), (146) 
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Spectral Function for Pj^yg 
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Figure 43: The pi/2 spectral strength as a function of momentum at fixed energies corre- 
sponding to -50, -150, and -250 MeV. The results demonstrate the increasing importance of 
high-momentum components with higher excitation energy in the A — 1 system. 



where the Na correspond to appropriate normalization constants |150j . 

The results for the spectral function S^j{k;u) for the pi/2 quantum numbers are shown at three 
different energies in Fig. H^l The long-dashed curve corresponds to -50 MeV, the full curve to -150 
MeV, and finally the short-dashed curve to -250 MeV. From these results it is clear that an important 
change in the momentum content of the sp strength occurs with increasing excitation energy in the 
A = 15 system. At higher excitation energy one finds more high-momentum components. Moreover, 
these high-momentum components are not observed in the quasihole states. This can be concluded 
by considering Fig. where the total momentum distribution, including the contribution from the 
quasihole states, is shown in the left panel. This requires the energy integration of the continuum hole 
strength for each k according to Eq. 0. This distribution is presented for various energy cut-offs. The 
quasihole part reflects the cross section for knockout reactions with small energy transfer, i.e. leading 
to the ground state of the final nucleus and excited states up to ~ 20 MeV. The curve denoted by > 
-100 MeV reflects the momentum distribution including all states of the final nucleus up to around 
80 MeV, etc. This result for ^^O is very similar to the observation for ^He made in Ref. |175j . where 
the contribution of the ground state to ground state transition exhibits also very few high-momentum 
components. As a consequence the high-momentum components of the momentum distribution due to 
short-range correlations can be observed mainly in knockout experiments with an energy transfer of 
the order of 100 MeV or more. The right panel of Fig. compares the total momentum distribution 
for ^^O obtained with different approaches, and different interactions, all of which properly include the 
effects of short-range correlations |184j . 

To understand this result, it is important to recall that the appearance of high- momentum compo- 
nents at a certain energy in the A — 1 system is related to the self-energy contribution containing 2hlp 
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Figure 44: The total momentum distribution of ^^O. The left panel also displays the quasi- 
hole contribution and the results obtained with various energy cut-offs in the integration 
of the spectral functions. The right panel compares the momentum distribution obtained 
with the Green's function (total) with that of other approaches: local density approxi- 
mation (LDA) |179j . Fermi hypernetted chain (FHNC) [210J and Variational Monte Carlo 
(VMC) |173j . It should be noted that different normalizations were used for the presentation 
of these panels. 



states at this energy. From energy conservation it is then clear that at low energy it is much harder to 
find such states with a high-momentum particle than at high energy. This same feature is observed in 
nuclear matter where the peak of the sp spectral function for momenta above kp increases in energy as 
k'^ |211j . As a result, the hole strength in nuclear matter as a function of momentum shows the same 
tendency as the result shown in Fig. 021 [HHIi ^-e. higher momenta become more dominant at higher 
excitation energy. 

In order to show the importance of the continuum part of the spectral functions as compared to 
the quasihole contribution and to visualize the effects of correlations, we have included in Tab. Q the 
particle numbers for each partial wave including the degeneracy of the states 

/£p poo 
dE j dkk'^Sij{k,E) , (147) 
oo Jo 

also separating the contributions originating from the quasihole states and those due to the contin- 
uum [as in Eq. ()146|) ]. As already noticed above only 14.025 out of the 16 nucleons of ^^O occupy the 
quasihole states in this calculation (while the experimental data suggests a much smaller number). An- 
other 1.13 nucleons are found in the 2hlp continuum with partial wave quantum numbers of the s and 
p shell, while an additional 0.687 nucleons are obtained from the continuum with orbital quantum num- 
bers of the d and / shells. The distinction between quasihole and continuum contributions is somewhat 
artificial for the Si/2 orbital since the coupling to low-lying 2hlp states leads to a strong fragmentation 
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E 


14.025 
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0.824 
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15.841 





Table 1: Distribution of nucleons in ^^O. Listed are the total occupation number fi for various partial 
waves [see Eq. ()147p ] but also the contributions from the quasihole (n'^'*) and the continuum part {n'^) 
of the spectral function, separately. The continuum part is split further into contributions originating 
from energies E below -150 MeV {n'^{E < —150)) and from energies below -100 MeV. The last line 
shows the sum of particle numbers for all partial waves listed. 



of the strength |212j . which is also observed experimentally |213j . It should be noted that the depletion 
of the occupation probabilities of the hole states, indicated in Tab^ is larger for the si/2 orbit. This 
feature can be ascribed to the closeness of the Si/2 Hartree-Fock energy to the 2hlp continuum which 
yields more leakage of strength to the continuum than for the pi/2 and P3/2 quasihole states. The sum 
of the particle numbers listed in Tab. Q is slightly smaller (15.841) than the number particle in ^^O 
this due to the fact that only partial waves up to £ = 3 were taken into account j214|. One must also 
keep in mind that the approach to the sp Green's function discussed here is not number-conserving, as 
the Green's functions used to evaluate the self-energy are not determined in a self-consistent way, as 
dicussed in Sec. 12.31 

The contributions to the total energy, as derived from the Migdal-Galitski-Koltun sume rule ()12|) . 
are shown in Tab. |21 for different angular momenta. The first two columns give the analogous results 
obtained from the solution of the Hartree-Fock (HF) and Brueckner-Hartree-Fock (BHF) terms. The 
latter result includes the 2plh correction to the nuclear-matter G-matrix (see Ref. |l15.(I] ) . The BHF 
result continues to describe the nucleus in terms of fully occupied sp states as in HF. However, as the 
sp states in BHF are more bound, the gain in binding energy from HF to BHF is accompanied by a 
reduction of the calculated radius of the nucleon distribution. The inclusion of the 2hlp contributions 
to the self-energy in the complete calculation reduces the absolute values of the quasihole energies 
(compare BHF and "Total" in Tab. |2I). Despite this reduction of the quasihole energies, however, the 
total binding energy is increased as compared to BHF. This increase of the binding energy is mainly 
due to the continuum part of the spectral function. Comparing various contributions to the intergral in 
Eq. ()12|) . one finds that only 37% of the total energy is due to the quasiholes Eq. p4H) . The dominating 
part (63%) results from the continuum part of the spectral functions although this continuum part only 
represents 1.8 nucleons (that is 11% of the total, see Tab.[T)). 

Summarizing, the calculation of the complete energy dependence of the hole spectral function demon- 
strates that the presence of high-momentum components in the nuclear ground state will only show 
up unambiguously at high excitation energy when probed by (e, e'p) reactions. These deeply bound 
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Ij 


e 


HF 

t 


AE 


e 


BHF 

t 


AE 


e 


Total 
t 


AE 


si/2 qh 


-36.91 


11.77 


-50.28 


-42.56 


11.91 


-61.30 


-34.30 


11.23 


-35.98 
















-90.36 


17.09 


-22.89 


P3/2 qh 


-15.35 


17.62 


9.08 


-20.34 


18.95 


-5.59 


-17.90 


18.06 


0.37 


P3/2 C 














-95.19 


35.19 


-9.96 


Pi/2 qh 


-11.46 


16.63 


10.34 


-17.07 


18.46 


2.76 


-14.14 


17.19 


5.47 


Pl/2 C 














-103.62 


35.94 


-5.84 


/ > 1 c 














-98.87 


63.17 


-12.27 


E/A(MeV) 
(r)(fm) 


-1.93 
2.59 


-4.01 
2.49 


-5.12 
2.55 



Table 2: Ground-state properties of ^^O. Listed are the energies e and kinetic energies t of the quasihole 
states (qh) and the corresponding mean values for the continuum contribution (c), normalized to 1, for 
the various partial waves. Multiplying the sum: ^(t + e) of these mean values with the corresponding 
particle numbers of Tab. I, one obtains the contribution AE to the energy of the ground state [as given 
by the Migdal-Galitski-Koltun sumrule, Eq. fT^ ]. Summing up all these contributions and dividing 
by the nucleon number yields the energy per nucleon E/A. Furthermore, the rms radius for nucleon 
distribution is shown. Results are presented for the Hartree-Fock (HF), Brueckner-Hartree-Fock (BHF) 
and the complete calculation of Sees. 15.21 and 15.31 (Total). The particle numbers for the qh states in 
HF and BHF are equal to the degeneracy of the states, all other occupation numbers are zero. All the 
energies are given in MeV. 

nucleons, not only generate the enhancement of the momentum distribution for momenta > 400Mev/c, 
depicted in Fig. but they are essential in understanding the binding of nuclear systems. More 
discussion of these high-momentum components is given in Sec. 15.91 

5.4 Faddeev approach to the treatment of collective excitations 

The results for the spectroscopic factors of Sec. 15. ll demonstrate that the understanding of the occupation 
of quasihole states, can only be achieved with a complete description of the coupling to the low-energy 
collective motion of the system. In this section, we review a method that has been recently developed 
in Refs. j^H 1196] , with the aim of pursuing a complete description of these features. 

The degrees of freedom that contribute to long-range correlations can be very complicated. However, 
when one is interested in the dynamics at small missing energies, close to the Fermi level, the problem 
can be simplified by focusing only on the excitations at low energy, generally a small number. If the 
other phonons away from the Fermi surface are high in energy they can be expected to mix weakly 
or not at all with the low-energy configurations. It is important to observe that this picture can be 
too optimistic if one works in terms of particle and hole orbitals as generated by an IPM. For example 
the direct diagonalization of the nuclear hamiltonian in a shell-model basis requires the inclusion of 
configurations at rather high energy (several hu), while still obtaining sizable mixing of strength in 
the lowest excited states (see for example Ref. j215|| ). In this respect, the Green's function approach 
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appears to account more easily for the shift in energy of the nuclear strength, as discussed in relation 
to the results of Fig. IHHl The coupling to 2plh fragments is responsible for moving sizable amounts 
of strength into the particle domain, even at of several tens of MeVs from its original location. This 
mixing comes into the sp propagator through the solution of the Dyson equation and therefore it is 
included in the description of the low-energy excitations. One thus expects that the dressed propagator 
represents a more accurate approximation of the sp excitations. As a result, the decoupling between 
low-energy excitations and those at higher energy is more effectice when the full fragmentation of the 
spectral strength is taken into account, as in the SCGF theory. To some extent, the use of fragmented 
propagators corresponds to expanding the wave function in terms of an already correlated single particle 
basis, as in the Monte Carlo shell model of Ref . |216j . From a more practical point of view, the energy 
distribution of the spectral function and the quenching of the spectroscopic factors act as to renormalize 
the interaction between quasiparticles and phonons (the real degrees of freedom in the system) allowing 
for a better convergence of the results. A second motivation to follow the approach of SCGF comes from 
the observation that a fully self-consistent calculation guarantees the fulfillment of the basic conservation 
laws and sum rules (see Sec. 12. 3|) [23 tHHI • 

In general, one may need to allow the sp motion to couple to both pp(hh) and ph motion and at the 
same time achieve a satisfactory description of these collective phonons. It has been known for a long 
time that the excitation spectrum of ^^O, as obtained in TDA approximation, is particularly inadequate. 
In order to account for the coupling to collective excitations that are actually observed in this nucleus 
it is necessary to at least consider an RPA description of the isoscalar negative parity states |161j . 
To account for the low-lying isoscalar positive parity states an even more complicated treatment will 
be required. Sizable collective effects are also present in the pp and hh excitations involving tensor 
correlations for isoscalar and pair correlations for isovector states. Each of the diagrams c) to /) of 
Fig. IH7I represents one of these couplings in a separate way. Moreover one must also remember that the 
diagrams coupling the ph and sp channel are affected by a violation of the Pauli principle already at the 
2plh/2hlp level, due to neglect of the exchange between the freely propagating hole or particle with 
the one propagating inside the ph phonon. More serious inconsistencies show up if one naively sums the 
above diagrams, with the intent of accounting for both pp(hh) and ph collective modes. In this case a 
double counting of the second order diagrams. Fig. HTh ) and b) , would come into play and there is no 
simple way to correct for this except to introduce spurious poles in the self-energy jST]. Therefore the 
inclusion of both the pp (hh) and ph collective phonons does require a complete treatment by summing 
these excitations to all orders. Such an expansion sums a series of diagrams like the one shown in Fig. 1451 
The interplay between the various phonons would ultimately generate an expansion of the 2plh/2hlp 
propagator R^u) [see Sec. 12. 2 j that automaically resolves the above issues. In Ref. jSI], it was shown 
that this all order summation can by achieved by extending the technique of the Faddeev equations 
for the three-body problem |2171 12181 1219j to the case of three interacting quasiparticle excitations. A 
corresponding approximation for the self-energy can then be obtained from Eq. ()19|) . 

To describe the formalism of the Faddeev equations as applied to the many-body case, we consider 
the Bethe-Salpeter equation for the 2plh propagator [which is analogous to Eq. (PTj) for the ph case]. 
In the time formulation, one has 

RfiuX,al3-y{tl,t2, t') 

= g,j.a{ti - t')gup{t2 - t')g^x(t' - ta) - guo,{t2 - t')g^ii3{ti - t')g^x{t' - ts) 
+9^tf.'iti - i!^gvv'{t2 - t2)gx'x{ts - h) 

where here and in the rest of this section we use the convention of summing over all repeated indices and 
integrate from — oo to -|-cx3 over all repeated time variables, unless specified otherwise. The interaction 
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Figure 45: Example of diagrams that are 
summed to all orders by means of the 
Faddeev equations. 




vertex is given by 



K^uX,Ql3"/(tl,t2, ts; ^4, ts, tg) 

+ K%%{t,M-MM)9~l{h-t,) + i^S/37(^l'^2,t3;t4,t5,t6). (149) 

In Eq. p49|) -R"^^^^ represent the ph irreducible vertex that appears in Eq. (EH) while K^pp^ and K'^pp^^ are 
the pp and 2plh irreducible vertices. It should be noted that in Eq. (fT^ the propagator R^yx^apy is only 
required at two times and therefore the complete knowledge of its time dependence is not necessary 
to obtain the self-energy. On the other hand, the dependence on the time variables ti, ^2 and ^3 is 
employed in the Bethe-Salpeter Eq. ()148p . thus an exact solution of the 2plh motion would require at 
least a 4-time object. 

Equation ()148|) can be reduced to a set of coupled equations in a way similar to the method proposed 
by Faddeev to solve the three-body problem |217| I218j . The inclusion of pp and ph RPA phonons in 
a consistent way requires an approach that provides a natural framework for correctly iterating in the 
self-energy quantities that have already been summed to all orders, like pp and ph RPA phonons. The 
calculation reviewed in Sec. I5.5l neglect the contribution of the irreducible /^(pp^) term in Eq. ()149|) and 
therefore require only three Faddeev components. Following standard notation in the literature |219j . 
the component R^^^.^ ajSy related to the sum of all diagrams ending with a vertex between legs j and 
/c, with {i,j,k) cyclic permutations of (1,2,3). We will employ the convention in which the third leg 
propagates in the opposite direction with respect to the first two. The Faddeev components R^'^^ can be 
written in terms of the 2plh propagator R and the contribution of the three dressed but noninteracting 
sp propagators. This definition is given in detail here for all three components, omitting explicit reference 
to the time variables for convenience of notation 



R 



(1) 



9..9px R^ 



2 i.9p.a 9ul3 9-1X — 9ua 9fJ.l3 9^x) 



R 



R 



(2) 

fiuXyO/Bj 
(3) 

'fj.uX,af3'y 



9p.e9pX -^ep no- ^ 



9ul3 9lX - 9ua 9fil3 9ix) 



9p.e9up K^p% Rr,aX,al3'y + ^ {9 pa 9ul3 9lX " 9va 9pl3 9-yx) , 



(150a) 
(150b) 
(150c) 
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where | is a symmetry factor. With these definitions the full propagator R{uj) is given by 



1=1,2,3 

The Faddeev equations now take the following form 



R^il\,ap-i - 7;i9fia gup g-yx- gua QfifB g^x) ■ (151) 



1 

fiiyXjOp-y 2 



^nlx,al3^ — o ^Sfia gu/B g^yX " fi'i/a fi'^/? A'ta) 



+ fl'A'A r^VA',/t'V"A" (-^/^"u" X" ,al3'r + ^ ^l" u" X" ^afi-y) ' — 1, 2, 3 (152) 

where the F^*^^ vertices obey the following symmetry relations and are defined by 

^^1^,0/37(^1' ^2, ^3; ^4, ^5, h) = r[,^A,/3a7(^2, ^1, ^3; h, ^4, ^e) = g ^lih " ^4) f [,^'^^^(^2, ^s; ^5, U) , (153a) 



r^tA,a/37('^i> ^2, ^3; ^4, ^5, ^e) - T^f j;^^^^^ (^2 , i^i, ^3; ^5, U, te) — g^l{h - h) r[fj^i/3(ti, t2] t4, ^s) • (153b) 

In Eqs. p53|) . (^^^(r) is the inverse sp propagator and matrices F^pp) and F^p'^) correspond to the proper 
time orderings of the four-point reducible vertex r^~^* that appears in Eq. (fT^ . These quantities solve 
the Bethe-Salpeter equation for the pp and ph motion, respectively 

^^^11(3(^1^^2] t3,t4) = -ft'^5^2_g(tl,t2;t3,^4) 

+ ^i'Utu h] t[, Q g,,{t[ - 4) g^^it', - Q K^l'lf^itl tl; t,, U) , (154a) 



+ f 52.(^1, t2] t'l, Q g,M - Q gUt', - Q K%%{t',, t',; t3, U) . (154b) 

Apart from neglecting the K^p'^^^ vertex, Eq. (jl52j) is otherwise an exact equation for the 2plh 
propagator, that involves quantities which depend on several times and for this reason still intractable. 
In order to construct a manageable approximation scheme one first requires a reduction to a set of 
equations that involve only two-time quantities but still include the relevant physical ingredients of 
interest. Here we give only a brief overview of this issue. More datails can be found in Refs. PH I220j 

The calculations of Sec. 15.51 employ a bare interaction Va/s^ys for the vertices K^^^^ and K^p'^'^ so 
that the Bethe-Salpeter equations (jl54|) reduce to the usual dressed RPA (DRPA) equations [3^ I221j . 
The solutions of these equations depend only on two times. These pp and ph phonons correspond to 
the dressed version of the phonons that are considered in Ref. jl45j (see also Fig. EZj) and represent 
the minimum step that maintains the simultaneous inclusion of both pp and ph collective low-lying 
excitations in the self-energy. However, it is possible to extend these vertex to include correlations that 
go beyond the RPA level [see Sec. 15. 6j . A second approximation requires to construct the Faddeev 
equations ()152|) for only two time variables. In doing this, the sp line that propagates freely in the 
Faddeev vertices p53|) remains blocked and can propagate only in one time direction. This means 
that the set of Eqs. p52|) split up in two separate expansions for the 2plh and the 2hlp components. 
For this reason, the pp and ph phonons will be summed only in one time direction in a TDA way 
contributing separately to the 2plh and 2hlp components of the self-energy. However, the collective 
RPA correlations in the pp and ph channels have already been computed through Eqs. ()154j) and 
therefore remain properly included in the this approach. An example of the diagrams that are included 
in the pp channel, Eq. ()153b|) . is shown in Fig. UHl 
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Figure 46: Diagrams that are included in the definition of the vertex for the pp channel. 
Here AF^ and AF^ are the forward- and backward-going part of the energy dependent 
contribution to the pp DRPA vertex (jl54ap . The contribution of these three diagrams can 
be factorized in an expression of the form G^'^ F'^'^^ G^^ only after having redefined the 
propagators G^^ and F*^^) to depend also on the particle and hole fragmentation indices 
{n,n',k). 

The remaining complication, related to the use of dressed propagators, concerns the interactions 
vertices (|153|) . In order to keep track of which quasiparticle or quasihole is coupled to the DRPA 
phonons, the F*^*) and the propagators R^^^ need to be redefined in such a way that their matrix elements 
also depend on the indices (n, n', k), which label the fragments of the sp propagator, Eq. 0. This 
implies that the eigenvalue equations will involve summations on both the sp indices {a, (3, 7) and 
the ones corresponding to the fragmentation, (n^, n/3, k^). The particular expression to be employed 
for the vertices F*^*-* depends on the approximation chosen for the pp and ph phonons and it can be 
crucial to control spurious solutions in the Faddeev equations. Here we only note that the diagrams 
of Fig. |3ni are those required for the DRPA, as discussed in Ref. jSl- The 2plh propagator and its 
Faddeev components, as defined in Eqs. ()150|) . are recovered only at the end by summing the solutions 
over all values of (n^, n^, k^) and (n^, n^, kx). 

Putting together all the above considerations, the resulting approximation to the Faddeev equa- 
tions (|152p can be rewritten in a way where all the propagators involved depend only on one energy 
variable (or two time variables). The forward-going part of this expansion can be written as follows 
(after Fourier transformation to the energy representation) 



d(«) ( \ 

1 



2 \^ fin^vnv\kx,cmaPnij'yk-yV^) ^ vnv^n^\kx^anaPnij'^k-^\ i 
"I" ^ tjn,,nn,.\k^ 1 1' n' 1/ n' \' k'. i^) T, ,/„/ ,,/„/ \iu \" 



vn„^Jm^\kx,^J.'n'^v'K\'k'^\'^ ) u'nl^l'n'^\'k'^,^l"n'|[u"n',^X"k^ 



\"k'l,cxnal3np'yk^i^) + ^fi"n'i^u"n'^X"k'l,ancl3npik.y(^)j ' "^ — 1,2,3(155) 

In Eq. p55|) . G^^ is the forward-going part of the 2plh propagator for three dressed but noninteracting 
lines. This quantity is therefore defined as follows (without implicit summation): 

^ im^un„Xkx,anaf3n,S"/k^,i^) ~ ^nf_,,na ^ni,,np ^kx,k^ 

A completely analogous set of equations apply for the backward-going (2hlp) expansion of i?(a;). 
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Equations ()155|) reduce the general "Faddeev-Bethe-Salpeter" expansion, Eq. ()152|) . to a tractable 
set of equations involving only two-time objects. It is important to note that these equations are still 
expressed in terms of the self-consistent solution gapil^) and include both the collective pp and ph 
phonons in a correct way. Thus they maintain all the features relevant for the physics that have been 
discussed at the beginning of this subsection. 



5.5 Single-particle spectral function of ^^O 

A self-consistent solution of the spectral function of ^^O was obtained in Ref. |196j by solving the set of 
Faddeev equations (jlSSj) . The Dyson equation was solved within a model space consisting of harmonic 
oscillator sp states. An oscillator parameter h = 1.76 fm was chosen (corresponding to hu = 13.4 MeV) 
and all the first four major shells (from Is to 2plf) plus the lgg/2 where included. Inside this model 
space, the interaction used was a Brueckner G-matrix |222j derived from the Bonn-C potential j2(J6j . 
For the solution of the Faddeev equations a G-matrix evaluated at a fixed starting energy uj =-25 MeV 
was chosen as a suitable average of the energy the most important 2hlp states that couple to the 
experimentally observed quasiholes [198]. However, the energy dependence of G{uj) was taken into 
account in the calculation of the mean field part of the self-energy, resulting in the Brueckner-Hartree- 
Fock approximation (this can be compared to the HF case in Eq. (jTH)) ), 



Sfr(^) = E / ^G^^M^ + ^')9.s{co') (157) 

The self-energy was then computed as in Eq. (fT^ . with Eq. (jl57p instead of T,^^. 

The normalization (spectroscopic factor) of each sp fragment was computed by means of Eq. (jSHl), 
which we rewrite here in the form 



, „^ y',- (158) 

a,l3 



where ~ (\E'^~^|cQ,|\EfQ ) represent the A;-th solution for the spectrocopic amlitudes and ~ Eq — 
E^~^ the corresponding eigenvalue. It is important to note that the Ti^^^iu) component of the self- 
energy also contributes to reduction the spectroscopic factor. This energy dependence, associated with 
the use of the G-matrix in Eq. p57|) . is not present in a standard HF self-energy and originates from 
the states at high energy that are excluded from the model space. In Ref. |195j . it was found that this 
procedure properly accounts for the depletion of spectroscopic factors that is induced by SRC, at least 
for the normally occupied shells in the IPM. 

The solution of RPA equations may lead to collective instabilities whenever interactions are employed 
that are too attractive. A particular sensitivity to the strength of the G-matrix interaction was found in 
the calculation of |196l I199j , only for the lowest isoscalar 0"*" solution. For this reason an artificial shift 
of its eigenvalue to the experimental energy of first excited state, 6.049 MeV, was maintained during the 
successive iterations. Note that the screening due the redistribution of the sp strength tends to reduce 
the instability problem in DRPA. Moreover, the origin of the instability for this particular solution has 
been identified and corrected later, in Ref. |199j . 

The Faddeev equations were solved first in terms of a IPM input propagator in both TDA and RPA 
approximations. For an IPM ansatz, the TDA calculation is equivalent to the one of Ref. |195j and yields 
the same results for the spectroscopic factors of the main particle and hole shells close to the Fermi 
energy. These correspond to 0.775 for pi/2 and 0.766 for P3/2 as reported in Tab. El The introduction 
of correlations reduces these values and brings them down to 0.745 and 0.725, respectively. This result 
reduces the discrepancy with the experiment by about 4% and shows that collectivity beyond the TDA 
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RPA 


1st itr. 


) in /~\ 1 T v 

ZllCl ILI . 


-4 V*H 1 "TV 


^Lll ILl. 




0.775 


0.745 


0.775 


0.777 


0.774 


0.776 




0.766 


0.725 


0.725 


0.727 


0.722 


0.724 






0.015 


0.027 


0.026 


0.026 


Total occ. 


14.56 


14.56 


14.56 


14.57 


14.58 


14.63 



Table 3: Hole spectroscopic factors (Za) for knockout of a £ = 1 proton from ^^O. These results refer 
to the initial IPM calculation in TDA and RPA and to the first four iterations of the DRPA equations. 
All the values are given as a fraction of the corresponding IPM value. Also included is the total number 
of nucleons, as deduced from the complete one-hole spectral function, for each iteration. 

level is relevant to explain the quenching of spectroscopic factors. Since the present formalism does not 
account for center-of-mass effects, the above quantities need to be increased by about 7% before they 
are compared with the experiment I197j . 

The RPA results have been iterated several times, with the aim of studying the effects of fragmen- 
tation on the RPA phonons and, subsequently, on the spectral strength. Since one is mainly interested 
in the low-energy excitations, it is sufficient to keep track only of the strongest fragments that appear 
— close to the Fermi energy — in the (dressed) sp propagator, Eq. (jHl), while the residual strength is 
collected in two effective poles, above and below the Fermi level, for each orbital |192t I196j . As shown in 
Tab. El only a few iterations are required to reach convergence for the values of the spectroscopic factors. 
The converged distribution of one-hole strength is compared with the experimental one in Fig. 071 The 
main difference between these results and the one obtained by using an IPM input is the appearance 
of a second smaller P3/2 fragment at -26.3 MeV. This peak arises in the first two iterations and appears 
to become stable in the last one, with a spectroscopic factor of 2.6%. This can be interpreted as a peak 
that describes the fragments seen experimentally at slightly higher energy. This result corresponds to 
the first time that such a fragment is obtained in calculations of the spectral strength. 

A second effect of including fragmentation in the construction of the RPA phonons is to increase the 
strength of the main hole peaks. The pi/2 strength increases from the 0.745 obtained with IPM input to 
0.776, essentially canceling the "improvement" gained by the introduction of RPA correlations over the 
TDA ones. The main peak of the P3/2 remains at 0.722 but the appearance of the secondary fragment 
slightly increases the overall strength for this orbital as well. This behavior is due to the fact that 
redistribution of the strength (and principally, the reduction of the main spectroscopic factors) tends to 
screen the nuclear interaction, with respect to the IPM case. The consequence is a reduced effect of RPA 
correlations when fragmentation is included in the construction of the phonons. This feature has also 
been observed in other self-consistent calculations of the sp spectral strength, for example in nuclear 
matter j72|. Obviously, this makes the disagreement with experiments a little worse and additional 
work is needed to resolve the disagreement with the data. Nevertheless, it is clear that fragmentation 
is a relevant feature of nuclear systems and that it has to be properly taken into account. 

Together with the main fragments, the Dyson equation produces also a large number of solutions 
with small spectroscopic factors. This strength extends down to about -130 MeV for the one-hole 
spectral function and up to about 100 MeV for the one-particle case. This background represents the 
strength that is removed from the main peaks and shifted up or down to medium missing energies. 
Another 10% of strength is moved to very high energies due to SRC Its location cannot be 

explicitly calculated in the present approach but the effects on the reduction of spectroscopic factors at 
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Figure 47: One-proton removal strength as a function of the hole sp energy = Eq — -E"^"^ 
for ^^O for angular momentum £ = 1 (left) and £ = 0, 2 (right). For the positive parity states, 
the solid bars correspond to results for ^5/2 and ^3/2 orbitals, while the thick lines refer to 
si/2- The top panels show the experimental values taken from |198j . The middle panels give 
the theoretical results for the self-consistent spectral function. The bottom panels show the 
results obtained by repeating the S"^*^ iteration with a modified ph-DRPA spectrum, in which 
the lowest eigenstates have been either removed or shifted to the corresponding experimental 
values. 



low energy is accounted through the energy dependence of the G-matrix. Accordingly the total number 
of nucleons deduced from the fully self-consistent spectral function is 14.6 as reported in Tab. El This 
gives an estimate for the overall occupancy of high-momentum states of about 10%, in agreement with 
direct calculations |183( ll5Uj as discussed in Sees. 15.21 and 15.31 

We note that the present approach also yields predictions for the strength distribution of quasipar- 
ticle excitations. Although little experimental information is available for the addition of a nucleon (as 
well as for the removal of a neutron). Other knockout and nucleon transfer reactions, employing nuclear 
probes, have been recently been developed or reanalyzed to extract hole spectroscopic factors of both 
stable and radioactive nuclei \2'2'6\ I224j . These tools may yield better accuracies than in the past and 
can, in principle, be used to extract information on the one-particle spectral distribution. 

A deeper insight into the mechanisms that generate the fragmentation pattern can be gained by 
investigating directly the connection between the spectral function and some specific collective states. 
To clarify this point we repeated the third iteration using exactly the same input but without keeping 
the lowest ph 0+ state at its experimental energy, it was discarded instead. The resulting p hole spectral 
function is shown in the lower- left panel of Fig. ^| In this calculation no breaking of the P3/2 strength 
is obtained but only a single peak is found with a spectroscopic factor equal to 0.75, which is the sum 
of the two fragments that are obtained when the O2 state is taken into account. This result can be 
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Figure 48: Examples of contributions involving the coupling of two independent ph phonons. 
In total, there are sixteen possible diagrams of this type, obtained by considering all the pos- 
sible couplings to a ph state. The two-phonon ERPA equations sum all of these contributions 
in terms of dressed sp propagators. 



interpreted by considering the P3/2 fragments as a hole in the ground state and in the first excited O"*" 
state of the ^^O core, respectively. The correct fragmentation pattern can be reproduced only when the 
latter two levels are close enough to one another in energy, so that the two configurations can mix. 

The other two low-lying states of ^^O that may be of some relevance are the isoscalar 1~ and 3~. 
These excitations are reproduced reasonably well by RPA type calculations [TBlj but are typically found 
at ~3 MeV above the experimental results. The lower panels of Fig. |^ show the results for the even 
parity spectral functions that are obtained when both the 3~ and the 1~ ph-DRPA solutions are shifted 
to match their experimental values. In this case, a ^5/2 hole peak is obtained at a missing energy 
-17.7 MeV, in agreement with experiment, and with a spectroscopic strength of 0.5%. It is interesting 
to note that the shifting of the 3~ and 1~ collective states does not produce any other noticeable change 
in the theoretical spectral function. 

5.6 Spectrum of ^^O and extension of RPA 

The above results suggest that the main impediment for further improvements of the description of the 
experimental spectral strength is associated with the deficiencies of the RPA (DRPA) description of the 
excited states. One important problem is the first excited state, with quantum numbers 0^, which is 
seen to produce important effects on the fragmentation of the p shell. The DRPA result do find this state 
at low energy but can describe at most one collective phonon for a given angular momentum, parity, 
and isospin combination J'^, T whereas several low-lying isoscalar 0''' and 2+ excited states are observed 
at low energy in ^^O, as well as additional 3~ and 1~ states. A possible way to proceed would be to first 
concentrate on an improved description of the collective phonons by extending the RPA to explicitly 
include the coupling to two-particle— two-hole (2p2h) states. Such an extended RPA procedure has 
been applied in heavier nuclei with considerable success |225| 122b] . In order to be relevant for ^^O, this 
approach requires an extension in which the coherence of the 2p2h states is included in the form of the 
presence of two-phonon excitations |227j . This corresponds to adding two-phonon configurations to the 
kernel K^^ of the BSE for the ph propagator ()21|). as shown in Fig. HHl ,199^. 

The formalism to include two-phonon configurations has been presented in Ref. |199j . where it is 
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referred as two-phonon extended RPA (ERPA). In that work the ph-DRPA equation has been solved 
first, using the self-consistent sp propagator discussed in Sec. 15.51 The lowest DRPA solutions for 
both the 0+, 3~ and 1~ channels were shifted down to their respective experimental energies and then 
employed to generate the two-phonon contributions for the ERPA calculation. The spectrum obtained 
for ^^O is displayed in Fig. Upland Tab. HI together the total ph strength of each state, 

Zn. = E . (159) 

a/3 

Tab. m also reports the relative strength of ph and two-phonon admixtures in the wave function for 
each solution of the ERPA equation. Due to the screening effects associated with dressing the sp 
propagators, the solution for the first isoscalar 0"^ state in DRPA is found much higher in energy 
at ~17 MeV. It is important to realize that this solution has a sharp ph character and therefore it 
cannot be identified with the experimental state. The latter state is known to be dominated by Ahu 
configurations |229[ 1^0, 215j . On the other hand inelastic electron scattering experiments |231j clearly 
excite this state. From the point of view of the SCGF approach, the one-body response is completely 
described by the polarization propagator, Eq. (|^. and therefore, the total one-body strength must be 
represented by Zn„ ()159|) . This indicates a strong coupling to the (dressed) ph propagator ((221) • ^^e 
basis of this similarity, it still seems plausible to shift the ph-DRPA solution down in energy in order 
to study its mixing with other configurations. 

The corresponding results yield an isoscalar 0"*" state with a predominant ph character at ~17 MeV, 
as it was found in DRPA. Although this solution is now characterized by a partial mixing with two- 
phonon configurations. Tab.EJshows that this is the result of mixing with the lowest solution in the same 
channel, which ends up at ~11 MeV. The latter is predominately a two-phonon state. It is also seen 
that in both cases the relevant configuration comes from the coupling of two 0^ phonons themselves, 
in Ref. |199j it was found that the wave functions for these two states contain several relevant ph 
configurations, obtained from different quasiparticle fragments in the pf {sd) shells combined with 
quasihole fragments of the p (s) shells. Therefore, the situation is more complicated than the simple 
picture of only two levels interacting with each other. 

The calculations of Ref. |232j have shown that the bulk of the 4p4h contributions to the first excited 
state of ^^O may come from the coupling of four different phonons with negative parity (3~ and 1^). 
However, the self-consistent role of coupling positive parity states was not considered in that work. 
Both this effect, the RPA correlations and the inclusion of the nuclear fragmentation allow for the -at 
least partial- inclusion of configurations beyond the 2p2h case even when no more than two-phonon 
coupling is considered, as discussed here. Still it would be helpful to further study the importance of 
three- and four-phonon configurations within this apporach. 

The low-lying negative parity isoscalar states (3~ and 1~) are only slightly affected by two-phonon 
contributions and remain substantially above the experimental energy at 9.23 and 10.90 MeV, respec- 
tively, and more correlations will be needed in order to lower their energy. We note that these wave 
functions contain a two-phonon admixture obtained by coupling the low-lying O"*" excitation to either the 
3^ or 1^ phonons. According to the above interpretation of the first O"*" excited state this corresponds 
to the inclusion of 3p3h and beyond. At the same time two additional negative parity solutions, with 
two-phonon character, are found at higher energy in accordance with experiment. These configurations 
are not reproduced by the simple DRPA. 

The two-phonon ERPA approach also generates a triplet of states at about 12 MeV, with quantum 
numbers 0^, 2^ and 4+. A similar triplet is found experimentally at 12.05, 11.52 and 11.10 MeV, which 
also corresponds to twice the experimental energy of the first 3~ phonon. The ERPA solutions for 
this triplet are almost exclusively obtained by 3~ 3~ states and therefore have a 2p2h character, in 
accordance with Refs. |23H I234[ I235j . Further interaction in the pp and hh channels will be needed in 
order to reproduce the correct experimental splitting and experimental strength ^232^ I233j . 
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Figure 49: Results for the DRPA and the two-phonon ERPA propagator of ^^O with the 
dressed input propagator computed in Ref. |196j . middle and last column respectively. In 
solving the ERPA equation, the lowest 3~, 1~ and O"*" levels of the DRPA propagator where 
shifted to their experimental energies. All other DRPA solutions were left unchanged. The 
excited states indicated by dashed lines are those for which the (E)RPA equation predicts 
a total spectral strength Zn„ lower than 10%. The first column reports the experimental 
results fh^ . 



5.7 Results for the ^^0{e, e'pp) cross section 

Data from two-nucleon knockout reactions have been measured recently at NIKHEF 12361 12371 1^^^^ 
and MAMI |239t l24Uj facilities for different targets. Theoretically, these cross sections can be computed 
with the models developed by the Pavia |24H I242j and Ghent [243^ groups for the emission of both two 
protons or a proton-neutron pair. Calculation for the ^^0{e,e'pp) case have shown that the transition 
to the ground state of ^^C is dominated by the high-momentum components in the two-nucleon overlap 
function. Comparison with the data from NIKHEF for this reaction have provided a clear signature of 
SRC effects |2SZ1 12111 • 

An important ingredient in the description of two-nucleon knockout reactions is the two-hole spectral 
function. This quantity can be extracted from the imaginary part of the two-body propagator g^^{u), 
compare Eqs. (j27|) and (j54|) . and contains the nuclear structure information on the correlated pair before 
the reaction takes place. The propagator g^^iuj) can be obtained by solving pp(hh) RPA equations as 
done in Ref. |245j . However, typical calculations of this kind yield results within a model space that is 
too small to contain the short-range effects. In these cases one needs to compute the high-momentum 
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Table 4: Excitation energy and total spectral strengths obtained for the principal solutions of DRPA and 
two-phonon ERPA equations. The total contribution of ph and two-phonon states of the ERPA solutions are 
shown. For states below 15 MeV, the columns on the right side give the individual contributions of all the 
relevant two-phonon contributions. 

components separately. In the following we outline the approach taken in Refs. |245| I241j and show a 
few results for the ^^0{e,e'pp) cross section. 

The sp spectral function of ^^O obtained in Ref. |195j was employed as the starting point of a 
hhDRPA calculation (see Eq. (j29j) ), eventually obtaining a two-hole spectral function of the form |245j . 

5'jr(Pi,P2,Pi',P2',u;) = 

E E Y.'^:i''ip^'^P2'){Xc,jTrx'^,j^<i>i^^ (160) 

k abcdeV MT3 

where J and T refer to the total angular momentum and isospin of the residual nucleus and M andTs 

______ JT jA. 2 

are their projections. In Eq. p60|l . the energies i?^ ' are the the DRPA solutions for excitation 
energies of the residual nucleus and the quantities 

X'aUT-i'^j'^-'llMjlK) (161) 

represent the components of the solutions for the two-body overlap amplitudes. The sums extend 
only over those orbitals that are occupied in the model space V and the spectral function S^^{uj) 
is transformed into the momentum representation by employing two-body harmonic oscillator wave 
functions $^^*^(pi, P2). 

The description of the high-momentum components due to SRC requires the inclusion of a very 
large number of basis states, at least up to lOO^u; in a harmonic oscillator basis |222j . However, 
the description of long-range correlations by solving a Bethe-Salpeter or an RPA equation ()29|) is not 
feasible within such a large space. For this reason the complete basis is split into a model space V, 
and a complementary space Q = 1 — V. For ^^O the space V was chosen to contain the orbitals up to 
the pf shell only, which is large enough to accommodate long-range effects. The justification for this 
procedure is that SRC are caused by close encounters of two nucleons, which mainly depend on the 
nuclear density and therefore are not very sensitive to details of the long-range structure. The latter, on 
the other hand, may be calculated within the space V with a suitable effective interaction in which the 
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SRC are incorporated at least in ladder approximation. This effective interaction can then be obtained 
by following Brueckner's individual pair approach j24ti| I247ij by solving the Bethe-Goldstone equation 
(BGE). Using the technique of Ref. |2224 the equation for the correlated pair wave function 

l^afe) = |<fab) + — ^V^l^af,) , (162) 

is solved in the finite nucleus. In this equation the Pauli operator Q prohibits scattering into orbits of 
the finite shell-model space P, in which long-range correlations are trated in an DRPA approach. In 
Eq. ()162|) . represents the same uncorrelated shell model wave functions appearing in Eq. ()16U|) . u is 
the propagation energy of the pair and Hq is the Hamiltonian without residual interaction. In this case, 
the energy uj refers to the propagation of two holes which precludes the vanishing of the denominator 
in Eq. (fT^ . 

From the solution of Eq. ()162|) one obtains the defect wave function as the difference between the 
correlated and uncorrelated pair wave functions 

\X) = 1^) - l<f) , (163) 

which represent the components of the two-body wave function that are orthogonal to the space P, and 
the G matrix as the effective interaction in V according to 

{^'I\G\^it') = {^ir\v\^in ■ (164) 

The essential step taken here, is to approximate the spectral function p60|) by the expression 
5'jt(Pi,P2,Pi',P2',u;) = 

E E E^c/^'(Pl^P20(X.W)*^aW^r(Pl,P2)5(o.-(^^^ (165) 

k abcd&V MT3 

where the summation over orbits is still limited to the finite shell-model space V but the uncorrelated 
wave functions are replaced by the correlated ones (see Eq. ()l(i2|l ). This is in line with the argument just 
given, that hard binary collisions, treated in the BGE and giving rise to high-momentum components, 
proceed independently of the long-range correlations. The latter are taken into account in the shell 
model amplitudes X within the limited space V. 

In the spectator model with a plane-wave approximation for the outgoing protons, the contribution 
of the spectral function to the cross section is given by the superposition S{pi,p2, E), Eq. This 
function is plotted in Fig. EOl at the kinematics corresponding to the ^^0{e,e'pp) measurements at 
NIKHEF |236j . for the transitions to the 0^ ground state and to the excited 2"*" levels of ^^C. In these 
kinematics the momenta pi and p2 of the protons are almost antiparallel to each other and coplanar 
with the momentum q transferred by the electron. As a consequence, high and equal values of their 
magnitudes correspond to high relative momenta in the plots of Fig. EHl The effects of including the 
SRC is shown in the lower plots, where the correlation functions ()163|) have been computed with both 
the Reid-Soft-Core and Bonn-A potentials. Figure EHl also shows the plots labeled 'No SRC that 
correspond to neglecting the correlation functions and computing the spectral functions with Eq. (|16Up . 
As expected, one observes that the latter gives no contribution at higher momenta. The high-momentum 
part of the spectral function is about a factor of two larger for the Reid-Soft-Core potential than for 
the Bonn-A potential in the momentum range around 3-5 fm~^. The short-range correlations give rise 
to a spectral function which is clearly distinct from the one without SRC. 

The Ghent-model |243j accounts for SRC in a different way. There, single-particle wave functions 
obtained from Hartree-Fock calculations are used and the many-body wave function, is approximated 
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Figure 50: Superposition of spectral functions appropriate for the removal of two protons 
with final momenta pi and p2 from ^^O leading to the final 0"'' ground state or first excited 2"*" 
state in ^^C. The plots are given for a kinematic setting used in experiments at NIKHEF |236j . 
The momentum vector q is fixed along the z-axis, with length 313 MeV/c. The momenta 
Pi and p2 are in the same plane with q at —49 and 123 degrees angles with respect to 
this transferred momentum, respectively. The upper plots correspond to harmonic oscillator 
wave functions without the inclusion of short-range correlations (see Eq. ()160p ). In the 
lower plots these SRC are incorporated by the defect functions computed the Reid and the 
Bonn-A potential (Eqs. (fT?K^ and (fTHK|l ). 
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Figure 51: The cross sections for the transition to the ground state in ^^C as a function 
of the missing momentum, measured for the reaction ^^0(e, e'pp)^^C at three values of uo. 
The curves are obtained from calculations of Refs. |241j (top) and |243j (bottom). The solid 
curves represent the calculated cross sections; the dashed and dotted curves correspond to 
the contributions of the one- and two-body hadronic currents, respectively. 



to first order in the correlation function, 

A 

^(n, . . . r^) = n (1 - 9{rij)) $(ri, . . . va)/^ , (166) 

i<j=i 

where is a normalisation factor. The effects of SRC are included in the factor (1 — g{rij)) by means 
of the correlation function g{rij) (which zero in the absence of correlations). The calculations presented 
here employed the g{rij) of Ref. |64 , also shown in Fig. EHl 

Figure |^ compares the cross sections for the calulations of Ref. |24H I243j with the experimental 
data, as a function of the missing momentum at the three values of uj |244j . For this calculation 
final state interactions for the outgoing protons were taken into account by using distorted waves but 
neglecting their mutual interaction. A pair of (pi/2)^ or (^3/2)^ protons can be coupled to total angular 
momentum J=0, leaving the ^'^C nucleus in the ground state, through either a ^Sq or a ^Pi relative 
state. The coupling scheme used in the calculation is as follows. The ^Sq state is always associated 
with an angular momentum of the center-of-mass motion L = 0, and the ^Pi state always with L = 1. 
From Fig. 1^ it is clear that for all three values of the energy the missing-momentum dependence of the 
measured cross sections is similar. In Refs. |23H I237j , it has been pointed out that such a momentum 
distribution reflects an angular momentum L = for the center-of-mass motion of the pair, and thus 
suggests a dominant role for the knockout of a pair driven by SRC. The theoretical cross sections 
are represented in Fig. |^ by the solid curves; the contributions of the one- and two-body currents 
are given by the dashed and dotted curves, respectively. The calculated cross sections agree well with 
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the data at all three values of uj. The curves, representing the contributions of the one- and two-body 
currents to the ground-state transition, indicate that at u;=180 and 210 MeV the reaction is dominated 
by one-body currents, and that the contribution of two-body currents increases with increasing energy 
transfer. Conceptually both models are quite similar, and neither of the two contains free parameters. 
The effects due to SRC, though, are accounted for in different ways as discussed above. The Bonn-A and 
Reid Soft Core potentials, adopted in the calculations, are both realistic NN-potentials and successfully 
used in many nuclear-structure calculations. Hence, the data agree with the theoretical results obtained 
independently with the two models. Furthermore, both models predict that the largest contribution to 
the cross section stems from one-body hadronic currents driven by SRC. This justifies the conclusion 
that evidence is obtained for SRC although the data do not distinguish details of the actual form of 
the SRC. Future experiments will include the removal of proton-neutron pairs and therefore may also 
yield information on tensor correlations |242[ 1248] . An improved description of the mutual interaction 
between the outgoing protons |^9 , 250j will help in further developing the usefulness of this reaction 
to elucidate details of short-range nuclear dynamics. 

5.8 Extraction of spectroscopic information from (e, e']?) reactions at high 

The analysis of the (e,e'p) data has relied on the Distorted Wave Impulse Approximation (DWIA) 
for both the Coulomb distortion of the electron waves (in heavy nuclei) and the outgoing proton ^01 
Uni I25H I252j . Some aspects of this analysis have been discussed in Sec. El The proton distortion is 
described in terms of an optical potential required to describe elastic proton scattering data at relevant 
energies ^3]. There is some uncertainty related to this treatment since elastic proton scattering is 
considered to be a surface reaction and no detailed information is obtained related to the interior of 
the nucleus. This uncertainty gives rise to an estimated error of about 10%. Such an estimate may 
be inferred by considering the difference between the relativistic and nonrelativistic treatment of the 
proton distortion. It is shown in Ref. |253j that this difference is essentially due to the reduction of 
the interior wave function in the relativistic case. This feature can also be generated by including a 
reasonable amount of nonlocality in the optical potential |253j . 

A serious challenge to the interpretation of (e,e'p) experiments was recently published in Refs. |254[ 
I255j . This challenge consists in questioning the validity of the constancy of the spectroscopic factor 
as a function of the four-momentum (squared), Q^, transferred by the virtual photon to the knocked- 
out nucleon. In Ref. ||25i] a conventional analysis of the world's data for (e,e'p) experiments on ^^C 
at low generated results consistent with previous expectations. Data at higher values of were 
then analyzed within the framework of a theoretical model which employs Skyrme-Hartree-Fock bound- 
state wave functions for the initial proton, a Glauber-type description of the final-state interaction 
of the outgoing proton, and a factorization approximation for the electromagnetic vertex |2541 1255j . 
Within the framework of this theoretical description, spectroscopic factors were obtained which increase 
substantially with increasing for the ^^C nucleus. 

The spectroscopic factor defined in Sec. |21 is a many-body quantity defined without reference to 
a probe. The discussion in Sec. |21 shows that this spectroscopic factor is obtained by comparing a 
calculation of the cross section with the actual data. The wave function of the removed nucleon is 
the solution of a Woods-Saxon potential at the appropriate removal energy. This potential is adjusted 
to generate an optimum description of the shape of the experimental cross section. This resulting 
theoretical representation of the cross section is then multiplied by a constant factor to coincide with 
the experimental cross section. This constant factor is then interpreted as the spectroscopic factor. 
Another important ingredient in this analysis is the choice of the electron-proton cross section which 
must be considered off-shell j256] . This leads to a small additional uncertainty in the analysis of low 

data as discussed recently in Ref. |257j . 
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Further clarification of tliis intriguing situation witli different spectroscopic factors at different 
lias been obtained in Ref. |258j . In this paper a study of recently published (e,e'p) data for ^^O at 

= 0.8 (GeV/c)^ |259j has been presented. An unfactorized approach was used, as it is required at 
low and as it has been recently advocated also for high reactions in Ref. |26Uj . However, the 

higher energy of the outgoing proton requires a description which contains different elements than the 
conventional low analysis. For the electromagnetic current operator the approach of Ref. |261j was 
followed, where a relativistic current operator was used in a Schrodinger-based calculation, avoiding 
any nonrelativistic reduction and including the effect of spinor distortion by the Dirac scalar and vector 
potentials. As for the final state, a recently developed eikonal description of the final-state interaction 
(FSI) of the proton with the nucleus that has been tested against DWIA solutions of a complex spin- 
dependent optical potential |262l I263j was employed. The absorption of the proton was described 
theoretically by linking it to the corresponding absorption of a nucleon propagating through nuclear 
matter. The relevant quantity is the nucleon self-energy which is obtained from SCGF calculations of 
nucleon spectral functions discussed in Sec. 14.61 This description of FSI was combined with previous 
results for the bound-state wave functions of thep-shell quasihole states in ^^O [150J discussed in Sec. 15.21 
In Ref. |2U9j these wave functions have been used to analyze low data for the ^^0(e,e'p) reaction |198j . 
These very same wave functions produce a good description of the shape of the coincidence cross section 
for the p-shell quasihole states at high |259j using the same spectroscopic factors obtained from the 
low g2 data [121. 

We first establish the quality of the original assertion pertaining to low by referring to Fig. 021 
for the p3/2 state in ^^O. A similar result is obtained for the pl/2 state |258j . The corresponding 
spectroscopic factors to obtain agreement with these low data are given by 0.644 and 0.537 for 
the pl/2 and p3/2 states, respectively. In Fig. 1^21 the same reaction is considered in a very different 
kinematical regime, namely at constant {q,uj) with = 0.8 (GeV/c)^ |259j . The data here refer to 
a five-fold differential cross section, avoiding any ambiguity in modelling the half off-shell elementary 
cross section aep. Results for the transition to the ground state pi/2 have been multiplied by 20. The 
theoretical calculations are displayed with solid lines for the results with the quasihole bound state and 
dashed lines by employing the wave function of Ref. |264j . The preference for the first choice is evident. 
In any case, it is remarkable that the calculations reproduce the data by using the same spectroscopic 
factors as in the previous kinematics, i.e. 2'opi/2 = 0.644 and 2'op3/2 = 0.537. Therefore, contrary 
to the findings of Ref. |255j . no need for a dependence of the spectroscopic factors over a wide 
kinematical range has been established. This is also in agreement with the results obtained by a mean- 
field description in the context of relativistic DWIA |265j . This outcome is reassuring and is further 
confirmed in Ref. |266j . In this paper it is also shown that a similar statement holds when the analysis 
of the data is performed using a relativistic description of the distorted proton waves combined with 
a corresponding mf description of the bound states. Also in this case, the slightly larger spectroscopic 
factors of 0.708 and 0.602 for the pl/2 and p3/2, respectively describe both the low and high data. 
The difference between the spectroscopic factors for these different analyses is about 5% which may 
be considered as a measure of the uncertainties associated with the assignment of spectroscopic factors 
until an improved microscopic desciption of the FSI becomes available. 

Finally, one may conclude that the treatment of the bound-state wave function is not responsible for 
the dependence found in Ref. |255j . From this observation one may infer that it is useful to extend 
the analysis of the high data to other nuclei using the eikonal description supplemented with the 
nuclear matter damping description while foregoing the use of microscopic quasihole wave functions. 
Clearly an analysis of the data considered in Refs. |254j and |255j will further clarify the validity of the 
present analysis. It is already encouraging to note that the damping employed in the calculations is 
adequate to describe the observed absorption in the NE18 experiment |267j as discussed in Ref. pSBj . 
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Figure 52: Cross section for the ^^0(e,e'p)^^N reaction at = 0.8 (GeV/c)^ in perpendicular 
kinematics [259J. Data for the pl/2 state have been muhiphed by 20. The sohd hues represent 
the resuh of the calculation [258^ . The dashed lines are obtained by replacing the quasihole 
states with the bound state wave functions of Ref . |264j . In all cases, the results have been 
rescaled by the same spectroscopic factors as used for the NIKHEF kinematics (see Fig. 021 
for example), namely Z=0.644 and Z=0.537 for the pl/2 and p3/2 states, respectively. 



5.9 Recent experimental data at higher missing energies 

The bulk of experimental data obtained during the last decades by means of the (e,e'p) reaction have 
focussed almost exclusively on the shells in the immediate vicinity of the Fermi energy, especially for 
heavier nuclei. As a consequence the amount of protons seen experimentally has been considerably 
smaller than the total number of protons Z inside the nucleus. 

More recently, a welcome addition to our knowledge of the sp strength distribution has been provided 
by the analysis of one of the last (e, e'p) experiments performed at NIKHEF. In this experiment the 
spectroscopic strength in a large range of missing momentum (up to 270 MeV/c) and missing energy 
(up to 100 MeV) was measured for 208p|^ |143| I144j . The analysis of this experiment for the first time 
identifies the occupation numbers of all the mean-field protons in this nucleus. These results are shown 
in Fig. inland compared with a theoretical calculation for nuclear matter by Benhar et al. [147^ . The 
analysis for these deeply bound levels yields a global depletion of 15% for all these proton orbits. This 
confirms the expectation that nuclear matter results yield a reliable estimate of the depletion due to 
SRC [nni- Recent theoretical calculations of second generation spectral functions in nuclear matter [72] 
confirm the robustness of this observation as shown in Fig. 1201 These new results have been obtained 
from self-consistent Green's functions which incorporate full off-shell propagation and include the effects 
of short-range and tensor correlations as discussed in Sec. 14.61 

The discussion in Sees. 15.21 and 15.31 for the case of finite nuclei, emphasized another consequence 
associated of the presence of short-range correlations in nuclear systems. Besides the sizable amount of 
strength that is removed from the mean-field region, these same correlations also supply a corresponding 
amount of strength in the ground state in the form of high-momentum components. A recent experiment 
at JLab, more fully discussed in Refs. |187| 1269] , has been able to identify these correlations for different 
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Figure 53: Experimental occupation number for the all the mean field shells of ^osp]-, ji43j 
The dotted line gives the value of occupation probability n^^^ computed for nuclear mat- 
ter [TiTj . 



nuclei at high missing energy |186j . Figure |^ shows first results from this experiment for the spectral 
function of ^^C. The location of these high-momentum components conforms with the mechanism that 
admixes these correlations with 2hlp states. This can only be accommodated at such high missing 
energies due to the following observation. Two-hole states in the nucleus tend to have a small total 
momentum. Momentum conservation then requires that a high momentum admixture is accompanied 
by a particle with equal and opposite momentum inside the corresponding self-energy contribution. The 
resulting energy for the presence of these high-momentum components is thus given by {s2h) — /'2m. 
These missing energies can therefore easily be as large as 100 to 200 MeV. As one can see from the 
picture, the maximum of the theoretical spectral function follows this expectation (indicated by the 
arrows) but the experiment clearly finds the maximum at slightly smaller energy with respect to the 
naive application of the above argument. A spreading of the peak in energy when moving to high 
momenta appears to be in agreement with the trend shown by Fig. |^ 

The spectral function at high momenta has been computed in finite nuclei with the approach of 
Sees. 15.21 and 15.31 jl83l 1130112141 and in local density approximation |179j . The experimental momentum 
distribution obtained for carbon is found to be in between these theoretical predictions |187j . We note 
that one limiting issue in taking measurements at higher missing energies comes from the necessity 
of transferring a large amount of energy from the probe to the struck nucleon. The outgoing proton 
in an experiment that probes this high-energy region can have kinetic energies of several hundreds of 
MeV, leading to strong rescattering effects, excitation of nucleon resonances, etc. In certain kinematical 
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Figure 54: Reduced spectral function for ^^C in the short-range correlated region obtained 
with (e, e'p) reaction |187j . At missing energies above 150 MeV, sizable contributions from 
the excitation of a A resonance are seen. The arrows indicate indicate the location of the 
maximum of the spectral function expected according to a simplified picture. 



conditions the associated final state effects can overwhelm the direct process. This has been the case for 
some previous experiments |27Ul 12711 1272j . The issue of how to minimize the FSI has been addressed 
in Ref. |185j . There, it is shown that FSI in exclusive (e, e'p) cross sections are dominated by two-step 
processes that become particularly relevant when perpendicular kinematics are employed to probe the 
regions of small spectral strength. A study of the kinematic conditions shows that for perpendicular 
kinematics the rescattered nucleons move spectral strength in the k-E plane, from the top of the 
ridge toward regions where the correlated strength is small, therefore submerging the direct signal in 
a large background noise. In Ref. |185j . it was suggested that the contribution of rescattering can 
be diminished using parallel kinematics. The Jefferson Lab data discussed here was therefore taken in 
these conditions p[B"6l ll87j . First calculations of the rescattering process appear to confirm the expected 
trend with respect to the kinematics \2bH\ I2ti9j 

We also note that other important contributions to the yield of this experiment are expected to come 
from (e, e'A) reactions followed by the decay of the A resonance. The broad peak due to the excitation 
of the A is clearly visible on the right hand side of Fig. |^ For heavier nuclei or when perpendicular 
kinematics are used, the above rescattering effects are seen to become more and more relevant and 
at the same time the strength from the delta region tends to come at low missing energy, filling the 
gap that separates it from correlated region. Therefore a proper treatment of final state interactions 
becomes necessary to study the behavior of the strength in the short-range region when the size of the 
nucleus is increased j269j . First calculations in this direction have been reported in Refs. |2ti8| I269j . 

The relevance of the experiments discussed in this section lies in the fact that they extend the exper- 
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imental knowledge of correlations to the regions of the momentum-energy plane that were previously 
inaccessible. Being able to identify high-momentum components [187| in addition to locating all the sp 
strength associated with the mean-field orbits [143^ .144|j completes the search for the missing protons. 

6 What a proton does in the nucleus and comparison with 
other fermion systems 

From the results discussed sofar it is clear that a crucial role in obtaining pertinent information about 
the properties of protons in the nucleus has been played by coincidence experiments involving electron 
beams. In the following we will therefore summarize on what has been learned in general terms about 
the properties of these protons. Before embarking on this overview, it is important to note that recent 
work on single-nucleon knockout with fast radioactive beams j224j proposes to extend this information 
for valence protons to neutrons and unstable nuclei. Some promising work in this direction has been 
reported in Refs. |273| I274j . Indeed, there is a long tradition with hadronic probes to study spectroscopic 
factors in nuclei. Due to the inherent complexity of hadron-induced reactions it has been more difficult 
to establish absolute spectroscopic factors from corresponding experimental data. Nevertheless, it is 
possible to generate a consistent analysis of (e,e'p) and (d,^He) experiments as shown in Ref. |223j . 
The original discrepancies between these different experiments disappear when a new analyis of the 
(d,^He) experiments is performed using finite-range DWBA procedures combined with bound-state 
wave functions obtained from the analysis of (e,e'p) reactions. The resulting spectroscopic factors then 
appear to be quite consistent with the spectroscopic factors obtained from the (e,e'p) reaction |223j . 
This observation demonstrates that with proper care it is possible to generate valuable information from 
hadron-induced nucleon knockout experiments. It may therefore be possible to extend the extraction 
of spectroscopic factors to nuclei far off stability. This exciting new development will allow the study 
of the properties of nucleons in the nuclear medium in different regions of the periodic table and may 
provide new challenges for our theoretical understanding. 

The present situation concerning the knowledge of the properties of protons in the nucleus may 
be summarized as follows. The consequences of short-range correlations in nuclear systems on the 
properties of protons appear to be theoretically well understood and becoming available for experimental 
scrutiny as discussed in Sec. El These results show that the effect of SRC is two-fold. First, it involves 
the depletion of spectroscopic strength from the mean-field domain as discussed in Sec. 15.91 for the 
experimental data obtained for ^°'^Pb. These data show that this depletion in heavy nuclei corresponds 
to about 15% for all the deeply bound proton levels. This result was predicted more than ten years ago 
based on nuclear-matter calculations briefiy reviewed in Sec. 14.21 For lighter nuclei all theoretical work 
suggests that this amount may be closer to 10% as discussed in Sec. 15.21 Accompanying this information 
is the realization that valence shells near the Fermi energy will not contain substantial amounts of high- 
momentum components. This has been experimentally confirmed and clarifies the other role played by 
SRC in nuclei, i. e. the admixtures of high-momentum components at high missing energy to account 
for the missing protons removed from the mf location. Recent attempts to quantify the amount of 
this strength and the corresponding location have been presented in Sec. 15.91 It now appears that 
about the right amount of high-momentum strength can be identified in the expected region although 
additional work is needed to complete the analysis. The location of these high-momentum components 
nevertheless broadly conforms with the mechanism that admixes these correlations with 2hlp states at 
large missing energies as discussed in Sec. 15.91 

Being able to identify high-momentum components in addition to locating all the sp strength associ- 
ated with the mean-field orbits |143| I144j completes the identification of the properties of protons in the 
ground state of the nucleus. The latter understanding is illustrated in Fig. EH Several generic diagrams 
are identified in Fig. EHl which have unique physical consequences for the redistribution of the single- 
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particle strength when they are taken into account in the solution of the Dyson equation as discussed 
in this paper. The middle column of the figure characterizes the mean-field picture that is used as a 
starting point of the theoretical description. The right column identifies the location of the sp strength 
of the orbits just below the Fermi energy when correlations are included. One may apply this picture 
for example to the 3si/2 proton orbit in ^°^Pb. The physical mechanisms responsible for the correlated 
strength distribution are also identified. The strength of this orbit remaining at the quasihole energy 
is about 65%. Long-range correlations are responsible for the loss of 20% of the strength due to the 
coupling to nearby 2plh and 2hlp states. This loss is symmetrically distributed above and below the 
Fermi energy and is physically represented by the coupling to low-lying surface modes and higher-lying 
giant resonances. The resulting occupation number of this orbit therefore corresponds to 75%. More 
deeply bound nucleons have higher occupation numbers corresponding to 85%. As discussed in Sec. 15.91 
this is true for all these deep-lying orbits and is consistent with a global depletion due to SRC of 15%. 
The corresponding location of this strength is identified at very high energy in the particle domain 
and is due to the short-range and tensor correlations induced by a realistic nucleon-nucleon interaction 
as discussed in Sec. |3] In the left column of Fig. EHl the generic diagram that is responsible for the 
admixture of high-momentum components in the ground state is depicted. The energy domain of these 
high-momentum nucleons is at large missing energies as discussed above. 

This rather complete picture of the properties of a strongly interacting proton in the nucleus signals 
a unique position for the field of nuclear physics. Indeed, unlike other fields with strong correlations 
effects, like particle or condensed matter physics, it is now possible in nuclear physics to state that the 
properties of the constituent protons inside the nucleus are identified experimentally and understood in 
global terms theoretically. For atoms and molecules it is also possible to extract this kind of information 
by employing the corresponding (e,2e) reaction |275t I27tij . This reaction generates the best possible 
information on the properties of individual electrons in these systems. Indeed, it was shown in 1981 that 
one can "measure" the square of the Is wave function of the Hydrogen atom in momentum space |277j . 
Similar results for electron wave functions in the medium have been obtained for a wide range of atoms 
and molecules p,TE\ . Similar wave function results for nuclei are given in Fig.^Jof Sec. El This technique 
may also become successful in identiying the properties of electrons in solids .278) . In nuclear physics 
the improved analysis of two-nucleon knockout reactions is expected to provide detailed information 
about the short-range interaction of nucleons in the nuclear medium. With the additional improvement 
of the analysis of experiments that probe high-momentum nucleons in the nucleus, one may therefore 
look forward to an even deeper understanding of nucleon properties and their interactions in the future. 



7 Summary and Outlook 

Several recent developments have been reported that are of general importance to the understanding 
of nuclei. We single out here the improved understanding of the nuclear-matter saturation problem 
and in particular the essentially complete understanding of the properties of protons in nuclei. These 
achievements have been accomplished by applying the method of self-consistent Green's function. This 
method has been presented in this overview together with a summary of the calculations that have 
been performed in nuclear physics. One of the motivations to employ this method is that the quantities 
that form the ingredients of the calculations have a very clear relationship with experimental data. 
This is true for the one-body propagator which contains information about one-nucleon removal (and 
addition) quantities but also for other propagators discussed in this work. These include the polarization 
propagator which yields information about excited states and corresponding transition probabilities 
but also the pp(hh) propagator which for the hh part can be related to two-nucleon knockout data. 
The appearance of all these quantities as ingredients of self-consistent calculations provides us with 
a multitude of opportunities to compare theory with experimental data. One of the strengths of the 
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method is therefore its abihty to identify relevant physical correlations that need to be included in the 
calculation based on this comparison with experiment. As an example, the fragmentation of the removal 
strength observed in nuclei immediately identifies the need for a dynamic content of the self-energy of 
a nucleon below the Fermi energy. The comparison of the theoretical sp strength distribution obtained 
with the second-order self-energy with experimental data then immediately points to the need for a 
better description of the intermediate states in the self-energy. This in turn leads to the consideration 
of the coupling to collective motion at low energy. The process is illustrated by the development of the 
Faddeev technique to the many-body problem discussed in Sec. El which allows the inclusion of these 
collective features of both ph as well as pp(hh) character. In turn, these microscopic phonons yield 
information that can be directly compared with the excitation spectra in the corresponding nuclei with 
A and A ±2 nucleons. This comparison suggests further improvements of these phonons which in turn 
can be taken into account in describing the sp properties and so on. 

Another important advantage of the SCGF method is associated with its ability to handle the 
relevant physical correlations. As has been discussed throughout this review the treatment of SRC 
can be handled by the summation of ladder diagrams whereas for LRC it appears that the Faddeev 
technique is necessary in finite nuclei. The treatment of SRC is most advanced for nuclear matter where 
the method is applied by several groups but for finite nuclei improvements should be studied in the 
future. The discussion of the results for the saturation properties of nuclear matter with the inclusion 
of SRC as given in Sec. 14. 71 suggests that it may be important to actually establish the validity of one of 
the underlying assumptions of the nuclear-matter problem. This assumption is that the best equation 
of state for a given Hamiltonian will automatically yield information relevant for finite nuclei. This 
assumption has been questioned in this review since it is unclear that long-range pionic modes which 
contribute significantly to nuclear-matter binding actually play a significant role in finite nuclei. We 
therefore urgently suggest investigations which clearly identify whether pion-exchange in ^osp]-, jg j-g^hy 
similar to the corresponding process in nuclear matter. If the answer is negative, nuclear physicists may 
have worried too much in the past about their inability to describe the empirical saturation properties 
of nuclear matter. 

It has been argued in this review that it has been possible to track the properties of all the protons 
in closed- shell nuclei. We emphasize that this achievement reflects the combined effort of experiment 
and theory. The arrival of new facilities which aim to probe the properties of nuclei far off stability will 
also provide us with new challenges to theory. The inclusion of the continuum in SCGF calculations is 
a technical challenge but the method appears sufficiently flexible to include this difficulty. Knowledge 
of the sp properties of these exotic nuclei will widen the scope of our understanding of the properties 
of nucleons in nuclei. Additional work is also necessary to improve the description of the spectrum and 
sp strength in ^^O. It appears likely that the Faddeev method can be applied to study the low-energy 
properties associated with particle addition (elastic scattering) which also has relevance to astrophysical 
problems. The flexibility of this method is not restricted to nuclear physics and may also be applied to 
other many-body systems like the electron gas. 
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